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ABSTRACT 


A large-scale cycle of available potential energy in the Northern Hemisphere over a 


period of about two weeks during late December 1958 and early January 


investigated in some detail. 


1959 has been 


During this cycle the zonal available potential energy first 


built up strongly to a maximum, and then when it began to decline, increases in eddy avail- 
able and eddy kinetic energy took place. These changes in the energy parameters wer¢ 
well related to variations in the poleward heat transport, large values of which signify 


substantial conversions from zonal to eddy available potential energy, and to variations 


in the conversion between potential and kinetic energy. 
the generation of available potential energy 


potential energy variations. 


Furthermore some estimates of 


show good consistency with the available 


Examination of this evecle of available potential energy on a 
d ] } 


regional basis indicates that it was almost completely dominated by developments over 


North America and vicinity. 
also illustrated 


1. INTRODUCTION 


Satellite observations of solar and terrestrial radiation 
will soon be available in sufficient quantity to provide 
the first direct measurements of the global heat budget. 
As radiation is measured on a continuing day-by-day 
basis it will be of great interest to study not only the 
large-scale spatial and temporal variations in this basic 
Planetary heating, but also the manner in which the 
atmosphere reacts to it. Since this heating is indeed 
basic in nature, any clear-cut reactions of the atmosphere 
to it would most likely be found in some of the more 
fundamental atmospheric parameters that directly meas- 
ure the large-scale energetics of the circulation. In 
particular it seems appropriate to investigate the role of 
tadiational heating with respect to the zonal and eddy 
‘components of kinetic and available potential energy as 
defined, for example, by Lorenz [4]. 


These various energy components and the typical flow 


'This 
tration 


search has been supported by the National Aeronautics and Space Adminis- 


The synoptic events associated with this energy cycle are 


the environment and between 


various forms are illustrated schematically in figure 1 in 


of energy from and to 


the manner presented by Lorenz [3]. Here it is seen 
that diabatic heating of the atmosphere (actually differ- 
ential heating between high and low latitudes) builds up 
the zonal available potential energy, which is basically a 
the latitudinal 


essentially a measure of the strength of latitudinal thermal 


function of variance In temperature, or 


eradients. As the zonal available energy keeps increasing, 
the strengthening thermal gradients lead to baroclinic 
instability and the growth of eddies. Thus eddy available 
thermal differences 
the the zonal 
available energy as northward heat transport increases 


potential energy (increasing within 


latitude circles) increases at expense of 


to relieve the excessive thermal gradients. Some of this 
increased eddy potential energy is then converted to eddy 
kinetic energy by means of direct thermal circulations 
(upward motion in the warm air, downward motion in 
the cold air) occurring in the eddies, or major cyclone 
waves. Differential heating relative to these eddies along 
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Proposed energy cycle for the earth’s atmosphere (after 
Arrows indicate prevailing directions of energy 
Dashed lines indicate that there is apparently 


Ficure 1. 
Lorenz [3]). 
transformations. 
no prevailing direction. 


latitude circles may either generate or dissipate the eddy 
available energy—a predominant trend in one direction 
or the other has not as yet been definitely established. 
Eddy kinetic energy is partially dissipated by friction 
and partially converted to zonal kinetic energy (much of 
the conversion being associated with the poleward transfer 
of momentum). The zonal kinetic energy is mainly sub- 
ject to dissipation by friction, but a portion of it may be 
converted back into zonal available potential energy. 

This typical energy flow picture has been generally 
substantiated by empirical evaluations of energy conver- 
sion terms (White and Saltzman [13]; Wiin-Nielsen [14]; 
Saltzman and Fleisher [10]; Jensen [1]). Although this 
mode of energy flux is representative of conditions that 
apparently must prevail on the average, it is realized that 
the generation and transformation of atmospheric energy 
may very well vary in both amount and direction at 
various times. This is the case both in the seasonal sense 
and also over time periods as short as a week or even a day, 
particularly when the energy in only one hemisphere or 
less is considered. 

To investigate more closely the behavior of the energy 
cycle over various time and geographic scales a program 
for regular daily computation of several atmospheric 
energy parameters of the Northern Hemisphere has been 
initiated (Winston [16]). As a substantial sample of such 
calculations accumulate, it will be possible to document the 
Moreover, as radiation data 


energy cycle more precisely. 
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observed by satellites also accumulate, there will be an 
opportunity to calculate the generation of available 
potential energy (primarily the zonal component) and to 
study the reactions of the energy cycle to variations in 
this radiant energy. These are areas of study for the 
future, however. 

For the present it is of interest to illustrate an actual 
case of a pronounced cycle of available potential energy 
which occurred during a pilot period of 42 days of energy 
computations in the winter of 1958-1959 (Winston [16)). 
This case seemed to fit the idealized energy flow picture to 
such a marked degree that it was of interest to examine it 
in some detail. Presentation of this case is of special 
significance because it is, to our knowledge, the first time 
that a cycle of large-scale available potential energy has 
been investigated for the actual atmosphere. In general 
circulation model experiments, on the other hand (Phillips 
[8]; Smagorinsky [12]), there has been extensive ‘“day-by- 
day” calculation of potential and kinetic energy and 
energy transformations. Yet it is surprising that there 
has not been a like body of data computed from the real 
atmosphere to check similarities of the models and the 
observed time variations of the atmospheric energetics. 
Occasionally the energy cycle in the models has been com- 
pared to the index cycle of the atmosphere, which was 
described in detail by Rossby and Willett [9], and indeed 
the energy cycle and index cycle in a gross sense are often 
similar phenomena. However, as some experience has 
been gained with energy parameters, it has been found that 
the index cycle and the energy cycle do not always occur 
at the same time and with the same degree of intensity. 
To a great extent this is due to physically real differences 
in the behavior of the zonal winds as contrasted with the 
energetics of the flow. But also it is in part due to the 
inadequacy of the zonal index as a general circulation 
In fact, if one stresses the broader connotation 
has attained, particularly 


parameter. 
that the term 
in terms of cold air penetrations to lower latitudes and 


“index cycle”’ 
expanded circumpolar westerlies (Namias [5]; Rossby and 
Willett [9]), then it is highly probable that the dynamic 
and thermodynamic consequences of the energy cycle are 
truly more descriptive of these typical “index-cycle” 
phenomena than any index variations themselves. Thus, 
the expedient of verifying energetics of general circula- 
tion models in the real atmosphere by means of the index 
It must be pointed 
out, however, that there have been comparisons of the 


cycle is a patently crude procedure. 


average energetics of models with average values of energy 
transformations computed for the atmosphere. 

The first published calculations of available potential 
energy for the Northern Hemisphere in the actual atmos- 
phere, which were made at about the same time as | hose 
calculated by us, are some average values for the month 
of February 1959 (Saltzman and Fleisher [10]). Their 
results consist not only of monthly average values of 
total, zonal, and eddy available potential energy, but also 
the components of the eddy available energy for indiv dual 
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wave numbers around latitude circles. It is believed that 
the present study of one particular energy cycle provides 
some interesting material which is complementary to that 
of Saltzman and Fleisher and which also aids in the 
physical comprehension of the atmospheric energy cycle. 


2. DEFINITIONS AND METHODS OF COMPUTATION 


The average available potential energy per unit area 
over the entire earth’s atmosphere may be expressed by 
the following approximate equation, as derived by Lorenz 
(4] 

i 1 "Po — Tre—1 77 oe 
A=35 (T,—T)-'T-(T")*dp, (1) 


where bars within the integral refer to averages over an 
entire isobaric surface, 7 is temperature, Tg is the dry 
adiabatic lapse rate, T is the lapse rate, 7” is the devia- 
tion of the temperature at any given point on the isobaric 
surface from its average value in that surface, p is pres- 
sure, and po is the average pressure at the surface of the 
earth. It should be noted that over short time periods 
i.e., of the order of days, weeks, or perhaps a month) 
variations of available potential energy are essentially 
dependent upon variations in the spatial variance of 
temperature on isobaric surfaces, (7”)*, since the mean 
temperature and mean lapse rates within an isobaric sur- 
face are relatively constant with time. In a longer-period 
sense (e.g., over several months) the time variations in 
available potential energy are, however, also influenced 
to a moderate degree by changes in mean lapse rate and 
to a lesser extent by changes in mean temperature. 

Available potential energy may also be expressed in 
terms of two components, a zonal and an eddy component 
Lorenz [4]). The expression for zonal available potential 
energy Is 

A,=* (" (r.—-T)—'T-"[T yp, (2) 

where the bracket refers to a latitudinal average and 
T}' is the departure of the latitudinal average temperature 
from the overall mean temperature. The expression for 
eddy available potential energy is 





— 1 Po — 


Ag (T4- I) iT \(T*)*dp, f2 


~~ 


where 7* is the departure of the temperature at any 
point from its latitudinal average value. Thus, zonal 
available potential energy is mainly dependent upon the 
the latitudinal 
essentially upon the strength of the north-south thermal 
gradients, while eddy available potential energy is mainly 
i function of the variance of temperature within latitude 
circles, 


variance of average temperatures, or 


having a minimum value when isotherms are 
parallel to the latitude circles and having maximum 


— —___. 


?Temperature data recorded at 5° latitude and 10° longitude intervals (diamond grid) 
by the tended Forecast Section of the Weather Bureau were used. 
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values when the isotherms display a wave pattern of 
large amplitude. 

In attempting to calculate available potential energy in 
the actual atmosphere from synoptic data it is necessary 
to restrict the computations to a limited area where 
In the caleula- 
tions presented here, the area considered was the Northern 
Hemisphere latitude 30° northward. 
atures ? at 700 mb. were used to represent the temperature 


analyzed temperature data are available. 


from Temper- 
field in the entire troposphere, and the calculations, which 
were made by means of a desk computer, were considered 
to apply to the layer between 1000 and 250 mb. No 
attempt has been made to weight these energy calculations 
for pressures less than 250 mb., since the available poten- 
tial energy would tend to be much smaller in most of the 
more stable stratospheric layers than in tropospheric 
layers of comparable pressure intervals. A few 
comparing available potential energy computed from 
data throughout the troposphere and stratosphere, with 


tests 


values using 700-mb. temperatures only, show that the 
use of 700-mb. temperatures alone provides a good ap- 
proximation to the available potential energy through 
the depth of the atmosphere. Another approximation in 
these calculations is the use of the standard atmosphere 
lapse rate for P. 

The average kinetic energy per unit area over the 
entire atmosphere may be expressed as follows (cf., 
Lorenz [4]): 


"PS «ne 


: V2dp, (4) 


where g is the acceleration of gravity, V is the horizontal 
wind speed, and the bar within the integral again repre- 
sents averaging over an entire isobaric surface. Kinetic 
energy may also be expressed in two components, zonal 
kinetic energy and eddy kinetic energy, the expressions 
for which are respectively as follows: 


= 1 7) ——— 

Kz =9 : \V Pdp, (5) 
and 

7 l "2. cee 

Kr=59"' | (V*)*dp, (6) 


where [V] is the average wind for each latitude circle and 
V* is the deviation of the wind from its latitudinal average 
value. Kinetic energy and its components were computed 
by means of a machine program which had previously 
been devised by C. L. Bristor. 
of 500-mb. geostrophic winds obtained from 500-mb. 
height data of the Joint Numerical Weather Prediction 
(JNWP) Unit on their octagonal grid, which fully covers 
the area between about latitude 20° N. and the pole. 
Expressions for time rates of changes of available po- 
tential and kinetic energy are also pertinent. Following 
Lorenz [4] we may write the following expressions for the 


rhis program makes use 


local changes of energy: 
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94 =—C+6 (7) 
ol 
and 
os dD, (S) 
ot 
where s 
C=—Rq"' |p 'Twdp, (9) 
G g 1 “rr, T) iT 'T’ (dp, (10) 
and 
= 
D=—g"'} — V-Fdp. (11) 


In these expressions (’ represents the conversion between 
potential and kinetic energy, G@ is the generation of avail- 
able potential energy, D is the frictional dissipation of 
kinetic energy, /? is the gas constant for dry air, w 1s equal 
to dp/dt (vertical motion in pressure coordinates), Q is the 
rate of addition of heat per unit mass, and F is the fric- 
tional force. 

Generation of available potential energy is greatest 
when there is a high positive covariance between heating 
and temperature (equation (10)). In other words, when 
the atmosphere is heated where the temperature is high 
and cooled where the temperature is low there is consider- 
able generation of available potential energy. The energ\ 
conversion term appears with opposite sign in equations 
(7) and (8). Energy is converted from available potential 
energy to kinetic energy when there is a correlation be- 
tween vertical motion and temperature such that upward 
the warm air and downward motion 
On the other hand, 


when the vertical motion field is reversed relative to the 


motion occurs in 


occurs in the cold air (equation (9)). 
temperature field, i.e., downward motion in the warm air 
and upward motion in the cold air, the kinetic energy is 
converted to available potential energy. 

It is obvious from these expressions that the place where 
radiational heating plays an important role is in the gener- 
This 


will be one of the most important terms to evaluate when 


ation of available potential energ\ (equation (10) 


satellite radiation data become available, and of course, it 
will be important to compare these evaluations of the 
veneration term with the observed variations in available 
potential energy. At present, however, some calculations 
of the energy generation term can be made from indirect 
estimates of heating on a hemispheric scale (Wiin-Nielsen 
and Brown [15]). 

The conversion term (equation(9)) is also of consider- 
able interest in understanding the time variations of both 
kinetic and available potential energy. From data kindly 
provided by A. Wiin-Nielsen, some evaluations of this 
term have been made. These data are based on initial 
vertical motions and thicknesses available as a by-product 
of the operational baroclinic numerical prediction calcu- 


lations of the JNWP Unit. 


These data were used by both 
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Wiin-Nielsen [14] and Saltzman and Fleisher [10] to obtain 
calculations of this energy conversion term. 

Expressions for conversion and generation of zonal and 
eddy components of available potential and kinetic energy 
are also given by Lorenz [4], but it is beyond the scope of 
this paper to present all of them here. It is pertinent, 
however, to write the expression for the conversion between 
zonal and eddy available potential energy, Cy, which is as 
follows: 





C;, RR i P T** 


|7'|’ 
dJo 7 0 ) dp, 


(12) 


denotes a deviation 


pera os 
gy titel ap) r, 4 


where @ is potential temperature, * 


from a latitudinal average value, and [ |’ denotes a devi- 
ation of a latitudinal average value from the mean value 
that 
basically this conversion involves the products of the 


over the entire isobaric surface. It will be noted 
horizontal transport of sensible heat and the north-south 
temperature gradient and of the vertical transport of heat 
and the vertical temperature gradient. 

In the discussion that follows only the northward heat 
transport itself is calculated. This transport, /7, across a 
given latitude circle, @, in a laver between two constant 
pressure surfaces, p,; and po, is obtained through the follow- 
ing approximate expression : 


Py ) AX - \ 


— 


¢,a COS Q, ( Pr 
H ; 


v2, 13 


where a is the radius of the earth, ¢, is the specific heat at 
constant pressure, v is the meridional wind component, and 
\ is longitude. In the evaluation of (13), p. is taken as 
850 mb. and p, as 500 mb., e* and 7* are average values 
in the laver 850-500 mb., and the summation is over 72 
intervals of AX, or 5 
Further treatment of the assumptions used in this 


of longitude, around the latitude 
circle. 
calculation, including the use of geostrophic flow, has 
already been given elsewhere (Winston [16]). 


3. DESCRIPTION OF THE ENERGY CYCLE AND VARIA- 
TIONS IN PARAMETERS ASSOCIATED WITH ENERGY 
CONVERSIONS 


The day-to-day variations of the available potential 
energy, and its zonal and eddy components for the period 


of interest are shown in figure 2. It 


may be seen that 
the major cycle of buildup and breakdown in the total 
available potential energ\ between December 24 and 
January 7 mainly reflected the variation in the zonal com- 
ponent of the available potential energy. This is essen- 
tially indicative of the buildup of strong thermal gradients 
between latitude circles up toa peak on January 1, whieh 
was then followed by a very substantial breakdown in 
these thermal gradients between January 1 and 10. The 


eddy available potential energy remained at a 1 arly 








(10'S ERG 


con: 
ren 
whe 
pre’ 
able 
on - 
only 


19 


Lem, 
‘| 


Con 


\ 


stro 
Was 
wor, 


] 
Cire} 








SEPTEMBER 1961 


60 ) .60 







mf 56 
TOTAL AVAILABLE 
POTENTIAL ENERGY 

of y 52 
48 |} : 
7 44 
40 | . 

) ZONAL AVAILABLE A : 

POTENTIAL ENERGY yo \ 
4 >) 
, a 
36 \/ \ i 
‘ \ 
; \ 

. ies ' \ 25 

: — \ F ‘ . 

\ sn f ‘ 
Ve 

_ \ 22 
Cc . 


Q OTAL KINETIC ENERGY SN ( 
24} 4 U 24 
VK ‘ ‘ 
/ 


X 
pd ; 
PO . eX 
\ ee v\ 
\ 





a / anh 
20 ] \ / cv 
EDDY AVA A \ \ 
Pi 22 didlos ie 
AL / v \ 
7 a 
16 \ ae EDDY KINETIC ENERGY ’ 
\ ani y Nes 
12 | wc Nw ho ie a 
ie / Oe te, \ ~ 
~ y ay” - a 4 See 
- a Y= ZONAL KINETIC ENERGY 
08 08 
22 24 26 28 30 3-5§ Ff oH 8 
DEC. 1958 JAN 1959 


Figure 2 Available potential energy, kinetic energy, and the 
zonal and eddy components of each for the period December 21, 
1958-January 13, 1959. 
OOOO GMT. 


All quantities computed once daily at 


constant, low level until December 29 and then began a 
general rise culminating in a peak value on January 4, 
when the zonal available potential energy was falling 
precipitously. Following this maximum the eddy avail- 
able energy dropped at a rapid rate to reach a minimum 
on January 7, but a rapid rebound brought it to values 
only slightly less than peak values between January 8 and 
12, when the zonal component reached its minimum value. 

This overall behavior of the eddy relative to the zonal 
component of the available potential energy suggests very 
strongly that the increased level of eddy available energy 
was derived from the zonal available energy. In other 
words, strong thermal gradients developed within latitude 
areles at the expense of the thermal gradients between 
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Figure 3.—(a) Total heat transport across all latitudes from 
20° N. northward, and (b) the conversion between availablk 


potential and kinetic energy (positive for conversion to kinetic 
period December 21, 1958-January 13, 1959 
Portion of 


energy), for the 
Both quantities computed once daily at 0000 Gaur 
energy conversion curve prior to December 26 is dashed since 


conversion data were unavailable for December 23 and 25 


latitude circles. This is supported by the variation of 
the total northward heat 
latitude circles) which is shown in figure 3. 


(summed over all 
Note that 


heat transport was at a relatively low level during the 


transport 


buildup of zonal available potential energy, but spurted 
rapidly upward at the time when the zonal available 
energy began to decline after January 1. In view of the 
direct role that heat transport plays in the conversion 
from zonal to eddy available potential energy (equation 
(12)) there is little doubt that this increase in heat trans- 
port after January 1 signified a marked increase in this 
energy conversion. 

Turning to the kinetic energy variations in figure 2, it 
is seen that the total kinetic energy also went through a 
cycle at this time, but it was almost completely dominated 


Note that 


the eddy (and total) kinetic energy generally increased 


by the variations in the eddy kinetic energy. 


for more than a week before reaching a maximum between 
about January 6 and 7, more than 2 days later than the 
Thus the 
next step in the energy flow, the transformation from 


peak in the eddy available potential energy. 
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eddy available to eddy kinetic energy (see fig. 1) evidently 
occurred, 

As for the final step in the typical energy flow, the con- 
version from eddy kinetic to zonal kinetic energy, there 
appears to be very little evidence that any appreciable 
increase in conversion to zonal kinetic energy was realized 
in this period since the zonal kinetic energy showed rather 
small variations. There 
maximum in the zonal kinetic energy on January 6, about 


time was, however, a weak 
the same time as the maximum in the eddy kinetic energy. 

The conversion between total available potential and 
total kinetic energy shown in figure 3 generally verifies 
the increased flow of energy from available potential to 
kinetic during the first several days of January. Com- 
parison of this conversion with the time changes in total 
kinetic energy that can be deduced from figure 2 shows 
that they were generally related in the direction indicated 
by equation (8); i.e., the conversion term tended to be 
larger when kinetic energy was markedly increasing and 
smaller when kinetic energy was markedly decreasing (on 
the assumption that D always acts to decrease kinetic 
energy and that it does not vary nearly so much as the 
conversion term 

It is interesting that the time variation of the conver- 
sion term, which is a function of the spatial covariance 
between vertical motion and temperature, appears to be 
total heat 
transport, which is a function of the spatial covariance 
between meridional This 
tends to verify the deductions of Kuo [2] that northward 


rather well correlated with the variations in 


motion and temperature. 


heat transport (and hence conversion from zonal to eddy 
available potential energy) and conversions from potential 
be closely related. However, 


to kinetic must 


despite this good correspondence in these quantities, the 


energy 


actual curves of eddy available and eddy kinetic energy 
(fig. 2) are imperfectly correlated, particularly near the 
time of maximum values in the two eddy components 
(January 3-9). In fact, a secondary minimum in eddy 
available potential energy occurred on January 7, very 


The 


explanation fer this lack of agreement probably lies in the 


close to the maximum in eddy kinetic energy. 


behavior of the generation of eddy available potential 
energy which apparently was responsible for the downturn 


> 


in eddy available energy between January 3 and 7. Some 


evidence for this will be presented in the following sections. 


4. THE GENERATION OF AVAILABLE 
POTENTIAL ENERGY 


It was indicated earlier that for a closed region the 
available potential energy can change due to either a 
conversion term, depending on a correlation between 
vertical motion and temperature, or a generation term, 
depending on a correlation between diabatic heating and 
temperature (equations (7), (9), (10)). The variations 


of the conversion term have already been shown (fig. 3) and 


Although readily computed through 


discussed above. 
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1958—January 


Quantities computed once daily at 0000 Gar. 


and eddy components for the period December 27, 


13, 1959. 


use of numerically calculated vertical motions its accuracy 
is not readily determinable. The problem of obtaining 
the diabatic heating, and consequently reliable estimates 
of the generation term however, is an even more difficult 
Recently Wiin-Nielsen Brown 
[15] have devised a method for computing the diabatic 
Their method 


problem. however, and 
heating hemispherically on a daily basis. 
essentially utilizes an equation derived from the thermo- 
dynamic energy equation where the vertical motion is 
estimated from the vorticity equation. 

By utilizing their method it is now possible to estimate 
the generation of available potential energy from synoptic 
data, and indeed this was also done by Wiin-Nielsen and 
Brown. Their computations of the average generation 
term for the month of January 1959 indicate a value of 
about 1.5 10° ergs em.~* sec.~', a value comparable to 
the average conversion from potential to kinetic energ) 
for the same month (Wiin-Nielsen [14]). This represents 
10° ergs em.~? sec.~' for the zonal vom- 
and —3.5> for the eddy 
Thus according to these calculations, the 


the sum of 5.0> 
ponent 10° ergs cm.~* sec.~! 
component. 


only source of available potential energy, at least diring 
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Figure 5.—Observed and computed available potential energy for 


the period December 27, 1958-January 13, 1959. 


latitudinal heating gradient. 
within latitude 


represents the heating associated with the long planetary 


1959, was the 


Apparently the heating circles, which 
waves and transient disturbances, served only to reduce 
the available potential energy. 

As part of the study of the cycle in available potential 
energy treated in this report, it was felt to be of interest 
lo compare the interdiurnal variation of the generation 
term and its zonal and eddy components with that of the 
potential For this Wiin- 
Nielsen has kindly made his computations for January 


available energy. purpose 
1959 available to us and these have been combined with 
similar computations made in the Meteorological Satellite 
Laboratory for the period December 27-31, 1958. 

The time variations of these energy generation terms 
are shown in figure 4. As can be seen, the zonal genera- 
tion was always positive, seldom dropping below 3 units, 
while the eddy generation was predominantly negative, 
iveraging frequently below —2 units. The total genera- 
tion, while usually positive, was occasionally negative, 
particularly between January 4 and 7, indicating that the 
heating field even in winter can at times act to destroy 
ailable potential energy. Of particular interest for 
this study was the similar behavior of the generation 
term and the total available potential energy (fig. 2) 


Both 


increased very markedly at the beginning of the period 


between December 27 and January 7. curves 


with the generation term reaching a maximum of 7 x 10° 
— about two days prior to the maximum 
At this time the eddy 
generation was small so that the total was nearly equal 


t 


0 the 


_ | 
ergs em.~? see. 
in the available potential energy. 
zonal 


reflecting a predominantly 


latitudinal heating distribution. 


component, 
Also significant is the 
lact that the generation term remained above 2 units for 
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Figure 6.—Available potential energy in four zones extending 
northward from 30° N. within the indicated longitudinal bound- 
aries for the period December 21, 1958-January 13, 1959. Quan- 


tities computed once daily at 0000 Gur. 


about six consecutive days. Considering the magnitude 
of the energy conversion term (fig. 3) this was sufficient 
to cause an increase in available potential energy during 
this period. 

After January 1 the intensity of the total generation 
of available potential energy was drastically reduced and 


Much of 


this decline was due to the sudden drop in the eddy gen- 


within a few days it had even become negative. 
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January 1, 1959, (c) January 4, 1959, and (d) January 10, 1959. 


Contour interval is 400 ft. 


eration term which fell nearly 6 units during the same 
interval. This was the result of more intense cold air 
outbreaks and the consequent heating of polar air masses 
in middle and lower latitudes, as will be illustrated in the 
next section. The sudden dip in the eddy available po- 
tential energy between January 4 and7 (fig. 2), which was 
discussed earlier, was probably a response to this strong 
heating and a direct consequence of the large negative 
value of the eddy generation. This was a period during 
which the computed values of both the (total) generation 
and (total) conversion terms contributed toward a decline 
in total available potential energy, agreeing with the ob- 
served decline. It is of interest in this regard to compare 
in more detail the computed behavior of the sum of C, the 
conversion term, and G the generation term, with the 


observed variation in available potential energy during 
this period. This will now be considered. 


Equation (7) may be integrated to give 


A, A, i+(G—C)ét, 14) 


is the value of available potential energy on a 


given day n, A,_; is the value a day earlier, and 6f is an 


where A, 


interval of 1 day. Specifying the value of Ap, which has 
been chosen here to be the observed value of the avail ible 
potential energy on December 26, one can compute A; 
and all succeeding values of A, for as many days as one 
pleases, providing values of G—C are available. These 
computed values of available potential energy are shown 
in comparison with the observed values for our perio: of 
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Ficure 8.—700-mb. isotherms for (a) December 26, 1958, (b) January 1, 1959, (c) January 4, 1959, and (d) January 10, 1959. Interval 


is 5° C. with area below —25° C. stippled. 


interest in figure 5. It is quite clear that the increase in conversion computed by this method should be increased 
available potential energy between December 27 and by about 40 percent or more. The weak values of the 
January 1 is explained very well by the sum of the genera- energy conversion probably also reflect the lack of in- 
tion and conversion terms, although the computed values clusion of moist adiabatic processes in the vertical motion 
do not rise to as high a maximum as the observed. The computations. Smagorinsky [11] has indicated that when 
most serious discrepancy between computed and observed these are included the upward motions are found to be 
values occurred after January 3 during the period of much larger than those obtained by assuming dry adiabatic 
strong decline in available potential energy. This motions. Since moist adiabatic effects would be expected 
period was represented by only a weak decline in the com- to be larger in warm air masses, the conversion from po- 
puted values, which suggests that the computed conver- tential to kinetic energy would also be expected to be 
sion term particularly after January 1, while indicating — larger. 

the correct trend, was too small. This is most likely Despite these discrepancies, the encouraging fact is that 
attributable to the fact that the vertical velocities used in such realistic estimates of the variations in available po- 
computing C are too small. Similar conclusions were — tential energy were obtained from the generation and con- 
were reached by Palmén [7], who suggested that the energy version terms, particularly since these terms were com- 
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puted on the basis of rather restrictive and incomplete 


models of the atmosphere. 


5. REGIONAL AND SYNOPTIC FEATURES 
OF THE ENERGY CYCLE 


The variations of total available potential energy shown 
in figure 2 were also examined from another point of view 
(as contrasted with the zonal and eddy components pre- 
viously considered). The available potential energy 
was partitioned into four regional components representa- 
tive of sectors 90° of longitude wide, stretching from the 
pole to 30° N. 
arbitrarily, but roughly in such a fashion as to cover the 


The longitudinal bound- 


These four sectors were chosen somewhat 


two oceans and two continents. 
aries of each zone and the portion of the available po- 
tential energy contributed by each of these zones are shown 
in figure 6. The most striking feature of the variations of 
these regional components was that the major cycle of 
buildup and decline in available potential energy showed 
up only in the component of energy in Zone 1 (North 
America and vicinity). Thus the variations in the avail- 
able potential energy for the whole Northern Hemisphere 
were basically dependent on the events in only this one 
particular region. The dominance of one area over 
another is not too unusual, of course, since differences in 
the types of flow patterns over various parts of the 
Northern Hemisphere have frequently been noted (ef., 
Namias [6]). 
of an energy cycle at this time in three of the four zones is 


However the virtually complete absence 


quite surprising, particularly since the total energy over 


all longitudes showed such a large variation. 

Further insight into this energy cycle is attained 
through inspection of synoptic charts, which showed 
some interesting developments over North America and 
also over other portions of the hemisphere during this 
period. To illustrate these events, 500-mb. height and 
700-mb. temperature fields for four selected days are 
presented in figures 7 and 8. 

Confining our attention at first mainly to the North 
American area, we see that on December 26, when avail- 
able potential energy was at a relatively low value, both 
North America (Zone 1 in fig. 6) and over the 
there was a rather ill- 


over 
hemisphere as a whole (fig. 2), 
defined circulation pattern (fig. 7a) with weak thermal 
gradients (fig. 8a). The high values of available potential 
energy which prevailed by January 1 were characterized 
by a deep cyclonic vortex over the Canadian Arctic 
(fig. 7b) with an associated extensive pool of cold air (fig. 
8b). By January 4 a cyclonic vortex was located in 
western Canada, some 20° latitude south of the one in 
figure 7b, and the broad cyclonic flow to its south domi- 
nated most of the United States (fig. 7¢). Accompanying 
these circulation developments, an extensive tongue of 
cold air covered western Canada and was also invading 
the United States (fig. 8c). This date marked the 
beginning of the overall drop in available potential energy 
6) and also the time of the 


over North America (fig. 
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maximum in eddy available potential energy for the 
hemisphere and of a very rapid drop in the zonal available 
energy (fig. 2). Thus the southward penetration of cold 
air over North America was one of the major instru- 
mentalities in the breakdown of the strong supply of 
available potential energy that had built up over North 
America and the hemisphere as a whole. Undoubtedly a 
major part of its contribution to the energy breakdown 
was in its role in the conversion of zonal available into 
eddy potential energy and thence into eddy kinetic energy, 
all of which occurred at this time (fig. 3). 

By the final day in this sequence, January 10, when 
available potential energy over North America (fig. 6) and 
the zonal available energy for the hemisphere (fig. 2) had 
dropped to relatively low values again, this cyclonic 
vortex was located near the Gulf of St. Lawrence with its 
main trough extending down the east coast of the United 
States (fig. 7d). The cold air associated with this system 
was now definitely weaker than before; in fact, most of 
the colder air (e.g., temperatures below —25° C.) over 
North America had disappeared (fig. 8d 
reflective of two processes, one the 


This warming 
was undoubtedly 
descent of the cold air which contributed to conversion 
from potential to kinetic energy and the other the heating 
of the cold air mass as it proceeded off the east coast out 
over the Atlantic. Very large amounts of heating were 
computed by Wiin-Nielsen and Brown [15] on the days 
when the cold air moved off the coast (January 5—6) and 
it is evident that the strong heating of this cold air con- 
tributed to the strong negative values of computed eddy 
generation on these days (fig. 4). 

Turning to other features of the flow and thermal fields 
it is of interest to note that a considerable amount of cold 
air was present over eastern Asia on each of the days 
shown in figure 8. This was very likely associated with 
the relatively high level of available potential energy in 
Zone 2 between December 26 and January 10 (fig. 6). 
Also, while a portion of this cold air was located farther 
west in Siberia (figs. 8a, 8b), the level of available energy 
in Zone 3 was fairly high, but as most of it moved east of 
Lake Baikal (figs. 8c, 8d) a noticeable drop in the energy 
in Zone 3 occurred. This more narrow longitudinal 
organization of the cold air near the Asian coast was 
accompanied by the sharpening of the Asian coastal 
Although the 
total energy for Zone 2 remained virtually the same, 


trough subsequent to January 1 (fig. 7). 


computations of longitudinal contributions to the north- 
ward heat transport (Winston [16]) showed that sub- 
stantial heat transports were effected by the Asian coastal 
trough subsequent to January 1. This signified that 
increased transformations from zonal to eddy available 
energy occurred in this zone during the period of develop- 
ment of increased eddy energy in January (fig. 2). Thus, de- 
spite the lack of a cycle in the overall energy, it appears that 
eddy Zone 2 (eastern Asia to the east central Pacific) did 
contribute toward the increase or maintenance of eddy 


energy. 
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January 10, 


Finally it is worth noting that the overall hemispheric 
pattern at 
10, fig. 7d) consisted of four major, large-amplitude waves 
middle latitudes. As in- 
this was a day of minimum zonal 
but with fairly large eddy 
This pattern is to be con- 
(although mainly 
zonal) flow pattern of December 26 (fig. 7a), 


flow the end of the energy cycle (January 


with large cyclone centers in 
dicated in figure 2, 
available potential energy, 
available potential energy. 
trasted with (1) the more heterogeneous 
before zonal 
available energy started increasing; (2) the zonal flow at 
higher latitudes surmounting large-amplitude waves at 
lower latitudes of January 1 (fig. 7b), when zonal available 
energy was at its peak and eddy available energy was 


> 


building up; and (3) the pattern of increasing wave 
auplitudes on January 4 (fig. 7c) when zonal available 
energy was dropping rapidly and eddy available energy 
Was at a peak. 

The latitudinal profiles of zonal wind and temperature 
for the whole hemisphere for three of these days sum- 
the differences at these different stages 

Note that at the beginning (December 
10) of the cycle there were rather 


mMarize some of 
further (fig. 9). 
26) and end (January 
broad zones of westerlies with maximum speeds located 
near latitude 35°. The flow at the end of the cycle was 
generally weaker than at the beginning except south 
of 30°. 
were concentrated in a narrow zone in middle latitudes 
peak near 50°) with a definite minimum in the flow near 
r 25°. These profiles 


By way of contrast the westerlies on January 1 


55° and a secondary maximum nez 
are roughly characteristic of the zonal wind distributions 
at various phases of the index cycle (cf., Namias [5}), 
but it is worth pointing out that the profile for January 1 


1959 (dotted 3 


does not represent a case of extremely high zonal index 
(in latitudes 35°-55° N.) 
The thermal profiles show definitely that the 


was marked by 


rather 
peak zonal available energy (January 1) 
rather low temperatures north of 55°, a cooling of up to 
5° © December 26. This high-latitude 
during a period when the westerlies were strengthening 
near latitudes 45°-55° 
development which Namias [5] described as containment 


from cooling 


was essentially the same type of 


of cold air in the polar cap during periods of strong 
westerly flow. By January 10, 
energy was at its lowest, temperatures were considerably 
higher north of 50° 
latitudes in connection 


when the zonal available 


The outflow of cold air into middle 
and the 
of zonal available into eddy available energy is strikingly 
clear. Or putting it the effects of the 
increased heat transports subsequent to January 1 were 
quite marked in the reduction of the northward thermal 
gradient. 


lower with conversion 


another way, 


6. SUMMARY 
This observed cycle of available potential energy for 
the Northern Hemisphere possessed many characteristics 
which had been deduced from previous evaluations of 
some long-period averages of energy and 
The buildup in zonal available 


parameters 
energy conversion terms. 
energy to a high value and then its subsequent breakdown 
with the accompanying conversion into eddy available 
potential and eddy kinetic energy have all been clearly 
demonstrated. There was, however, little if any reaction 
to this potential energy cycle in the zonal component of 
the kinetic The fact that 


energy. estimates of energy 
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conversion and energy generation obtained from ad- 
mittedly crude atmospheric models showed rather good 
correspondence with the observed short-period Variations 
in the various energy parameters was especially interest - 
ing. The partition of available potential energy on a 
regional (longitudinal) basis showed that this clear-cut 
hemispheric available potential was 


evele in energy 


dominated to a remarkable extent by events in only 


about one-quarter of the hemisphere. Investigations of 
the synoptic aspects of this energy cycle showed that this 
dominance by one region (North America and vicinity) 
was associated with the development of cold air over the 
Canadian Arctic then the subsequent southward 
penetration of this large mass of cold air into the United 


and 
States and out over the Atlantic. Obviously this signifies 
that studies of atmospheric energetics require close atten- 
tion to geographical and synoptic-seale influences. 
Although this evele displayed many characteristics of 
the idealized energy cycle, it should not be concluded that 
this observed eyvyecle is truly typical of most large-scale 
energy eyeles. Cursory inspection of an accumulating 
record of available potential energy over more than a 
vear indeed reveals that large-scale variations of the zonal 
and eddy components of the potential energy occur in a 
variety of ways. Study of the nature of other large-scale 
changes in the energy is contemplated, particularly with 
respect to probable differences in the manner in which the 


potential energy is generated. 
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A tabulation of percent frequencies of inversion and/or isothermal conditions, based 
below 500 feet above station elevation, for Weather Bureau radiosonde stations in the 
contiguous United States for four observation times, provides a sampling of daytime and 
nighttime stability conditions for all time zones. These data are analyzed with respect to 


nighttime surface wind speed and cloud cover. 


Radiosonde data are compared to vertical 


temperature gradient data obtained from meteorologically instrumented towers and valley- 
ridge stations. Finally, an attempt is made to delineate, geographically and climatologically, 


the percent frequency of low-level stability for the entire country. 


1. INTRODUCTION 


Air pollution surveys of entire industrial regions and 
metropolitan areas are in demand by public and private 
agencies for purposes of evaluating levels of air contami- 
nation, air zoning regulations, and identification of 
obnoxious sources. In recent years meteorology has as- 
sumed an important role in the atomic energy industry for 
engineering and operational applications common to the 
industry generally, but in particular because of its appli- 
cation to the study of dispersion of airborne radioactivity 
by the atmosphere. Weather may play an important role 
in either conventional or nuclear industries and the proper 
appreciation of certain meteorological factors is paramount 
for an intelligent analysis and evaluation of air pollution 
phenomena. 

In essence the dilution efficiency of the atmosphere 
depends on the wind and temperature gradients, both of 
which vary vertically, horizontally, and with time. Since 
the general spatial problems of slow atmospheric dispersion 
which affect large areas unquestionably arise at times when 
a stable stratified laver of air exists near the surface, a 
knowledge of the frequency of low-level stability for geo- 
graphical and/or climatological areas will be helpful for 
assessing the potential for air pollution inherent to a given 
region. This study involves the compilation of inversion 
frequencies and the analysis of related meteorological 
parameters in an attempt to delineate, geographically and 
frequency of low-level 


climatologically, the percent 


stability for the contiguous United States. 


2. INVERSION DATA 


The initial step for obtaining statistics of inversion 
“Support for this study was jointly provided by the U.S. Atomic Energy Commission 
U.S. Public Heaith Service 


occurrence was to select upper air data that most nearly 
Table 1 


gives the observation times for radiosonde soundings from 


corresponds to daytime and nighttime periods. 


which vertical temperature gradient data could be ob- 
tained. It was decided that a 2-vear sampling of both 
0300-1500 emt and 0000-1200 GmT observation periods 
would be sufficient to give representative frequencies of 
daytime and nighttime inversion occurrence for each time 
zone. Accordingly, the National Weather Records Center 
at Asheville, N.C., tabulated for selected stations a percent 
frequency of inversions and/or isothermal layers based at 
This level 


was arbitrarily selected to define a minimum vertical depth 


or below 500 feet above the station elevation. 


of air, and it also permitted convenient tabulation of data 
from radiosonde sounding plots. These frequencies, given 
to the nearest whole percent, are based on all observations 
actually taken for a given time and season (Winter; 
December, January, February; Mareh, April, 
May; etc.). These data are listed in the Appendix. 


Spring: 


3. PROCEDURE 


Rather than knowing only a percent frequency of 
inversions for any given observation time, it is desirable 
to have an estimate of inversion occurrence as a percent 


of total time sampled (seasonal or annual Continuous 


TARLE 1 Radiosonde observation times 
I tandard 
Date Green 
wicl 
Eastern | Central Mountair Pacific 
2300 1 Six 1700) lun iL 
Prior to Apr. 1, 1948 | 1100 On00 O500 oso 0300 
ate j 0300 220) 2100 ann luo 
April IMS—June 1957 | 1) ooo cago Oso aO700 
] _ | Oooo 1400 Is) 1700 LU 
une 1957 present | 1a Oo700 0600 O50 400 
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TARLE 2 \pprorimate number of nighttime hours (in 24-hour day) 


Winter, Spring | Summer Fall 
Nortt 15 11 v 13 
Centra 4 ll 10 13 
Sout! 13 12 11 12 


recording of vertical temperature gradient such as that 
obtained from towers gives such statistics, but the radio- 
sonde data available do not permit such a tabulation 
directly. Therefore, to prepare frequencies of low-level 
stability as a percent of total seasonal hours from radio- 
sonde data, the following hypothesis was used for analyz- 
ing data from inland stations: The entire nocturnal period 
was considered stable if an inversion was observed at any 
time during the night, and the maximum percent frequency 
of inversions for any one observation time within a 
nocturnal period was assumed applicable to the entire 
nighttime period; and, conversely, the entire daytime 
period was assumed to be unstable, and low values of 
inversion frequencies occurring at Observation times during 
the daytime heating period were not considered in com- 
putation of inversion frequency as a percent of total 
hours. 

A review of table 1 suggests that for the most part 
daytime observations are made either shortly after sun- 
rise (transition period of stability to instability) or shortly 
before sunset (instability to stability); consequently, the 


‘daytime”’ observations are not truly representative of 


the unstable solar heating period. In general, it is felt 
that ignoring the low frequency values for day time obser- 
vations compensated for the apparent overestimation 
arrived at from attributing the maximum frequency value 
to an entire nocturnal period, 

A rough determination of length of nighttime periods 
was made from the Smithsonian Meteorological Tables 
[9] and resulted in three station categories: (a) the north- 
ern tier of States; (b) an area roughly between 30° and 
40° N. latitude, comprising the largest number of stations; 
and (c) the southern regions along the Gulf of Mexico. 
The number of nighttime hours in each category by 


seasons is shown 1n table 2. 


TABLE 4. 





LS1T 
Gul Vie LS1 
Burrwood, La. (coasta 
Lake Charles, La n'and 
Southeastern Atlantic LS1 
Key West, Fla. (marine 
Miami, Fla, (coastal 
Jacksonville, Fla. (inland 
Northeastern and Middle Atlantic... LST 


Hempstead, L.L., N.Y. 
Lakehurst, N.J. (inland 
\ hington, D.C. (inland 


coastal 


( omparison of coastal and inland inver 
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TaBLe 3.—Example of inversion percent frequency determination for 


inland stations 


Bismarck, N. Dak, 


Maxi 
mum Night Fraction Inversion 
percent Observation time of frequency 
value time (LS1 period 24-hour percent of 
night- hours lay total hours 
time 
i B ixB 
15 
Winter 67 0600 15 4 625 | 67.625=42 
11 . 2 
Spring 61 0600 ll 4 458 61.45 s 
u 
Summer 77 2100 and 0600 i) 4 375 | 77X.375=29 
13 
Fall 69 2100 13 M2 | 69X.542=37 
24 


The method for computing a percent frequency of low- 
level stability for inland stations, defined as a percent of 
total seasonal (annual) hours, is exemplified in table 3, 
Dak. The 


results from this technique were compared to continuously) 


using radiosonde data frem Bismarck, N. 


recorded vertical temperature gradient data from towers 
(discussed later), and appear to be quite adequate for an 
overall climatological analysis. 

In analyzing and computing the inversion frequency 
data, it became apparent that there was a notable excep- 
tion to this technique. Inversion frequencies obtained 
from coastal stations reflect marine influences; that is, in 
coastal areas low-level stability may be either inhibited 
or enhanced by advection, in contrast to inland areas that, 
for the most part, exhibit radiation-type inversions. The 
the Atlantic 
coastal sections of the Southeastern States and the entire 


coastal regions could be divided into (a 


coastal region of the Gulf of Mexico, where relatively) 


warm water exists offshore, and (b) coastal sections 
of New England and the Pacific States 


waters are relatively cold. 


where ocean 
The effects of air-water-land 
temperature differential and associated land-sea breeze 
regimes on low-level stability are apparently reflected in 
the inversion frequency values of many coastal stations 
that are affected by a frequent flow of air from over the 


ocean. In general, for coastal areas along the South- 


sion frequen percent 





Winter Spring Summer Fa 
Bd ) } f 
10 66 77 8&8! 14 23 45 3 40629 21 «1341 2 
1 75 65 70 2 35 74 40 0 13 80 32 6 63 73 
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TABLE 5. 


Representative 


nighttime tion 
value (percent) 
Miami, Fla 1 

iter 39 2200 
ring 27 2200 
immer 27 2200 

I 41 2200 

Oakland, Calif 

Winter 77 0400 
Spring 45 0400 
summer 21 0400 
k 2 0400 

eastern States and Gulf of Mexico, a flow from over 


warm waters tends to inhibit inversion formation at night, 
The 
net result is lower inversion frequency for coastal areas 
than for inland areas in these regions. The cold waters 
off the Atlantic 
produce cool sea breezes which enhance low-level stability, 


while daytime sea breezes exhibit neutral stability. 


Pacific coast and northeastern coast 
resulting in higher frequencies of inversions, particularly 
during daytime hours, for these areas. Examples of these 
coastal regimes are given in table 4, where a list of in- 
version frequencies for coastal and nearby inland stations 
can be compared. 

A review of these comparisons suggests that locations 
only a few miles inland from the coast exhibit a conti- 
Lake 


(e.2.. 


nental-type frequency of low-level stability 
Charles and Jacksonville). Consequently, in analyzing 
the radiosonde data for stations in coastal areas, it was 
necessary subjectively to classify by season each station 
as “inland” or “coastal”? according to the marine effects 
reflected by the frequency data, as well as to consider 
station proximity to the water and prevailing wind direc- 
tions. For stations designated coastal, both daytime and 
nighttime inversion frequency values were considered for 
computing an inversion frequency as a percent of total 
hours, instead of neglecting daytime values, as was done 
for inland stations. An example of attributing both a 
representative daytime and nighttime percent frequency 
value to the respective diurnal periods for computing 
inversion frequencies is shown in table 5 for the coastal 
es of Oakland and Miami. 

lsopleths of inversion frequency, expressed as a percent 


eit 


of total hours, for the contiguous United States are shown 
1 A-E. 


lsopleths of the maximum inversion frequency value 


In figure 


reported at any of the four observation times are shown 
in figure 2 A-E. 
percent frequency of nights (or days) that had at least one 
hou 


these 


Conservatively, these values show the 
of inversion. However, it is reasonable to interpret 
statistics to show a percent frequency of nights (or 
davs, for some coastal areas) having several hours of inver- 
slOlis, since an inversion reported at any one hourly obser- 
vation during the night (or day) implies that such stability 
will continue (or has continued) for more than one hour 


fit already exists. This has been shown by analysis of 


LST observa- 
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Example of inversion percent frequency determination for coastal stations 


Night fraction || Representative | LST observa Day fraction Inversion per- 


of 24-hour day || daytime value tion of 24-hour day || cent frequency 
percent) for total hours 
1 I Ry 14 Rl 
0. 542 6 1000 0. 458 21+3=24 
500 2 1000 . 500 14+1=15 
458 2 1000 542 2+1=13 
500 3 1000 . 500 20. 5+1. 5=22 
2 10 1600 417 15+4=49 
458 14 1600 ‘42 21+8=29 
417 $4 1600 583 90+20= 29 
542 l 1600 458 $4+7=41 
cumulative frequencies of inversions obtained from 


meteorological tower data [1, 4, 8, 10] at inland sites. 


Such data for four sites are shown in figure 3. These 
data are plotted on log-probability scales to enhance the 
“breaks” in the curves (abrupt changes in probability) 


which reflect the length of stable nocturnal periods. 
These analyses of tower data also support the concept 
used to determine inversion frequeney for inland stations 
as a percent of total time. That is, a percent frequency 
of inversions at night for any one hour gives a value that 
is approximately representative for any hour during the 
nocturnal stable period. 

The inherent shortcomings involved in this type of 
analysis are obvious. The subjectivity in the analytical 
methods, differences of topography, geography, observa- 
tion time (relative to the diurnal cycle), and sampling 
periods are sources of possible error. However, a more 


exact quantitative assessment of low-level inversion 
frequency is not easily obtained from existing radiosonde 
data. It is hoped that this study may provide a means for 
comparing one section of the country to another, relative 
to inversion frequency, as well as furnishing information 
helpful for developing a diffusion climatology of the 


United States. 
4. CLOUD COVER AND LOW WIND SPEED 


Meteorologists, from experience, would expect clear 
skies and light surfaces winds at night to favor the forma- 
tion of nocturnal radiation-type inversions over land. 
Cloud cover and wind speed data for available stations, 
tabulated for a 5-vear period in Summary of Hourly 
Observations by the U.S. Weather Bureau were analyzed 
data. These 


with to the inversion frequency 


statistics are tabulated in the Appendix. 


respect 


The inherent subjectivity involved in nighttime obser- 


vations of cloud cover (e.g., thin cirrus), together with 
different exposures and elevations of anemometry, restricts 
any quantitative analysis. However, a statistical study 
was made of the relationship of nighttime cloud cover and 
wind speed to inversion frequency, resulting in poor 
correlation for all seasons. For example, figure 4 indi- 
cates a poor relationship, although a trend is evident, of 
observed maximum inversion frequency to occurrence of 


nighttime winds (<7 m.p.h.) and nighttime cloud cover 
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Figure 3.—Cumulative percent frequency of inversions from meteorological towers and valley-ridge stations: 4) Hanford, Wash., 


(B) Shippingport, Pa., (C) 


(<3/10) respectively, for winter. Although this indicates 
that precise frequency values of low-level inversions may 
not be obtained for a specific location from wind speed 
and cloud cover data alone, isopleths of percent frequen- 
cies of cloud cover and wind speed, shown in figures 5 and 
6 respectively, indicate that the general areas having a 
higher frequency of inversions are also characterized by a 
higher frequency of relatively clear nights with light winds, 
and vice versa. ‘This relationship can be further illus- 
trated by graphic addition of the cloud cover and light 
wind frequency values. The sums of these percent 
frequencies were divided by 2 to give conventional values 
within the range of 0-100. 
(cloud cover+wind speed)/2 are shown in figure 7. 


5. CURRENT OBSERVATION DATA 


Since two of the observation times (0300 and 1500 Gmr) 


These isopleths of nighttime 


used in computing inversion frequency are no longer in 


Middletown, Conn., and (D) Idaho Falls, Idaho 


use, it is desirable to compare the computed frequency 
values to those obtained from a current observation time. 
The 1200 Gut observation value was chosen for compari- 
son, since this is the current local observation time in all 
four time zones that most nearly represents nighttime 
conditions. Such a comparison serves two purposes: 
(1) if percent frequencies from four different local obser- 
vation times compare favorably to the computed frequency, 
this supports the hypothesis used for computing the 
percent values; and (2) additional 1200 GMT inversion 
frequencies, obtained either from towers or other upper 
air facilities, may be analyzed for relatively short sampling 
periods (months) to provide estimates of low-level 
inversion frequency for a given inland site. 

However, since the 1200 Gut observation value was 
used more often than any other value to compute \n 
inversion frequency, expressed as a percent of total tine 
(i.e., of the available radiosonde data the maximum 
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FIGURE 4. 
speed <7 m.p.h. 


percent value was observed most often at 1200 Gm, in 
any season, for at least 60 percent of the stations), it was 
necessary to analyze the unrelated data to justify the 
selection of the 1200 Gar observation value as being repre- 
sentative of the nocturnal period, for all time zones. 
This was done for all inland stations by determining the 
variation of computed frequency, based on the 1200 Gur 
value, from the computed frequency based on the maxi- 
mum observed value, for all cases when the maximum 
value was at an observation time other than 1200 Gmr. 
This comparison showed that for any season the computed 
frequency determined from the 1200 GMT value, when 
averaged for all inland stations, differed by less than 3.4 
percent from the computed frequency value determined 
from maximum frequencies at 0000, 0300, or 1500 Gmr. 
On the other hand, since the fraction (nighttime hours)/24, 
attributed to the length of the nocturnal period is factored 
into the computed frequency value, it can be estimated 
that the 1200 amr frequency value averages about 7 per- 
cent less than the maximum value, when the maximum 
frequeney occurs at 0000, 0300, or 1500 Gur. This indi- 
cates that of the available radiosonde data analyzed, the 
1200 Gur percent value is the most representative of the 
nocturnal period for all four time zones. Also, it was of 
interest to note that the 0300 emt value yields a fairly 
close approximation of nighttime inversion frequency for 
fall and winter seasons; however, the 1500 G@mtT value 
appears to be unrepresentative of the nighttime period, 
except for western sections of the country in winter. 


Relationship of maximum observed nighttime inversion percent frequency to (A 
Occurrence for winter. 
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nighttime cloud cover <3/10 and (B) wind 


The next step was to relate empirically the computed 
inversion frequency to the currently observed 1200 Gur 
percent value. This is done in figure 8, which compares 
the computed percent value on the ordinate (Y) to the 
Each point 
represents a station, and a least square fit provides the 
Coastal stations are not included in the 


1200 Gar percent value on the abscissa (X). 


regression line. 
regression analysis, but they are designated in the graphs. 
Their position relative to the regression line may provide 
an index to continentality. Since the ordinate is deter- 
mined from the abscissa value for about 60 percent of the 
cases, no correlations were attempted. However, despite 
the related data inherent in the statistics, good relevance 
is indicated. In general, the relationships are sufficiently 
good to permit use of the 1200 Gar radiosonde data, for 
inland stations in all time zones, for determining inversion 


frequencies as a percent of total time. 
6. TOWER DATA 


Temperatures obtained from meteorological towers and 
valley-ridge field stations and an instrumented TV tower 
[1, 2, 4, 6, 7, 8, 10, 12] provide inversion frequency statis- 
tics that can be compared to inversion frequencies com- 
puted from nearby radiosonde stations. For all nine such 


stations temperature gradient data are continuously 


recorded and thereby readily provide frequencies in the 
form of percent of total hours. These statistics, which 
can be directly compared to the computed frequencies, 


are shown in table 6. 
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TABLE 6 
a I oul 

Winter 
] \ Nev j i 4 
Yueca Flat, Nev. (v rid 3 
l Or rad nae 3 
I Fermi, Mich, (tow 3 
Mount Clemens, Mich. (radiosonde 15 
Detroit, Mict VV tower 23 
Albany, N.Y. (radiosond 27 
Middletown, Conn. (tower 33 
Nashville, Tenn. (radiosonde 31 
(ireensboro, N.C. (radiosonde 4 
Oak Ridge, Tenn. (towers valley-ridg 41 
Boise, Idaho (radiosonde Al 
Idaho Falls, Idaho (tower 52 
Salem-Portland, Oreg. (radiosonde 32 
Spokane, Wash. (radiosonde 37 
Hanford, Wash. (tower 45 
Omaho, Nebr. (radiosonde 42 
Sioux Falls, 8. Dak. (valley-ridge Pathfinder Site 37 
Pittsburgh, Pa. (radiosonde 24 
Shippingport, Pa. (valley-ridge 14 


Average of Tower Shielding Facility (TSF), Melton Valley (MV), and Oak Ridge National Laboratory site (X-10 Sta 


from valley floor; for MV, SZ =60 ft.; 5 ft. from the floor of a wooded valley; for X-10, AZ 


While most of the towers give vertical temperature 
gradients within a height interval of less than 500 feet, 
this is not a serious obstacle for comparing the two sets 
of data, because nearly all nocturnal inversions (based 
below 500 feet above the surface) reported by inland 
radiosonde stations are based below 100 feet above the 
ground. In fact, most are surface based. The inherent 
discrepancies involved in differences of (a) periods and 
lengths of record, (b) geography, and (c) instrumentation 
are obvious; consequently, it is surprising to find that 
the frequency values colnpare so closely. In general, the 
tower data suggest that the methods of computation for 
determining inversion frequency from radiosonde data are 
sufficiently reliable for general climatological analysis. 


7. OTHER CONSIDERATIONS 


Since most of the Weather Bureau radiosonde stations 
are located at airports which are at the outskirts of cities, 
the inversion frequencies for the most part are represent- 
ative of semi-rural or suburban areas. Consequently, one 
would expect higher frequencies of nocturnal inversions 
in the open country; but, lower frequencies for central 
city areas would be expected because of the mechanical 
turbulence and heat effects imposed by large metropol- 
itan areas [3]. In the Rocky Mountain and Appalachian 
regions most of the weather stations are located in valleys; 
consequently, for these areas the percent values reflect 
inversion occurrence for elevations lower than the average 
for the region. Less frequent inversion occurrence would 
be expected on mountain slopes and_ tops. 

The climatice-svnoptic aspects of the west coast, partic- 
ularly the coastal areas of California and southern Oregon, 
The maritime influences and 


deserve special attention. 
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Comparison of ¢ ompule d inversion freque ney to tower and valley-r idge vertical te m perature gradie nt data 






Percent of total hours 
Period of record Vertical height 
ft 
Spring | Summer Fall Annual 

$Y 37 ) 1 June 55-May 59 
32 3s Ww 43 January 50 March 60 AZ = 500 
2? 30 36 28 «1946-1950 
23 34 30 283) December 56-November 59 AZ = 100-25 
} my 30 23 | February 48-December 50 
as) 36 34 30 | December 56-November 59__| 4Z =300-20 
29 31 $6 31 September 53-August 55 
23 27 $2 29 | June 57—-May 59 

April 56-March 59 AZ =200-5 
27 32 40) 33 | June 55-May 59. 
36 3H) 2 41 | June 55-May 59 
42 4) 4 45 1955-1958 
3s 36 50 44 | June 55-May 59 
37 45 57 is January 56-June 58 AZ =250-5 
29 28 38 32 June 56-May 59 
2 31 40) 34 | June 56-May 59 
$2 27 42 37 1955-1958 = AZ = 400-50 
27 $3 $8 35 | June 56-May 59 
iW) 47 37 43 | May 58-April 60_. AZ =154 
31 27 34 29 | June 56-May 59 
19 23 28 21 | April 55-May 57... AZ =400 


For TSF, AZ =300 ft.; above a wooded ridge top, 500 ft 
135 ft.; 5 ft. atop a smal! ridge overlooking an industrial type valley, 90 ft. below 


a semipermanent high pressure pattern produce dominat- 
ing effects. The large-scale subsidence of air at altitudes 
above the ocean-cooled surface layer results in a stable 
laver which is semipermanent over the California and 
Oregon coasts during spring, summer, and fall months. 
The subsidence inversion acts as a “‘lid’’, limiting the 
vertical motions originating in the marine air below; con- 
sequently, it is an important contributor to the gross ac- 
cumulation of air pollution since the inversion base limits 
the vertical volume of air in which pollutants may 
accumulate. 

Because the subsidence inversion is often based between 
500 and 2000 feet above ground along the coast, it is ap- 
parent that the data used in this study (inversion based 
< 500 feet above ground) vield very conservative fre- 
quency values for the west coast area, relative to any as- 
sessment of atmospheric pollution potential for that area. 
Unlike studies for other areas of the country, a study of 
the diffusion climatology of the coastal areas of Oregon 
and California must consider a frequent and _ persistent 
occurrence of subsidence inversions based below 2000 feet 
above the ground, as well as radiation-type inversions 
based below 500 feet. It is important to note that the 
present study reflects for the most part only surface-based 
radiation-type inversions or low-level stability resulting 
from cool sea breezes. For further details of west coast 
inversions, the reader is referred to studies by Stanford 
Research Institute [11] and J. G. Edinger [5]. 


8. CLIMATIC-GEOGRAPHIC DELINEATION 


A review of the distribution of inversion frequency over 
the country suggests a general delineation into regions »\¢- 
cording to climatological and geographical factors, and 


this classification is shown in figure 9. The factors relit- 
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Ficure 5.—Percent frequency of nighttime cloud cover<3/10: | 
4) Winter, (B) Spring, (C) Summer, (D) Fall, (E) Annual i 
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FiGgureE 6.— Percent frequeney of nighttime wind speed <7 m.p.h.: | 
(A) Winter, (B) Spring, (C) Summer, (D) Fall, (FE) Annual. L 
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Figure 8.—Comparisons of computed inversion frequency 
frequency of inversions on the abscissa (X 


ing to inversion occurrence for each region are discussed 
below. 

Great Lakes. 
storm passages with their associated cloudiness and high 
winds, particularly from late fall to early spring, which 
result in relatively low frequencies of nocturnal radiation 


This area is characterized by frequent 


inversions. An analysis of climatological data for specific 
stations in this area shows a high occurrence of cloudiness 
during nighttime hours, which apparently reflects the ef- 
fects of the Great Lakes imparting moisture to the air in 
i water. Summer and fall 


over the lake 


slightly 


its trajectory 


show low-level 


months higher frequencies of 


stability than do the winter and spring months. 
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Generally, topography plays a small part in this area 
However, the western slopes of the Appalachians probably 
enhance instability in western Pennsylvania and New 
York when the prevailing west and northwest winds blow 
This is particularly true during the winter and spring, 
and is reflected by the frequent occurrence of snow showe! 
activity in these areas. In summary, inversions based 
below 500 feet above the surface may be expected to occu! 
about 20 to 30 percent of the time in any season in thi 
Great Lakes area, although higher frequencies can b 
expected in remote sheltered valley Ss. 

Appalachian Mountains.—Less storm frequency in th 
Southeastern States, compared to the area north o! Vir 
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FIGURE 9 


ginia, undoubtedly contributes to the higher inversion 
frequency found over the southern Appalachian area. 
Mountain ridges afford protection to valley areas; conse- 
quently, light winds during nighttime hours are common 
anticvelonic flow 


Frequent occurrence of 


over the Southeastern States, particularly during summet 


nh allevs. 


and fall months, also enhances nocturnal radiation in- 
version formation, since clear skies and light winds are 
associated with many of these high pressure systems. — In 
general, low-level stabilization occurs 30 to 45 percent of 
the time in this area, while more topographically protected 
areas have more frequent occurrences of inversions. 
Gulf-Atlantie Coast—The warm waters of the Gulf of 
Mexico and Gulf Stream tend to inhibit inversion forma- 
tion along the immediate coast in these regions. During 
the warmer months of the vear the pressure gradient rein- 
breeze over most of the coastal areas of 
This rein- 


lorces the sea 
the Gulf of Mexico and Southeastern States. 
loreement is sufficient to produce relatively strong winds 
during nights along immediate coastal areas, which helps 


0248461 4 


Climatice-geographie delineation of low-level inversion frequency for contiguous United States. 


to delay or inhibit nocturnal radiation inversion forma- 
tion. Extreme southern Florida depicts minimum con- 
tinentalitvy, while stations farther inland (e.g., Jackson- 
ville, Lake Charles, Charleston) show higher occurrences 
In general, from Long Island 
southward Atlantic Hatteras, 
along the Gulf coast and southern Florida, inversions 


Howey er, 


of mghttime stability 
along the coust to and 
may be expected 10 to 35 percent of the time. 
higher frequencies are indicated farther inland, as well as 
for the coastal area between Cape Hatteras and Jackson- 
ville. The higher frequeney of low -level stability along 
the New England coast, reflects the stabilizing effects of 
cold water off shore. 

Central Plains.—The mid-section of the country has a 
pronounced continental type climate, and as such, has 
inversion frequencies closely related to the diurnal cycle; 
that is, there is a definite tendency for nocturnal stabili- 
zation and daytime instability in the lower levels. In 
general, inversions occur 20 to 30 percent of the time 
during the spring and summer, while during the fall and 
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winter months inversions may be expected about 30 to 
45 percent of the time. The higher frequency during 
fall and winter probably is a reflection of minimum 
storminess in fall and maximum length of a stable noc- 
turnal period in winter. The opposite is true for the 
spring and summer months. During the colder months, 
particularly if snow cover exists, warm southerly advec- 
tion over an existing cold surface may enhance low-level 
stability in these areas; however, such advection regimes 
are usually the incipient stages of warm fronts and 
cyclones with their associated precipitation and stormi- 
ness, which ultimately produce less stable lapse rates 
below 500-foot elevations above ground. 

Rocky Mountains. 
nentalitv. With the exception of the northern regions 
during the winter and spring, storms occur infrequently, 
Unlike the 


Appalachian regions, dry air prevails most of the time, 


This area also has marked conti- 


particularly over the southwestern areas. 


and inversion frequency is directly related to the diurnal 
evele. Radiation inversions form near the surface on 
most nights, and they are dissipated shortly after sunrise. 
Since the longest nocturnal stable period and minimum 
daytime solar heating occur during the winter, it is this 
season that is potentially the most stable in this area. In 
general, inversions occur 40 to 55 percent of the time in 
fall and winter, and about 30 to 40 percent of the time in 
spring and summer. 

West Coast. 
including areas only a few miles from the coast, are most 


Along 


the immediate coast radiation inversions may merge with 


Radiation inversions over interior areas, 
prominent during the late fall and winter months. 


subsidence inversions, and low-level stability may persist 
until noon; or, On Occasion, inversions persist for several 
days during which fog often occurs [11]. For the most 
part, the immediate coast has sufficient surface winds to 
provide good ventilation; however, the persistence of the 
subsidence inversion in addition to radiation inversions, 
the daily reversal of wind direction (land-sea breeze), and 
topographical sheltering of the valley and basin areas 
along coastal sections combine to create a potentially 
adverse climate as far as air pollution is concerned. 
Pacific Northwest Coast. 
relatively high winds are dominant features of the climate 
This region is 


Precipitation, cloudiness, and 


of Washington and Oregon coastal areas. 
situated far enough north of the high pressure cell which 
persists off the California coast to permit most storm 
systems from the Pacific to pass inland over the area. 
Consequently, storm activity is frequent, particularly 
during the winter and early spring. The instability 
associated with these frequent storm passages results in 


Radi- 


a relatively low frequency of radiation inversions. 
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ation inversions are most frequent during the summer 
and early fall, when storm activity is at a minimum. 
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APPENDIX: CLIMATOLOGICAL DATA 
Peccent frequency of Station and period of record used 
Inversion Nighttime 
Station Season 
Observations at 0300 and 1500 GMT Observations at 0000 and 1200 GMT 
GMT Total Surface Cloud 
03°15 OO 12 time Wind Cover 
7Tmph <3/10 
on, Ohio Winter $2 22 
Spring 45 29 
Summer si) 52 
Fall . 5Y 45 
Albany, N.Y Winter 47 11 31 438 27 45 28 | Rome, N. Y., September 1953-August Albany, N.Y., June 1957-May 1959 
Spring 63 4 15 45 29 52 29 1955 
Summer 75 #1 27 82 31 65 40) 
Fall “oO 7 47 «66 36 63 4) 
uquerque, N. Mex Winter 67 72 8 82 4s 58 4 Albuquerque, N. Mex June 1955 Albuquerque, N. Mex., June 1957 
Spring 38 24 2 70 32 40 60 ay 1957 May 1959 
Summer 29 «8 4 7l 30 6 4 
Fall 72 4i 2 a 4] 4s 7t 
wn, Pa Winter 49 $6 
Spring 7) 35 
Summer 79 53 
Fall of 4s 
A rillo, Tex Winter 71 68 16 74 43 25 58 Amarillo Air Force Base, Tex., June Amarillo Air Termina lex Ju 
Spring 45 18 4 65 30) 15 ‘4 1955-September 1956; Amarillo Air 1957—-May 1959 
Summer ~~ 610 «3 (RR 34 23 55 Terminal, October 1956-May 1957 
Fall 71 «43 6 70 3s 26 71 
Asheville, N.C Winter 65 42 
Spring LD) 44 
summer oH 49 
Fall s3 
\ nt Gi Winter 65 31 47 738 3s st 4) Atlanta, Ga., June 1953-May 1955... Athe! G June 1957-M 950 
Spring 67 i 2 31 4! 49 
Summer 70 2 13 70 29 65 47 
Fall 70 10 51 «75 $8 A 
Augusta, Ga Winter O68 45 
Spring Os 5 
Summer S2 $s 
Fa sv) 55 
Austil eX_. Winter 45 4) 
43 2 
1 0) ‘ 
I i) 63 
I Md Winter $t 37 
Spring 4 4“) 
Su e1 47 5 
) 3 
I Rou I 4 43 
$ 53 
7] 70 
I Mont 4 20 
7 $5 
35 
$5 14 
I 4 38 
i) 44 
st 4 
hs 7 
I k. N. Dak t 7 37 67 42 14 3 Bis k. N. Dak J r M B k N. Dak J 7-Ma 
i) 4 2 6 2 tal 3 ' 4 
a 69 2 29 2 
6o 4 t 37 414 ul 
I Id 79 7 > s 5} 44 2 Boise, Id 1 55-M VT B I ine 1957-M 159 
t 4 2 84 38 $5 t 
68 33 7 36 j 6s 
” OS & ¥2 Ww te ti 
I \M 7 + 
s 35 
y is 
23 4s 
I e, 1 3 7 it 3 {2 Browns ‘ I June ; M B t I yu 57-M 
21 1 2 22 32 1957 59 
Q » 2 “3 lt 7 
Fa 57 23 ’ 2 3 4 69 
I N.Y Winter 27 18 17 24 17 15 18 | Buffalo. N.Y., June 1955—May 1957 Buffalo, N.Y., June 1957-M 459 . 
Spring SO 44 7 1S s2 
Summer! 7 l 5 SS 26 24 44 
“ 9 » 48 30) ) 39 
k, ¢ f 8] 59 
a3 2 
SS 65 
i) hs 
I m, lowa 2 37 
30) 4! 
Summer 7 57 
ill 3s 5S 
I wd, La Winter 63 36 46 69 40 Burrwood, La., June 1955-May 1957 Burrwood, La,, June 1957-M 1959, 
Spring 57 11 18 51 26 
summer 2H l 4 il 11 
Fall 38 3 «13 33 21 
( Hatteras, N.¢ Winter 0) 63 2 2 $2 Hatteras, N.C., June 1955-February Co Nf Jur 4 M 
Spring 5) 6 43 52 24 1957; Cape Hatteras, March 1957 
Summer 143 11 2 10 May 1957 
Fall 29 5 31 35 19 
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Spring 
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Inversion 


GMT 
15 OO 
27 «39 

7 24 

1 27 
9 49 
23 «5Y 

‘23 

2 ii 
qa os 
38 27 

{ 1 
24 20 
21 WV 
¢ 2 

i 4 
15 31 
7 4 
56 8 
19 
69 47 
s l 
1 2S 
| 3 
33 2 
11 1 
55 4 
91 22 
52 0 
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Percent frequency of 


Tot 
time 


il 


1 
él 


43 


410 
+ at) 


Nighttime 


Surface 
Wind 
7mph 





Cloud 
Cover 
<3/10 


Observations at 0300 ar 


Caribou, Maine 


Charlestor 


( 


olumbia, 


sc 


Mi 


June 


June 


June 


Station and period of record used 


1QAs 


OAS 


M 
May 
M 
M 
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Caribou, Maine, 


Cc 
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j 


Observations at 0000 and 1200 GMT 


June 1957—-May 1959 


harleston, 8.C., June 1957 


May 1959 


( imbia, M June 157-May 1959 
Da 1, O June 1957-May 1959 
Denver, Colo., June 1957-May 1959 
Dodg ( Kar June 1957-M 
1959 

Fl Paso. Tex.. June 1957-May 1959 
Ely, Nev., June 1957—-May 1959. 
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APPENDIX—Continued 





Percent frequency of Station and period of record used 
Inversion Nighttime 
Station Season 
Observations at 0300 and 1500 Gm Observations at 0000 and 1200 cmt 
GMT Total Surface Cloud 
03°15 0O 12 time Wind Cover 
7mph <3/10 
nt, Mich Winter 25 15 28 46 29 Mount Clemens Air Fores Jase, | Flint, Mich., June 1957-May 1959 
Spring 99 9 #4 65 30 Mich., June 1954-May 1956, 
Summer 72 & 2 7 27 
Fall 47 10 21 6&8 34 
t Worth, Tex Winter 58 38 13 55 34 Fort Worth, Tex., June 1955-May Fort Worth, Tex., June 1957-May 1959 
Spring 42 7 2 5i 23 1957 ‘ F 
Summer 45 2 2 8&8 22 
Fall 69 27 18 55 37 
t Wayne, Ind Winter 39 25 
Spring 41 34 
summer 75 53 
Fall 60 w) 
esno, Calif. Winter S4 39 
Spring 64 63 
summer 4s 91 
all sl 73 
veston, Tex Winter 20 Bh 
Spring IS 36 
Summer 22 a 
Fa 31 63 
ilasgow, Mont Winter 75 75 5 76 is Glasgow, Mont., June 1955-May 1957_| Glasgow, Mont., June 1957-May 1959 
Spring 68 38 7 72 33 7 
Summer! ) 14 4 S4 $2 
Fall 79 53 17 8&8 44 
irand Junction, Col Winter 81 86 8 OL 53 Grand Junction, Colo June 1955 Grand Junction, Colo... June 1957-M 
Spring 60 32 0 67 31 May 1957 1959 
Summer 67 17 1 WO 3s 
i ” 7 t 4 44 
irand Rapids, Mi Winter 413 
Spring 44 
Sumime 76 51 
ill 57 45 
i t Falls, Mont Winter 73 74 41 «76 4s 18 $1 Great Falls, Mont June 1955-May Great Falls, Mont June 1957-M 
Spring 32. Sh 1 6S $1 6 32 1957 1059 
Summer 34 2S 27 2s $2 51 
Fall i 70 & 70 38 27 15 
i n Bay, Wis Winter 431 2 OS 34 Green Bay, Wis., June 1955-May 1957 Green Bay , Wis., June 1957-May 1959 
Spring 3 8 9 71 $3 
Summer 766=«C«&@ .. 29 
Fall 59 16 22 61 33 
i nsboro, N.C Winter 73 #15 58 72 43 54 43 | Greensboro, N.C June 1955-May Greensboro, N.C., June 1957—-May 
Spring 70 3°13 66 32 a 414 1957 1959 
Summer 78 1 ll 66 33 79 is 
Fall 74 4 52 74 40 oR 5S 
H sburg, Pa Winter 60 32 
Spring ‘4 32 
Summer 85 52 
Fall 77 47 
Hartford, Conr Winter i) 34 
Spring 4 32 
summer sl 3y 
Fall te 4] 
Helena, Mont--. Winter 69 DH 
Spring 63 33 
Summer 67 55 
Fall 69 4s 
Houston, Tex Winter 22 31 
Spring 26 39 
Summer 1) 6s 
Fall 34 61 
ipolis, Ind Winter 24 30 
Spring 29 38 
Summer 57 7 
Fall 4) , 
Int itional Falls, Minn Winter 55 48 12 8 34 International Falls, Minn., June 1955 International Falls, Minn., June, 1957 
Spring 51 15 2 857 26 May 1957. May 1959 
Summer 70 4 2 65 26 
Fall 4519 7 47 25 
n, Miss Winter 54 32 31 58 34 47 38 Jackson, Miss., June 1944—January Jackson, Miss., June 1957- May 1959, 
Spring 9 10 7 27 62 49 1945, February 1956-May 1957 
Summer 6s 4 il 72 30 83 61 
Fall 70 24 2 58 38 70 62 
J nville, Fla Winter 69 27 44 59 4) 60 48 USNAS, Jacksonville, Fla.. June | WBAS, Jacksonville, Fla., June 1957 
Spring 53 4 $10 61 28 63 54 1955-October 1955; W BAS, Jackson- May 1959 
Summer 38 1 ll 56 23 75 44 ville, November 1955-May 1957, 
Fall 51 12 37 7 41 63 47 
J t, I Winter 46 25 29 56 33 Joliet, Il., March 1951-February 1953_| Peoria, Ill. June 1957-May 1959 
Spring 63 5 3 63 29 
Summer 82 7 4 65 34 
Fall 71 15 39 73 38 | 
kK s City, Mo Winter 39 47 
Spring 34 43 
Summer 47 7 
Fall 47 62 
K Vest, I Winter ; +s é¢o u Key West (Navy), Fla., June 1955 Key West (Navy), Fla., June 1957 
Spring 1 0 1 4 1 May 1957 May 1959 
Summer 0 0 1 0 l 
Fall 1 0 0 1 1 Data sparse throughout period, 
K e, Tenn Winter 55 34 | 
Spring 57 44 
summer 7s i 
Fall 73 53 | 
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Winter 
Spring 
Summer 
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Winter 
Spring 
Summer 


Spring 
Summer 
Fall 
Winter 
Spring 
Summer 
I ii 
Winter 
Spring 
Summer 
Fall 
Winter 
Spring 
Summer 
Fall 
Winter 
Spring 
Summer 
Fall 
Winter 
Spring 
Summer 
Fall 
Winter 
Spring 
Summer 
Fall 
Winter 
Spring 
Summer 
Fall 
Winter 
Spring 
Summer 
Fall 
Winter. 
Spring 
Summer 


Spring 
Summer 
Fall 
Winter 
Spring 
Summer 
Fall 
Winter 
Spring 
Summer 
Fall 
Winter 
Spring 
Summer 
Fall 

W inter 
Sprin 
Summer 


Fall 
Winter 






Fall 
Wirter 
Spring 
Summer 
Fall 
Wints 
Spring 
Summer 


all 





Spring 
Summer 
Fall 
Winter 


Spring 


Spring 
Summer 


Fall 


35 
13 
63 
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Percent frequency of 


Inversion 


Nighttime 


GMT Total Surface | Cloud 

15 OO 12 time Wind Cover 

7mph <3/10 

44 32 

3y 35 

HoH 44 

46 48 

22 41 

28 36 

47 1s 

36 MO) 

30 12 «38 36 43 37 

5 $ 5S 27 ) 44 

3 71 ST) 76 65 

12 » 61 41 72 63 
17 24 
2 21 
12 $2 
ww “) 4] 57 
37 $ 67 $1 
Is 4 Yu $4 
59 18 85 a) 

31 23 

J4 35 

s Me 

33 MO) 
y2 2 We ‘ 
54 0 SO ’ 

41 1 89 3 a 

ss Oo Ww 4 - 

17 32 46 20 0 3u 

11 5 Sl 27 Ss 

1 7 61 +) 83 a) 

9 24 64 42 77 ay 

a a t) 7t MO 

31 a 42 30 79 45 

14 a 19 “) 4] 

99 32 8 44 ”) 51 

is 6) 

4 4] 

SS 62 

67 me) 

35 33 

37 38 

65 uM) 

44 1) 

70 16 ¢ 44 87 38 

40 2 74 $4 77 0) 

$2 0 SO 33 st H4 

73 6 «76 4] 87 $8 

31 thal 

3Y 44 

hh 5Y 

on 61 

6 29 Op 24 47 SS 

2 10 47 1 $y +] 

2 21 65 13 75 49 

$ 19 75 22 4 is 
67 12 73 443 
Is) 1 65 $1) 
2 $ 4 7 
20 s 61 4 

-t $2 

$4 14 

ivi 5] 

Li) 1s 

35 3u 

+ 41 

it i) 

ot 7 

sl l4 

6s ”% 

7¢ 4 

sf i 

28 s7 

10 14 

67 4 

15 58 

3Y 37 

$1 4 

os ‘ 

7 54 

y 4 3s #2 1) 

F + © oa!) 73 1) 

6 10 72 30 Ww 4 

If $3 tid 42 s4 Mt 
15 28 4 29g 
uv 4 5 +0 
. =) ; 
WW 21 63 $4 
30 
$4 
$s 


Lake Charles, La., June 1955-May Lake Charles, La., June 1957—-May 1959 
1957 
Lakehurst (Navy N.J., June 1953 Station ed July 19. 
May 1955 
Lander, Wyo., June 1955-May 1957 Lander, Wyo., June 19 May 1959, 
Las Vegas, Nev., June 1955- May 1957.| Las Vegas, Ne June 1957-May 1959 
Little Rock, Ark., June 1955-May 1957_, Little Rock, Ark., June 1957-May 19 
Long Beach, ( f March 1954 Santa Moni Calif., June 1957-May 
February 195¢ 1959 
Medford, Oreg ine 19 M 57 Medfor Oreg., June 1957-May 1959. 
M ii, I , June 1955-May 1957 Miami, Fla., June 1957-May 1989 
M ind, Tex., June 1955-May 1957 Midland, Tex., June 1957-May 19 
Mavwell I mit , Ala., | Donnelly Field, Montgomery, A 
Decem be ve 1954 June 1957-Mavy 1959 
Selfridge Air Force Bast Mich., | Flint, Mich., June 1957-May 19 
June 1954-May 1956 


Observations at 0300 and 1500 GMT 


Station and period of re« 


ord used 


SEPTEMBER 





1961 





S. Navy 


Marine Climati 


Observations 





Atlas of the World Set 


it 0000 and 1200 GMT 


References 














SEPTEMBER 1961 











Station Season 
03 
Na ille, Tenn Winter 51 
Spring 0 
Summer 76 
Fall 74 
N ark, N.J Winter 
Spring 
Summer 
Fall 
N Orleans, La W inter 
Spring 
Summer 
ra - 
N York, N.Y Winter 21 
Spring $i 
Summer 45 
Fall is 
k, Va W inter 28 
Spring 31 
Summer m4 | 
Fall 29 
\ PI Net 6s 
r 74 
Vv 
Cali HH 
Spring 23 
Summer 4 
Fall s 
) Uta Winter 77 
Spring 63 
Summer 62 
Fall s4 
\k yma ( , Okla W inter 63 
Spring 70 
Summer SS 
Fall 85 
) ai, Net Winter 61 
Spring 5S 
Summer 70 
Fall 71 
i ia, Pa Winter 
Spring 
Summer 
Fall 
Ari Winter SS 
Spring 5S 
Summer 59 
Fall 5 
gh, Pa Winter 33 
34 
Idaho 
‘ Ma Vinter $ 
Spring rm 
Summer 7s 
Fall 70 
Oreg Winter 3th 
Spring 19 
Summer 14 
Fall iy 
R.I I 
Summer 
Fall 
Colo Winter 
Sprin 
Summer 
Fall 
N.( Winter 
Spring 
Summer 
Fall 
S. Dak . Winter 73 
Spring 64 
Summer 71 
Fall 77 
il, Va. > Winter 
Spr 
Sur I 
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GMT 
15 00 
24 35 

4 s 

l 6 
If 38 
7 18 
3 36 

1 2 

1 24 
28 40 
s 15 
3 ¥ 
21 30 
63 36 
Is 3 
9 4 
51 20 
83) tO«iW 
35 «14 
13 34 
62 #15 
ti 7 
23 ] 
28 2 
6y 3 
38 33 
15 6 
5 3 
27 16 
7 43 
23 Ss 
15 y 
ew 0 
Hs 0 
3s 0 
SS 0 
24 =O) 
6 - 
l 4 
15 25 
4 $s 
t 1s 
4 41 
16 O65 
il Ss 
“0 l 
3 0 
43 2 
7l 49 
29 3 
ll 12 
is 619 
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Inversion 





Total 
12 time 
3 31 
ti 27 
73 $2 
t4 un) 
24 14 
32 lf 
21 19 
30 ” 
4s 28 
Mw 23 
45 19 
At) 27 
74 43 
HY 32 
76 32 
67 43 
77 iy 
45 29 
4 tl 
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ied 
uea 
Station ind period of record used 
Observations at 0300 and 1500 Gm1 Observations at 0000 and 1200 cm 
Portland, Oreg., June 1954-May 195¢ Salem, Ore June 1957-Mavy 1959 
Ogden, Utah, June 154—-M 6 Salt Lake ( Ut June 57-M 
1959 
san Antoni rex June 1955-May Fan Antor rex une 1957-M 
Qh; 1959 
: Dieg« ( f June 4-May San Diego, ( June 1957-M ) 
ve 
Long ¥ Beacl ( M on4 s ( June y5T7-M 








Barksdale Air Force 
1955-December 1955 
port, Shreveport, I 


May 1957 





Spokane, Wash., June 1955 


Pampa, Fla., June 1956-May 1957 rampa, Fla., June 1957-May 1 
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Percent frequency of Station and period of record used 
In version Nighttime 
Station Season 
Observations at 0300 and 1500 GMT Observations at 0000 and 1200 amt 
OMT Total | Surface | Cloud 
08 15 00 12 time Wind Cover 
<7mph| =3/10 
toosh Island, Wash Winter 15 17 14 17 Tatoosh Island, Wash., June 1955 Tatoosh Island, Wast June 1957 
Spring 15 11 18 17 May 1957 May 1959 
Summer 20 19 37 27 
‘all 38 29 33 34 
edo, Ohio Winter 14 23 24 22 | Toledo, Ohio, June 1946—-December 
Spring 4 22 31 32 1950 
Summer 1 30) 63 51 
Fall 5 365 45 49 
eka, Kans Winter 39 «16 «48 28 39 45 Topeka, Kans., June 1955-May 1957 ropeka, Kans., June 1957-May 1950 
Spring 13 2 8 24 40) 2 
Summer 7 2 & 29 50 56 
Fall 24 «11 & 35 5O 60 
ilsa, Okla Winter 27 47 
Spring 28 4s 
Summer 41 60 
Fall 35 59 
icson, Ariz Winter 65 83 2 &9 52 Davis Monthan Air Force Base, Municipal Airport ucson, Ariz 
Spring 29 61 #1 8 41 Tucson, Ariz., June 1955-February June 1957-May 1960 (September 
Summer 19 15 4 74 31 1956; Municipal Airport, Tucson, 1958-June 1959 06Z & 18Z only 
Fall 77 70 #1 8 4s March 1956-May 1957 
Washington, D.C Winter 44622 30 48 28 37 43 | Washington, D.C., June 1955-May | Washington. D.C June 1957-May 
Spring 42 5 12 5i 23 38 38 1957 1959 
Summer 7 2 5 8&7 24 58 49 
Fall 56 8S 34 59 32 4 53 
Wichita, Kans Winter 19 49 
Spring 17 45 
Summer 13 5h 
‘all 13 63 
Winnemucca, Nev Winter 6 &2 48 Station opened May 1, 1956 Winnemucca, Nev.. June 1957-M 
Spring 0 &S8 41 1959 
Summer 0 92 38 
Fall 1 91 49 
Yakima, Wash Winter 76 26 
Spring 55 0) 
Summer 49 70 
Fall il 55 
stown, Ohio Winter 28 23 
Spring 10 33 
summer 72 57 
Fall 49 1s 
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Weather Note 


THE EFFECT OF DIFFERENTIAL FRICTION BETWEEN LAND AND WATER ON THE MOVEMENT OF DONNA IN THE 
VICINITY OF EASTERN NORTH CAROLINA 


EUGENE W. HOOVER 


Weather Bureau, Washington National Airport, Washington, D.C 


Manuscript received Marct 


The small-scale irregularities in the movement of hurri- 
cane Donna over eastern North Carolina in September 
1960 support the conclusion of Dunn et al. [1] on the 
effects of differential wind friction between land and water 
surfaces. They suggested that the increase in surface 
friction and the greater cross-isobaric flow over land con- 
tributed to an increase in pressure in the forward right 
quadrant. This increase in surface pressure caused a 
temporary deviation to the left and brief slowing in the 
forward speed of the hurricane. The duration of these 
small-scale irregularities was contingent upon the time 
required for the core of strongest wind to move inland. 
These vacillations were noted by Dunn et al. [1] in the 
tracks of hurricanes Connie, Diane, and Tone and by 
Dunn [2] in the track of Donna. 

The unusually large dimensions of the eve in the vicinity 
of North Carolina make Donna a remarkable hurricane to 
illustrate the concept of differential friction. The magni- 
tude of time and distance elements of deviations in the 
storm track are large enough to evaluate quantitatively. 

The track of the radar fixes of the center of Donna over 
eastern North Carolina (fig. 1) deviated inland to the left 
about one degree of latitude, which is approximately the 
same as the reported diameter of the eve. The forward 
speed (see data table in fig. 1) slowed from 27 kt. over 
water (segment E to F) to 13 kt. over land (segment 
B to C). 

The differential friction effect reversed as the hurricane 
moved off shore. The decreasing friction over water per- 
mitted the winds to increase and become more nearly 
parallel to the isobars off shore. According to the sug- 
gested hypothesis, this flow results in a depletion of air 
near the coast. This depletion of air causes the hurricane 
to be temporarily displaced to the right, and simultane- 
ously to accelerate forward briefly. The radar fixes of 
Donna in figure | illustrate the turn to the right and the 
forward acceleration to an average speed of 38 kt. between 
points C and D (see data table in fig. 1). Both the magni- 
tude and geographical position of the observed vacillations 


over eastern North Carolina are in agreement with those 
suggested by the differential friction hypothesis. 


27, 1961; revised May 1, 1961] 
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MOVEMENT OF HURRICANE DONNA 
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Track |Distance|E lapsed Time | Speed 
Segment] (N Mi) (Hr Min.) | (Kt) 
AtoD| 230 8:30 27 
AtoB; 110 3:50 29 
a BtoC 28 2:10 
Ctod 94 
EtoF 108 
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1.—Track of radar fixes on center of hurricane Donna, 2300 
GMT September 11 to 1930 GMT September 12, 1960 
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MOISTURE AND COTTON AT HARVEST TIME 
IN THE MISSISSIPPI DELTA 


J. A. RILEY 


Delta Agricultural Weather Project, U.S. Weather Bureau, Stoneville, Miss 
{Manuscript received February 28, 1961; revised May 18, 1961] 


ABSTRACT 


The moisture content of cotton at time of harvest has a direct effect on cotton quality. 
Weather, in particular relative humidity, strongly influences variation in cotton moisture 
content. This paper summarizes three years of observations of weather in and around 
cotton fields at harvest time in the Mississippi Delta. 

Detailed studies of the variation of relative humidity, temperature, wind, sunlight, dew, 
and evaporation have been made in leafed fields and in fields where leaves were removed 
artificially. Single weather variables and combinations of variables have been related to 
dew intensity and to the rate of evaporation. Surface weather maps are presented to show 
the average pressure patterns associated with extremes in dew and evaporation. A pre- 
liminary cotton picking guide, developed jointly by the Weather Bureau and the Cotton 
Harvesting Section of the U.S. Department of Agriculture, which gives the farmer an objec- 
tive estimate of when to start picking cotton on the basis of weather variation, illustrates 


‘ 
e 
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how weather information can improve farm efficiency. 


1. INTRODUCTION 


Drv weather is good cotton picking weather. In the 
humid Cotton Belt, over 8 million bales with a cash value 
of over a billion dollars were harvested in 1959. Some of 
this cotton was harvested too wet and the consequent 
reduction in quality resulted in a price drop of $10 to $15 
per bale; thus damp cotton in the humid belt is a problem 
of considerable economic importance. 

Wet weather in the form of excessive rain gives poor 
harvest conditions; mechanical equipment is inoperative 
when soils are water-logged; if the rain persists, maturity 
may be delayed until the plants are caught by frost; and 
frequent soaking of lint by rain stains the cotton. This 
type of wet weather damage has no control. Wet weather 
in the form of high humidity is damaging to quality if the 
cotton is picked while wet. If the cotton is not harvested 
until the lint has dried, higher quality results. This type 
of wet weather damage does have a control; a harvest 
program that restricts operations to periods of optimum 
moisture conditions of cotton on the stalk. 

When cotton is picked too wet or with a seed cotton 
moisture content above about 10 percent, staining is more 
By the 


middle of the harvest season, most gins are overloaded 


frequent and more trash adheres to the lint [1,2]. 


and many trailer loads of cotton are kept waiting at the 
gin vard. The delay accelerates bacterial activity which 
causes increased discoloration and general deterioration 
of fiber and seed. After wet cotton reaches the gin, fairly 
high temperatues are necessary for drying the lint to an 


acceptable moisture level. This makes the fiber brittle 
and staple length is often reduced. 

Research has shown that relative humidity exerts the 
greatest control over seed cotton moisture variation with 
the exception of actual rainfall. Wooten and Montgomery 
found that after a morning dew a relative humidity of 
approximately 50 percent is the dividing line between 
cotton dry enough to pick and cotton that is too wet. 
Under the usual harvest weather conditions, relative 
humidity can be used to give a realistic estimate of the 
seed cotton moisture content on the stalk [3,4]. 

Defoliation, or the removal of leaves by chemical 
application, is now practiced on over half of the cotton 
acreage in the humid belt. Its original purpose was to 
reduce leaf trash gathered by the ever increasing number 
of mechanical pickers [5]. In the last few vears, it has 
been recognized as an aid in reducing the moisture content 
of seed cotton [6]. This study illustrates the change in 
microclimate that results in cotton fields from defoliation. 

The objective of this study is the development of an 
understanding of moisture variation within cotton fields 
which will form the basis of a technique for forecasting 
“safe” picking conditions. 

The history of weather observation programs designed 
to increase the efficiency of the cotton industry dates back 
many vears. The Organic Act of the Weather Bureau 
enacted in 1890 states its mission in part as”. . . the report- 
ing of temperature and rainfall conditions for the cotton 
This Act is still in effect and its most 


” 


interests 











Hygrothermographs in (A) leafed and (B) defoliated 


cotton fields 


FiGurRe | 


recent implementation is the Delta Agricultural Weather 
Project, under which the work reported here was done. 
This study summarizes observations and tests made in 
cotton fields during the fall seasons of 1958 through 1960. 
All tests were made in cooperation with the Cotton Har- 
vesting Section of the U.S. Department of Agriculture and 
the Delta Branch of the Mississippi Agricultural Experi- 
ment Station. 

Observations were made in Delta Experiment Station 
cotton fields about one mile south of the Administration 
Building. The area is free from obstructions and typical 
of most Delta cotton fields. Deer Creek makes a wide loop 
around the area at a distance of about 0.4 mile north of the 
observation site, 0.8 mile to the east, and 0.3 mile to the 


south. Trees about 50 feet high are scattered along the 
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TaBLE 1|.—Cross section of leaf, stalk and lint in square feet per cubic 
foot of space occupied by defoliated and undefoliated cotton plants. 


Leaf cover Leal Stalk and rota! 

rea lint area area 

Defoliated 2.5 16. 5 19.0 
Undefoliated 116.5 16.5 127.0 


creek. Measurements were made in both defoliated and 
undefoliated cotton fields with similar crop cover extending 
several hundred yards in all directions from the instru- 
ment site. 

Because of the remarkably uniform growth character- 
istics of cotton grown in this experimental field and the 
scientific farming practices, cotton stands of great simi- 
larity were available for microweather observations during 
the three fall seasons. The average height of the cotton 
plants did vary slightly: in 1958-59 cotton averaged 5 feet 
high and in 1960, 4 feet. 
table 1 shows the 3-vear average values of leaf, and stalk 


Leaf density was very similar; 


and lint cross-section area in square feet for each cubic 
foot of space occupied by plants. The leaf area, or the 
area available for transpiration, in the defoliated field was 
only 2 percent of the leaf area in the undefoliated field. 
The total area of leaf, stalk and lint, or the area available 
for obstruction to wind and sunlight, in the defoliated field 


was 15 percent of the area in the undefoliated. 
2. RELATIVE HUMIDITY AND TEMPERATURE 


Relative humidity shows considerable variation in cot- 
ton fields and in turn is a key to the variation in seed 
cotton moisture. The air surrounding cotton plants re- 
ceives vapor from the transpiring plants, evaporation from 
moist soil and rainfall, and from mixing with moist air 
moving northward from the Gulf of Mexico. At night, 
the air frequently loses water vapor by dew formation 
but this is usually returned the following day by evapora- 
tion. Vapor losses of more importance occur with rising 
air currents, especially when accompanied by strong dry 
winds. Most of the diurnal variation in relative humidity, 
which is an important variable in designing a_ picking 
program, is caused by variation in temperature rather than 
in absolute content of moisture. 

Hygrothermographs were placed in instrument shelters 
at a height of 10 inches above the ground (see fig. 1) and 
at 54 inches above the ground (top of cotton averaged 60 
inches) in both defoliated and undefoliated cotton fields. 
Table 2 shows the average daily minimum relative hu- 
midity and maximum and minimum temperature and 
temperature range for 35 days without rain from Sep- 
tember 22—November 5, 1959. 

The minimum daily humidity deep within the plant 
zone averaged 9 percent lower in the defoliated field than 
in the leafed field. At the top of the plant zone, relative 
humidity was slightly higher in the defoliated field. 
Highest maximum temperatures, which are approximately 
concurrent with minimum relative humidity, occurred in 
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TABLE 2.—Average daily minimum relative humidity and mazrimum 
and minimum temperature and temperature range for 35 days 


without rain from Sept. 22—Nov. 5, 1959, Stonenlle, Miss. 


Minimum, Maxi- Minimum, Temper- 
Exposure relative mum temper- ature 
humidity | temper- ature range 
ature 
Percent F I I 
+-inch undefoliated 48.0 75.0 57.8 17.2 
(-ineh undefoliated 56.8 74.3 i) 16.2 
j-inch defoliated 48.8 79. 6 58.0 18. 4 
0-inch defoliated 47.7 80.9 57.4 23.5 


the defoliated field. The largest temperature range was 
at the low level in the defoliated field and the smallest 
range was within the plant cover of the leafed field. 

The vertical variation of relative humidity in the 
defoliated field was similar to the variation that occurs 
In this case, the relative humidity 
was slightly less near the ground than at the top of the 
plant. Geiger [7] calls this distribution the “dry type”’ 
and indicates that it is very rare, and indeed, he stated 


above bare ground. 


that it has been observed only in southern India, par- 
It is possible that 
the instruments used in this experiment are not precise 


ticularly over the black cotton soils. 


enough to detect a 1 percent difference; however, it was 
unusual that the relative humidity was not higher near 
the ground. 

The variation in relative humidity in a leafed field is 
The low level is sheltered 
from the sun and wind; thus moisture is not as free to 


mainly the result of exposure. 


disperse into the atmosphere as from the top of the plant 
zone. A smaller effect is the lower temperatures caused 
by shading. Using the 1959 data as an example (table 
2), if the maximum temperature of 74.3° F. observed at 
(the maximum at 


54 inches), the relative humidity of the air parcel would 


the 10-inch level were raised to 75.0° F. 
be reduced less than one percent. Since the difference in 
relative humidity between the two levels was actually 8.8 
the temperature 
account for the measured difference. The higher humidity 
at. the low level must have been caused by the lack of 


percent, difference in alone cannot 


mixing with drier air above the plant zone. 

Observations made during the fall seasons of 1958 and 
1960 conformed in pattern to the detail given above for 
1959. the humidity 
averaged between 8 and 10 percent higher within the 


For the 3-vear period, relative 
plant zone of the leafed field, and the humidity remained 
below 50 percent for one additional hour in the defoliated 
field. 


3. WIND AND SUNLIGHT 
Wind Air both the 


vertical and the horizontal plane to change the moisture 
Wind can completely 
change temperature and moisture characteristics of air 


variation. movement acts in 


level on and near cotton plants. 


surrounding cotton plants at times of air mass change and 


cause minor modifications at other times. Vertical air 
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Anemometer, hygro- 
top of 


Figure 2.—Weather tower in cotton field. 


thermograph, and Livingston atmometers located on 


tower. 


movement mixes the air near the plant with air of differing 
characteristics at some height above the ground. 

Vertical air 
however, the air stability 


movement was not measured directly; 
vertical 
feet 


was 


which is related to 


a hl 
rhe 


above the ground minus the temperature at 


movement was measured. temperature at 5 
35 feet 
this study. 


the measurement of air stability used in 


When the difference was positive, i.e., the temperature 
near the ground was higher than the temperature at the 
upper level, the air was more likely to undergo vertical 
motion by way of turbulent eddies. In this portion of 
the report, stability and wind are represented by daytime 
average (7 a.m. to 5 p.m. tst) and nighttime average 
(5 p.m. to 7 a.m. the next morning). From September 
30-November 8, 

1.45° and the nighttime stability, 
vertical motion was much more of a factor in the daytime 


1960, the daytime stability averaged 
2.18°; consequently 
than at night. 

The wind was measured at 40 feet above the ground on 
top of the observation tower shown in figure 2 and at 7 feet 
above the ground which was 3 feet above the defoliated 
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FIGURE 3 


cotton tops shown in figure 3. The wind on the tower 
averaged 5.27 m.p.h. and at the lower level averaged 
2.45 m.p.h. During the daytime, the upper wind aver- 
aged 6.50 m.p.h and the lower 3.65 m.p.h.; thus, the lower 
At night, the 


wind averaged 4.40 m.p.h at the upper level and 1.55 


level was 56 percent ol the upper level. 


m.p h. at the lower level. At night the lower wind was 
only 35 percent of the upper level. 


~ 
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FIGURE 4 
ibove the ground, above 4-feet high cotton, as a function of air 
stability at Stoneville, Miss., September 30-November 8, 1960. 


telationship between wind speed at 40 feet and 7 feet 


Anemometer, instrument shelter, and Livingston atmometer (observer refilling atmometer) in cotton fields 


The diurnal variation of wind with height above the 


cotton field is related to the diurnal variation of stability 


of the air column. Figure 4 shows the relation between 


winds at the upper and lower levels as a function of 


stabilitv. The top solid line represents 34 cases of sta- 


bility between +1° and +3° with an average stability of 
1.7°. The wind at 7 feet was a little more than half 


as strong as the wind at 40 feet. 

The dashed straight line represents 22 cases of al 
stability between +1° and 1°, with an average stability 
-0.2°. This average stability figure means that the 
temperature at the 5-foot level averaged nearly the same 
as the temperature at the 35-foot level. In this case 
there was practically no wind at the 7-foot level when 
the wind at 40 feet was 0.7 m.p.h. or less. The wind at 
7 feet was usually about half as strong as the wind at 
40 feet with a wind of 5 m.p.h. or more at the 40-foot level 

The lower solid straight line represents 24 cases when 
the air stability was between —1° and —6°, or an average 
of 3.2°. All of these cases occurred at night and 
represent very stable air; i.e., there was little or no 
vertical mixing of the air in the lower 40 feet above the 
ground. On this type of night, the wind was practicall 
zero at the 7-foot level with a wind of 2 m.p.h. at 40 fee 
With an average wind speed of 10 m.p.h. at 40 feet, tl 
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‘ABLE 3.—The average wind speed in defoliated and leafed fields at 
2 feet above the ground (cotton plants 4 feet high) and at 7? feet and 
{0 feet above the ground, Nov. 1-8, 1960, Stoneville, Miss. 


Undefoli- |Defoliated. 40 feet 
Time (LsT) ated. 2 2feetabove 7feetabove| above 
feet above| ground ground ground 
ground 
m.-5 p.m pees 0. 33 29 5. 22 9. 02 
).m.-7 a.m ncunbintah 0.01 0.04 1. 96 5. 66 


wind at 7 feet was only about 3 m.p.h. under stable 
conditions. 

The dashed curved line is a parabola representing the 
same conditions shown in the lower solid straight line 
(stable air). This curve is valid only for wind speeds of 
about 1 to 6 m.p.h. at 40 feet and is included to empha- 
size the stilling influence of very stable conditions; the 
air was practically calm at 7 
m.p.h. at the 40-foot level. 

During the period of November 1-8, 1960, wind readings 


feet with less than 4 or 5 


were made at a height of 2 feet above the ground in fields 
of defoliated and undefoliated cotton about 4 feet high. 
During this period, the defoliated field offered only about 
30 percent as much surface leaf, stalk, and boll obstruction 
to wind as did the leafed field. Table 3 shows the average 
daytime and nighttime wind speed in the defoliated and 
leafed fields at 2 feet above the ground and at the 7-foot 
and 40-foot levels. 

Practically no wind occurred within the cotton foliage 
at night unless the prevailing wind was relatively strong, 
at least 10 m.p.h. 
cotton averaged only 25 percent of the wind at 7 feet and 


Daytime wind speeds in the defoliated 


the wind in the leafed cotton was only 6 percent of the wind 
at 7 feet. At night, the wind speed in the defoliated 
cotton was only 2 percent of the 7-foot wind while the 
wind in the leafed field was less than 1 percent of the 
7-foot wind. 
Sunlight variation.—The difference between the amount 
of water evaporated from a white and a black Livingston 
atmometer gives an indication of sunlight or solar radia- 
tion [8]. 
white and 


Table 4 shows the total evaporation from the 
32 feet (unobstructed 
location) and at 26 inches (within a leafed field) and the 


black atmometers at 


percentage of increase in evaporation of the black over the 
white. 

At the unobstructed location, evaporation from the 
black atmometer totaled 57 
Within the leafed field, evaporation 
from the black atmometer totaled 16 percent more than 


percent more than from the 
white atmometer. 
from the white atmometer. The percentage increase in 
the leafed field 


centage increase at the unobstructed location. 


amounted to only 28 percent of the per- 
This 28 
percent can be assumed to approximate the amount of 
sunlight that reached the inner zone of a leafed cotton 
field. 
to the top of a 4- to 5-foot cotton plant was available to 
the middle of the plant. 


Less than a third of the sunlight that was available 
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TaRLeE 4.—Ebvaporation from Livingston atmometers and percent 
increases of black over white at an unobstructed area and within 
cotton foliage, Sept. 30—-Nov. 8, 1960, Stoneville, Miss. 





Evaporation Percent 
Increase 
Exposure of black 
White Black over 
whit 
32 feet unobstructed 1012 1591 7.2 
26 inches—within leafed cotton 300 349 16.3 


Both wind and vertical air currents are at a maximum in 
the daytime and are usually quite small at night. Wind 
increases with increasing height and with decreasing 
amounts of leaf obstruction. Sunlight is much reduced 
within the leaf Higher humidity and con- 
sequently higher seed cotton moisture content in well 
leafed fields, as compared with defoliated fields, is in part 
related to the absence of wind and sunlight in the leafed 


area. 


canopy. 


4. DEW INTENSITY 


Dew is a frequent and important cause of wet cotton 
and must be considered in designing a cotton harvest 
The daily probability of rain is under 20 per- 
cent while the probability of dew is over 80 percent in the 


program. 


Delta during the harvest season. 
Dew 
painted, 


thin, 
dew, 


the observation site on 
On 
approximately 2 cm.’ of water condensed on each 100 cm.? 
of surface. 
dew that formed on cotton plants and is the moisture 
equivalent of about 0.01 inch of rain. 
this report dew was classified as: none, light, moderate, 
Rain fell on 5 nights of the 40-day observation 
September 30-November 8, 1960 
intervals were omitted from this portion of the report. 


was collected at 


metal sheets. nights with heavy 


This is assumed to be fairly representative of 
For purposes of 


heavy. 
period, and these 

Huschke [9] lists the following conditions as favorable 
for dew formation: “ (a) a radiating surface, well insulated 
from the heat supply of the soil, on which vapor may 
condense; (b) a clear, still atmosphere with low specific 
humidity in all but the surface layers, to permit sufficient 
effective terrestrial radiation to cool the surface; and (ce) 
high relative humidity in the surface air layers, or an 
adjacent source of moisture such as a lake.”’ 

The availability of moisture is very rarely the limiting 
factor in dew formation in the Delta because of the many 
lakes and ponds in the area. During periods of extended 
drought, the availability of moisture is vastly reduced. 
Many of the lakes and ponds dry out and trees and crops 
under moisture stress do not supply much moisture by 
transpiration. During the drought years in the early 
1950’s the lack of local moisture was a limiting factor in 
dew formation; however, it has not been a factor in the 
3 vears of current Weather Bureau experiments in cotton 
fields. 

Relative 
determine the relationship between relative humidity and 


humidity in relation to dew intensity.~—To 
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TABLE 5.—Relative humidity at 5 feet above the ground, 1 foot above 
defoliated cotton, as related to dew intensity Sept. 5O0-Nov. 8, 1960, 


Stoneville, Miss 


Percent of Average relative humidity 


Average nights on preceding day 

relative relative 

Dew intensity humidity | humidity 
1 a.m did not 10 p.m.-12 7 p.m 1 p.m 
6 a.m reach 100 midnight ¥p.m 3 p.m 

percent 
lleavy loo 0 ON 92 42 
Mi t 100 0 6 AY 43 
Light v1 22 7s ti 53 
Nom SS 67 aS AS 44 


dew intensity on cotton fields, the humidity was measured 
at a height of 5 feet above the ground, which was 1 foot 
above the top of a defoliated cotton field. As shown in 
table 5, the relationship varied with the time of day. 
After midnight, there was a close direct relation between 
relative humidity and dew intensity. Between the pre- 
ceding afternoon’s relative humidity and dew occurrence 
next morning, the relation was inverse. There was very 
little relation between humidity during the early evening 
and dew intensity the following morning. 

During nights with heavy dew the air remained satu- 
rated; relative humidity equaled 100 percent from 1 a.m. 
to 6a.m. During 5 of the 7 nights with moderate dew, 
the humidity was 100 percent continuously during that 
period, but with light or no dew it was never continuously 
100 percent for the 5-hour period. 

On nights with heavy dew the average time of saturation 
was 12 p.m. (midnight) and the latest time of saturation 
With moderate dew, the air reached satura- 
On nights of light or no dew, 


was | a.m. 
tion about 14% hours later. 
the air was saturated for only very brief periods if at all. 

Table 5 shows that the highest humidities for nighttime 
periods before midnight occurred with moderate to heavy 
In the case of light dew and no dew, higher humidi- 
At such times, fairly high 


dew. 
ties occurred with no dew. 
humidity in the evening was accompanied by cloudiness 
but the clouds persisted through the night and prevented 
the ground from cooling and thereby prevented saturation. 

The last column of table 5 shows that low afternoon 
humidities often intensity. This 
merely indicates that clear days were usually followed by 


preceded high dew 


clear nights and the large diurnal temperature change 
gave low relative humidities in the day and high values at 


night. 


TABLE 6 Relation of percent opaque sky cover, excluding fog, as 
measured by the Gree nwood, Miss., FAA station during the hours 
of 7 p.m. through 6 a.m., to dew intensity at Stoneville, Miss. 


Sept. 30-Nov. 8, 1960. 


Percent of 
Dew intensity opaque 
sky cover 
Heavy ; s 
Moderate 5 
Light 4 


Nom ‘2 
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TABLE 7.—Relation of dew intensity and average wind speed in m.p.h 
from 5 p.m. to ? a.m. at two levels above cotton fields at Stoneville 
Miss.. Ne pl. 30 Nov. 8, 1960. 


Average speed Average 
7 feet above speed 
Dew intensity the ground 40 feet 
or 3 feet above ibove the 
defoliated ground 
cotton plants 
m.p.h mph 
Heavy ; 0.42 2. 46 
Moderate su 2 99 
Light 1.75 5. 46 
None 3.10 6.81 


Cloud cover in relation to dew intensity.—The variable 
used to represent cloud cover in this experiment was the 
percent of opaque sky cover, excluding obstruction by fog, 
from the hours of 7 p.m. through 6 a.m. as measured by 
the Greenwood FAA station 43 miles to the east. The 
distance between the two observation points introduces a 
small error; however, the difference is usually quite small. 
Table 6 shows that an average of 72 percent opaque sky 
cover occurred on nights with no dew formation. The 
effect of lesser cloud amounts was indeterminate. 

Wind speed in relation to dew inte nsity. Table 7 shows 
the average wind speed from 5 p.m. to 7 a.m. as recorded 
at 7 feet above the ground and at 40 feet above the 
ground. The speed at 7 feet is indicative of conditions 
3 feet above defoliated cotton plants and the speed at 40 
feet compares with wind observations at airport stations 
in the area. 
average winds at 40 feet of less than 3 m.p-h., light dew 


Moderate to heavy dew was associated with 
with average winds of 5-6 m.p.h., and no dew with an 
average of about 7 m.p.h. 

Table 8 shows the relation of dew intensity to air 
stability. On nights with no dew, the average stability 
was not far from zero; this means that the air in the first 
30 feet above the cotton plant was near neutral stability; 
i.e., Vertical motion was possible most of the time. On 
nights with dew, the air was normally quite stratified and 
little vertical motion occurred. Occasionally heavy dew 
occurred on cotton plants while the top of the tower was 
dry, indicating a very shallow temperature inversion. 
Since the inversion height is normally above the tower top, 
the stability measurements of this experiment can differ- 
entiate between dew and no dew but cannot be used to 


forecast dew intensity. 


TABLE 8. Relation of dew inte nsity to air stability as measured by 
the temperature at 5 feet above the ground minus the temperature at 
35 feet above the ground from 5 p.m. to? a.m. Sept. 30-Nov. 8, 1960, 


Stoneville, Miss. 


Dew intensity Average 
stability 

Heavy 2.49 
Moderate . —2. 63 
Light 3.90 
None 27 
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Combined effect of cloud cover and wind on dew intensity. 
Cloud cover and wind speed normally do not undergo 
extreme changes during the night, but there are exceptions. 
For example, a thundershower or a front may move 
through the area during the evening giving gusty winds 
and cloudy skies up until midnight, while after midnight, 
Clear calm weather 


Both cloud 


cover and wind speed are represented by a nightly average 


clear calm conditions may prevail. 
normally produces a moderate to heavy dew. 


in this experiment; thus it is not possible to account for 
the rare occurrences of radically changing weather during 
the night. 

Figure 5 shows the four classes of dew intensity as a 
function of cloud cover and wind speed at 40 feet. In 
table 6, the lower range of cloud cover shows little relation 
to dew intensity. Consequently, in figure 5, the classes 
of dew intensity for cloud cover of less than 20 percent 
are related almost exclusively to wind speed. 

As both wind and cloud cover are regularly forecast, 
it might be possible to use the relation shown in this 
graph to make regular forecasts of dew intensity. This 
has not been tested on independent data. 


5. EVAPORATION 


Evaporation from a water surface in an enclosed 
container is a relatively simple precess compared with 
evaporation from a field of cotton. Inside the container, 
an equilibrium is established between water molecules 
escaping from the water surface and those returning to 
the surface. In an open cotton field equilibrium is a 
rare condition; normally moisture is being taken away by 
evaporation or added by dew formation. 

Evaporation Was measured at 3 levels in both defoliated 
and undefoliated fields during the period 7:45 a.m. to 5:00 
p.m. on 35 days without rain, September 22—~November 5, 
1959. Table 9 shows that evaporation at the 10- and 
21-inch levels in the defoliated field was nearly twice that 
field. Even at the 54-inch level, 


which was only 6 inches below the top of the plant, 


in the undefoliated 
evaporation Was greater in the defoliated field. 

The evaporation observations during 1960 included the 
addition of two black Livingston atmometers. Because 
of the different radiation absorption characteristics of the 
white and black atmometers, the difference in evapora- 
tion between the two atmometers gives a measure of 


solar radiation. Observations in 1960 also included 
Tarte 9 lrerage et aporation from Livingston atmometers in 
oot-high defoliated and undefoliate d cotton fie lds on 35 days 
hout rain during the pe riod 7:45 a.m. to 5:00 p.m., Se pl. 22 

Vov. 5, 1959, Stoneville. Miss. 


Evaporation (cm 
Height above 
ground (inches 


Undefoliated Defoliated 


10 9.8 18.1 
21 10 
4 23. 1 24.9 
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Figure 5.— Dew intensity in relation to percent of opaque cloud 
cover from 7 p.m. through 6 a.m. and to average wind speed in 
m.p.h. from 5 p.m. to 7 a.m. at 40 feet above the ground, Stone- 


ville, Miss., September 30-November 8, 1960. 


evaporation readings, both black and white, made on the 
weather tower shown in figure 2. 

Table 10 summarizes evaporation readings made during 
the period 7 a.m.—5 p.m., September 30-November 8, 1960. 
Evaporation was directly proportional to the height above 
the ground and inversely proportional to the amount of 
obstruction. Within the two fields, 
evaporation was almost twice as great In the defoliated 


tvpes of cotton 
field, in agreement with measurements made the year 
before. 

The variation in evaporation between the 5-foot and 


TABLE 10 lreraqe evaporation fron Livingston atmometers at 
different heights and different er posures during the per od ? a.m 
5 p.m., Se pl. 30—-Nov. 8. 1960. Stoneville, Miss. The latte {wo 


columns relate evaporation at individual locations to evaporation 


at 5-foot level 


AV £ ( 
Height Type Exposur evapora Estimation of 5-foot itio 
tion sporatl yweth 
t 
32.2 feet W hite Unobstructed 25.3 Y=1.314+1.21\ 07 
$2.2 feet Black Unobstructed 39.8 Y =4.514+1.78\ ws 
5.0 feet W hits 1 foot above 19.8 
defoliated 
2, inches Whit Defoliated 15.9 ) 44+-0.78 \ i7¢ 
18 inches Whit Defoliated 14.0 Y=0.28+0.69 \ 
26 inches W hite Undefoliated 9 ) 1.744-0.53.\ 3 
18 inches W hite Undefoliated 7.5 ) 1.77+0.47 \ ) 
18 inches Black Undefoliat 8.7 ) 2.03+0.54\ HIS 
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32-foot levels was primarily due to stronger winds at the 
upper level. The variation between readings in fields of 
different leaf 


humidity, amounts of sunshine, and wind speeds. 


cover was caused by differing relative 

The fifth column in table 10 contains regression equa- 
tions which estimate the evaporation at the 5-foot level 
The last column 


relation 


as a function of the individual location. 


gives correlation coefficients which show the 
between the evaporation at the various observation points 
and the evaporation at the 5-foot level. The evaporation 
at each of the various observation points appears to be 
closely related to the evaporation at the 5-foot level. 
within the leafed field, although highly 
showed the poorest relation. This lower 


correlation is believed to have been caused by a slight 


Readings 
correlated 


change in leaf cover that occurred during the experiment, 
anv fundamental difference in relation to 
A light frost occurred October 21 
They dried up during 


and not by 
measured variables. 
and a few top leaves were killed. 
the last of October and subsequent frosts killed more 
leaves early in November. Besides offering less resistance 
to wind and sunlight, the dried leaves late in the season 
did not transpire as much moisture; thus the sheltered, 
moist climate that extisted through the first two-thirds 
of October was modified considerably by the end of the 
period. 

To obtain an estimate of how the evaporation within 
the leafed field increased during the experiment, the 10-day 
total evaporation at the 5-foot level was divided into the 
10-day total evaporation from the various levels within the 
cotton fields. Table 11 shows these measures. 

Evaporation in the undefoliated cotton field using the 
5-foot reading as a standard shows some increase during 
the third period and a sharp increase during the last 10-day 
period. This undoubtedly was due to the reduction of 
leaf cover, and to a lesser extent to the aging of the existing 
leaf cover as cold weather advanced. One interesting 
point is that the black atmometer in the leafed field had 
less evaporation than the white atmometer 8 inches 
higher during the first 30 days of the period; but, during 
the last 10 days, the condition reversed. This was caused 
by the increase in sunlight within the leafed field as the 
leaf cover gradually became less. For comparison pur- 


poses, evaporation at the same heights in the defoliated 


TARLE 11 10-day total evaporation (en from different instruments 
within cotton fields divided bu the 10-day total at the 5-foot level 


Defoliatel cotton Undefoliated cotton 
Periods White White White Black 
26 in IS in 26 in IS in 1S in 
Sept. 30-Oct. 9 0. S10 0. 719 0. 354 0. 268 0. 341 
Oct. 10-Oct. 19 Bu) 783 34 208) 340 
Oct. 20-Oct. 29 S16 719 487 391 459 
Oct. 30-Nov. 8 752 66S 544 SOS 572 
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field are included. Variations occurred there too, but 
were smaller and in the opposite direction to the changes 
in the leafed field; this tends to emphasize the drying 
effect that occurred in the leafed field with the advance 
of the season. 

Evaporation is measured by a number of different 
methods and each tends to give a different estimate of the 
same elusive weather element. Pan evaporation was 
measured regularly at the Stoneville weather station from 
April 1 to October 31. For the month of October, the 
evaporation from the 5-foot Livingston atmometer was 
related significantly, at the 1 percent level, to pan evapora- 
Total evaporation from the Livingston atmometer 
It is much 


tion. 
was 710 cm.’ and from the pan, 3.44 inches. 
easier to read 2-cm.® evaporation from the Livingston 
atmometer than it is to read the comparable evaporation 
of 0.01 inch from the pan. 
atmometer can be moved more easily. 


In addition the Livingston 
The aerody namics 
of the two different evaporimeters are different and it is 
believed that Livingston atmometers are much _ better 
suited to this type of investigation. 

To evaluate weather elements that influence evapora- 
tion, and thus drying power near cotton plants, various 
measurements of sunshine, relative humidity, wind, and 
temperature were related to the daily evaporation at the 
5-foot level for the period September 30—November 8, 
1960. 

Sunshine intensity in total langlevs per day and _ the 
average relative humidity from 1 p.m. to 3 p.m. at the 
5-foot level were found to be related significantly at the | 
percent level to the evaporation. Average wind speed 
from 7 a.m. to 5 p.m. at the 7-foot level was related signifi- 
cantly at the 5 percent level to the evaporation. No 
measure of temperature including maximum, minimum, 
and afternoon average, and no value of air stability was 
found that gave a reliable estimate of evaporation during 
the 40-day period in 1960. 

The following multiple regression equation relates cer- 
tain weather measures to evaporation: 
0.232.X; 


NX, = 13.781+-0.196.X; 0.034.N, 


where .Y, 1s total evaporation (em.*) at 5 feet from 7 a.m. 
to 5 p.m.; X, 
feet from 7 a.m. to 5 p.m.; .X; is average relative humidity 
at 5 feet from 1 


is average wind speed (m.p.h. 10) at 7 


(percent, with decimal point ignored 
p.m. to 3 p.m.; Yy is sunshine intensity in total langleys 
per day. 

The coefficient of multiple correlation for the above 
equation is 0.913, the standard error of estimate is 3.58. 
The standard error of estimate for each pair of variables is: 
wind and relative humidity, 4.04; wind and_ sunshine, 
4.17; 
combination of wind and relative humidity gives the 
}\ine 


and relative humidity and sunshine, 4.71. ‘The 


smallest error, while the combination wind and suns 
is only slightly higher. Sun and relative humidity give 
The reason is the interrelation of the 


the largest error. 
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Figure 6.—Hourly weather observations in or near cotton fields, Stoneville, Miss., October 12, 1960 A) Total hourly evaporation 
in em B) Total hourly solar radiation ir langleys hr.-! (C) Average hourly wind speed in m.p.h D) Temperature in °F. and 


relative humidity in percent. 


two variables; i.e., sunny days are generally dry (low rela- 


tive humidity) and vice versa. The partial correlation 
coefficients are: wind, 0.65; relative humidity, 0.52; and 
sunshine, 0.46. Wind is shown to be the most important 
single factor in the multiple regression equation, despite 
the fact that it had a smaller simple correlation coefficient. 
TI 


and sunshine constant 


e effect of holding the mutually related relative humidity 


is to increase the importance 
of wind. 

Figures 6 A-D show detailed measurements for one day 
during the experiment, October 12, 1960. High pressure 
with dry mild air centered over Tennessee extended into 
the Delta. 


sunrise cleared by 8 a.m. 


A few tenths of altocumulus clouds around 
Small, scattered, fair weather 
cumulus prevailed from 10 a.m. to5 p.m. Dew was heavy 
at 7 a.m., disappearing from defoliated cotton about 9 
Light 
dew formed again on exposed surfaces on both fields about 


a.m. and from undefoliated cotton about 10 a.m. 
6:30 p.m. Such weather is typical of a fairly good cotton 
picking day in the Delta. 

ligure 6A shows hourly evaporation totals at different 
different As indicated 
evaporation increased with increasing exposure (lack of 


heights and exposures. earlier, 
leaf and stalk obstruction) and with increasing height 
Besides the variation in total amount, 


it is of interest to note the time lag of maximum evapora- 


above the ground. 


tion that occurred within the cotton fields as compared 
The _ first 
evaporation occurred on the tower during the hour 8-9 


with more exposed locations. measurable 


a.m. 
inch level in the leafed field until 2 hours later, 10-11 a.m. 


Measurable evaporation did not occur at the 18- 


Evaporation at the exposed locations showed an early 
peak during the hours of 11 a.m. to 1 p.m., with another 
peak from 2-3 p.m. The early peak, although greatly 
damped in amplitude, extended to the 26-inch level in the 
leafed field. At the 18-inch level in the leafed field, there 
was a single maximum at 2-3 p.m., but the black at- 
mometer did show a rather sharp increase from 11 a.m 
to 1 p.m. 

Figure 6B shows the greatest hourly solar radiation, 63 
Seattered small 
cumulus clouds caused sunlight to fluctuate considerably 
after 11:20 a.m. 
all levels occurred from 2-3 p.m., which was several hours 


langlevs, occurred from 10 to 11 a.m. 
The maximum hourly evaporation at 


after the time of maximum solar radiation. 

Figure 6C shows an early maximum of wind from 11-12 
a.m. The greatest wind was measured from 2-3 p.m., 
which was also the period of greatest evaporation at all 
levels. 

Figure 6D shows the temperature and relative humidity 
variation throughout the day. Temperature reached a 
maximum and relative humidity a minimum about 11 
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FIGURE 7 Average sea level pressure in millibars at 12 a.m. 4) Heavy dew. B) No dew with southerly winds. (C) No dew because 


of rait D) No dew with northerly winds 


a.m. and only minor changes occurred until after 4 p.m. 
Leaf temperature readings as shown in figure 6D, were 
made when the sun was not obscured by clouds. Small 
electrically recording temperature probes, normally used 
to Ineasure temperature inside fruit, were taped onto the 
bottom portion of leaves. Leaves that were measured 
were alive and green but rather dry and apparently trans- 
piring very little. The maximum difference between the 
leaf temperature and the air temperature at 5 feet 
occurred shortly after noon, during a period of maximum 
sunlight. At that time, the leaf temperature was 11 
higher than the surrounding air temperature. Diiferences 
between leaf and air temperature on this day were typical 
of readings made on several other fairly clear days. The 
leaf temperature was a little below the air temperature 
early in the morning until the dew was nearly all evapo- 
, then the leaf was 
warmer than the air until just before sunset. When the 
leaf no longer received direct solar radiation (usually 4-5 
During the period 


rated (usually between 8 and 10 a.m. 


p.m.), it became colder than the air. 


of dusk, with clear skies, the leaf was about 5° cooler than 


the air, which of course is a necessary condition for dew 





Stoneville Miss., September 30—-November 8. 1960. 


formation. On the day shown in figure 6D, dew had 
begun to form on well exposed leaf and = stalk surfaces 


by 6:30 p.m, 


6. RELATION OF DEW AND EVAPORATION 
TO PRESSURE PATTERNS 


Dew intensity and sea level pressure patterns. Figure 7 
shows the average sea level pressure pattern for 15 cases 
of heavy dew. High pressure is the significant feature; 
however, there are two different locations of the high 
center that give dew in the Delta. The small knob of 
high pressure that appears over central Mississippi in 
figure 7A reflects four cases where a high pressure system 
just to the south of the Delta was associated with heavy 
dew. The larger and more pronounced high pressure 
system to the north represents the average center on the 
nine other nights of heavy dew. 

Figure 7B shows the average sea level pressure pattert 
for four nights of no dew. Figure 7C shows the pattern 
of five nights on which rain fell and no dew formed and 
figure 7D shows the average pattern for the two remaining 


nights of no dew. 
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The nights comprising figure 7B were grouped together 
The 
vind at 40 feet averaged 4.3 m.p.h. and the cloud cover 
North wind prevailed during the 


vecause the wind was from the south in the Delta. 


veraged 73 percent. 
ights grouped in figure 7D; here the wind was much 
stronger, 11.8 m.p.h., and the cloud cover was about the 
same, 72 percent. 

It is interesting to note that the three maps of no dew 
seem to form a series showing a trough of low pressure, 
or a front, moving from northwest through the Delta. 
With an approaching front, increasing cloudiness prevents 
dew. When the front is overhead, rain combined with 
cloud cover prevails, and immediately after the front has 
passed, gusty winds prevent dew. 

Strong frontal passages as far south as the Delta are 
rare in September and the first half of October; during 
that period cloudiness and wind associated with weak 
fronts or instability lines are the usual conditions limiting 
dew formation. Later in the vear, active fronts and more 
pronounced low pressure systems are the limiting factors. 
In most cases, cloudiness is the limiting factor in dew 
formation before the middle of October. It is very rare 
After 


the middle of October, wind is the main limiting factor 


to have a night with no dew before mid-October. 


and nights with no dew become somewhat more frequent. 
Figure SA 
shows the average sea level pressure pattern at 12 a.m. 


Evaporation and sea level pressure patterns. 


on the days when evaporation totaled 30 cm.’ or greater. 
Figure 8B shows the same data for six days when the 
The first map, 
showing strong drving conditions, features a low pressure 


evaporation totaled less than 10 em. 


trough from the Texas Panhandle through southwestern 
Nebraska to Lake Superior, a high pressure center just 
south of the Delta, and relatively strong pressure gradient 
just north of the Delta that would bring dry air and 
Figure 8B, showing 
poor drying weather, features a weak trough from southern 


rather gusty winds from the west. 


Texas through western Arkansas to Lake Michigan, high 
pressure over Colorado and the Virginias, and a weak 
pressure gradient over the Delta that would bring light, 
moist south to southeast winds over the Delta. 

The charts indicate that the strength of the westerly 
wind over the area just west of the Delta is the most 
important synoptic feature. Strong westerly winds give 
good drying conditions, while southeasterly winds give 


poor drying conditions. As a measure of the strength 
of the westerly winds, the 12 a.m. sea level pressure at 
(Pats) the 


Sioux City (pgyx) and the following regression equation 


Galveston was subtracted from pressure at 


developed: 


Kvaporation (at 5 feet) =19.49—0.07 (pgux— Pets) 


The correlation coefficient relating the sea level pressure 
difference between the two stations to the evaporation 
has been about the same magnitude as the correlation 
coeflicient relating the local wind speed at 7 feet to the 
evaporation reading. 
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FiGuRE 8.—Average sea level pressure in millibars at 12 a.m. (A 
Six days with evaporation from Livingston atmometers at 5 feet 
above the ground (one foot above defoliated cotton) of 30 em. 
or more. (B) Six days when evaporation was less than 10 cm.? 
Stoneville, Miss. September 30—-November 8, 1960. 


7. ESTIMATING PICKING TIME 


Figure 9 shows a cotton picking guide designed jointly 
with the Cotton Harvesting Section of the U.S. Depart- 
ment of Agriculture. The graph is based on average 
conditions measured during the falls of 1958 and 1959. 
The three major variables are: time of vear, leaf cover, 
Minor modifications are made by 
By using the guide, 


and dew intensity. 
considering cloudiness and wind. 
an individual farmer in the Mississippi Delta can make 
an early morning estimate of when his cotton will be 


dry enough to harvest. 
8. CONCLUSION 


Observations taken within cotton fields during the 
harvest season indicated that relative humidity variation 
was strongly influenced by the amount of leaf cover. 
Daytime humidities were higher deep within a leafed 
plant zone than near the top. Relative humidity at 
midday averaged 9 percent lower in defoliated fields 
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September October November 
20 C (@) 20 31 10 20 30 
10 a.m. 10 + 3 u - + + 10 a.m. 
Heavy dew -- Undefoliated wan — 
| a om 1” - 
9 a.m. t ae T 9 a.m. 
—e —_ 
Pee iid Heavy dew -- Defoliated a 
a ao diag a —_ 
8 a.m. — 8 a.m. 
qua 
ona cg al 
ae Light dew -- Undefoliated 
—_—_ 
—_ 
7a. -—— 4 7 a.m. 
Light dew -- Defoliated 
6 a.m. ! u 6 a.m. 
Rain past 24 hours -- Use heavy dew lines 
Gusty winds -- Subtract 30 minutes 
Cloudy sky — Add 30 minutes 
Raining — DO NOT PICK 
Figure %.—Daily guide for beginning mechanical cotton picking, Mississippi Delta Area, designed jointly by the Cotton Harvesting 


Section, ARS, and 


than in leafed fields. It normally stayed below 50 per- 


cent, an estimate of “‘safe’’ picking conditions, for an 
additional hour in defoliated fields. 

Wind and sunlight within a cotton field were found 
to be inversely proportional to the amount of leaf cover. 
The decrease of wind and sunlight within the cotton 
field was related to the variation of relative humidity 
at the same level. Sunlight was much reduced within 
a leafed cotton field. Less than a third of the sunlight 
measured at the top of 4-foot-high cotton penetrated 
to within 14; feet of the ground. Wind increased with 
increasing height above the ground and decreasing amounts 
of obstruction. At 2 feet above the ground in the cotton 
fields, daytime winds in a leafed field averaged only 6 
percent of the wind at 7 feet above the ground and in a 
defoliated field they averaged 25 percent of the wind at 7 
feet above the ground. 

Dew is the most frequent wetting factor in cotton 
fields at harvest time. 
skies, light winds, stable air, and high nighttime humidities 


Heavy dew occurred with clear 


In the cases studies here, a large high 


in cotton fields 


the Weather Bureau. 


pressure system at sea level over, or north of, the Delta 
was associated with heavy dew; low pressure, especially 
in the form of a northeast-southwest trough in the area, 
was associated with no dew. 

High values of evaporation were favored by strong 
sunlight, low humidities, and strong winds. Evapo- 
ration was twice as great in a defoliated field as in a 
leafed field. In addition to the difference in total eva- 
poration, there was a time lag in evaporation in the 
leafed field. A strong westerly flow of air, as indicated 
on the sea level map, was the best synoptic predictor 
of high rates of evaporation in the Delta. 

The cotton picking guide designed jointly with the 
Cotton Harvesting Section of the U.S. Department of 
Agriculture gives the farmer an early morning estimate 
of when his cotton will be dry enough to harvest. Further 
studies relating weather variations to seed cotton moisture 
variation and objective techniques for forecasting dew 
and evaporation are expected to refine the guide with 
an end in mind of giving the farmer a technique for plan- 
ning his picking activities at least 12 hours in advance. 

















EPTEMBER 1961 
ACKNOWLEDGMENTS 


The studies reported in this paper are part of a joint 
project with the Cotton Harvesting Section of the U.S. 
Department of Agriculture and the Delta Branch of the 
\ississippi Agricultural Experiment Station. Grateful 
acknowledgment is expressed to Engineers E. B. William- 
son and O. B. Wooten for design and implementation 
of the project. Valueble suggestions were made by 
Dr. Wm. L. Giles, Station Superintendent Vice 
President for Agriculture of Mississippi State University), 


(now 


and by the Cotton Quality Improvement Committee 
of the Delta Council. Timely insecticide applications 
by the Entomology Section saved the experimental 


site from possible damage on several occasions during 
the 3-vear project and valuable assistance was received 
from the Cotton Physiology Section. Mr. L. L. Means, 
Chief, Public and Agricultural Forecasts Section, Weather 
Bureau, provided encouragement and numerous helpful 
suggestions throughout the project. 


REFERENCES 


1, O. B. Wooten, and R. A. Montgomery, “The Effects of Relative 


Humidity and Spindle Moisture on Machine-Picked Cotton,” 
Cir. 204, Mississippi State University, April 1956. 


MONTHLY WEATHER REVIEW 


9 


“~. 


353 


O. B. Wooten, R. A. Montgomery, and J. A. Riley, “Relationship 
of Moisture to Machine Picking of Cotton,’’ Paper presented 
at Southeast Meeting, American Society of Agri- 
cultural Engineers, Feb. 1959. 


Sectional 


. A. Riley and E. B. Williamson, ‘Relative Humidity in Cotton 
Fields at Harvest Time,’’ Bul. 581, Mississippi State Univer- 
sity, July 1959. 

E. B. Williamson and J. A. Riley, ‘‘The Interrelated Effects of 
Defoliation, Weather, Mechanical Picking on Cotton 
Quality,” Paper Winter Meeting, American 
Society of Agricultural Engineers, Dec. 1960. 

G. B. Crowe and H. R. Defolia- 

tion,’ Bul. 552, Mississippi State University, June 1955. 


and 
presented at 


Carns, ‘‘Economics of Cotton 


a Modification 
Modified Micro- 


Phy siology 


. A. Riley and E. B. Williamson, ‘ Defoliation 
of Cotton Harvesting Conditions Through a 
climate,” Cotton Defoliation 
Conferrence, Greenville, 8. C., Jan. 1961. 


Proceedings, and 


R. Geiger, The Climate Near the Ground, Translation by M. N. 
and others of the 2d German edition of Das Klima 
1956, The Harvard 


Stewart 
der Bodennahen Luftschicht issued Aug. 
University Press, 494 pp., 1957. 

B. E. Livingston, ‘““Atmometers of Porous Porcelain and Paper, 
Their Use in Physiological Ecology,” Ecology, vol. XVI, No. 
3, July 1935. 

R. E. Huschke (Editor), } 
Meteorological Soc iety, Boston, Mass.. 


American 
163 


Glossary of Meteorology, 


1959, 638 pp. (p 





304 MONTHLY WEATHER REVIEW 


Weather Note 
THE “DOUBLE EYE’ OF HURRICANE DONNA 


C. L. JORDAN 


LT. FRANK J. SCHATZLE, USN 


On the afternoon of September 6, 1960, hurricane Donna 
presented a striking picture on the radarscope of a U.S. 
Navy hurricane reconnaissance aircraft! as illustrated 
by figure 1. The structure of the inner portion of the 
storm as indicated by the radar was very unusual in that 
not only was the exact center of the storm clear of precipi- 
tation echoes but there was also a relatively wide ring in 
the storm interior with little or no radar return. Descrip- 
tions of the core of this storm given to the press following 
the flight led to newspaper stories of a “‘double eve,” and 
at least one newspaper headline referred to “two-eyed”’ 
Donna. These accounts were misleading since they failed 
to mention that the “eves” were concentric. In the 
transmitted weather reports, the storm core was described 
in terms of two concentric wall clouds and the inner 
diameter of each of the radar rings was given. As will 
be discussed below, this “double eve” structure persisted 
for at least several hours and connections between the 
inner and outer rings, if present, were very weak. At the 
time of the photographs presented in this note, hurricane 
Donna was located at approximately 21.5° N., 69.5° W., 
and was moving toward the west-northwest at a speed of 
about 8 kt 

A low-level penetration of the storm core was attempted 
from the south-southwest shortly after the time of the 
photograph presented in figure 1. The aircraft ? was 
at an altitude near 1,000 feet and penetrated to a point 
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about 32 miles from the center of the inner ring or just 
outside the edge of the outer ring. At this point the Fieure 1.—A photograph of the scope of an APS-20E radar in a 


penetration attempt was abandoned because of excessive 


engine cooling due to the heavy rain. Surface wind speed ? , 
. the center of the photograph 


estimates based on the state of the sea reached 100 kt. Turk 
at 45 miles from the storm center and the maximum value of hurricane Donna. The aircraft 
was entered as “120+”. The central sea level pressure 1500 ft. and on a heading of 


of the storm at this time was near 940 mb. 
A photograph of the scope of the vertically scanning 
radar taken during the attempted low-level penetration 





1 Of the Airborne Early Warning Squadron Four then located at the Naval Air Station, 
Jacksonville, Fla 


2 The aircraft was of the Super Constellation type, U.S. Navy designation WV-3 ference from other electronic equipment. 


U.S. Navy hurricane reconnaissance aircraft 
The aircraft position (at 
10 miles south of Grand 


mately 2100 GmT, September 6, 1960. 


Island and about 105 miles west-southwest of the cente 
an altitude of about 
Some of the echoes to the 
south and southwest of the aircraft 
were from the mountains of Hispaniola. 


mately 1.5° in the horizontal 
symmetrical spiralling pattern of light lines was due to 


taken at approxi- 


a range of 75-150 miles 


APS-20E is ai 


S band (10.4 em.) radar with a peak power of 2.0 megawatts 
The beam width, as defined by the half-power points, is approx 
the vertical. Th 


ter- 
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Figure 2.—A photograph of the range-height scope of the APS-45 

radar in a U.S. Navy hurricane reconnaissance aircraft taken 
about 30 minutes after figure 1. Height is shown on the vertical 
scale and slant range from the aircraft on the horizontal scale. 
The aircraft was at a relative position given approximately by 
“x”? on figure 1 and was at an altitude of about 1,000 ft. The 
view is directly across the storm center toward the north-north- 
The APS—45 is an X band radar (3.2 em.) with a peak 
power of 450 kilowatts. The beam width, as defined by the 
half-power points, is approximately 3.0° in the horizontal and 
1.0° in the vertical. 


east. 


fig. 2) gives additional information on the structure of 
the storm core. This range-height view, taken from a 
relative position indicated by the “x’’ in figure 1, 
directed 


north-northeast. 


was 


across the eve from south-southwest toward 
The cross-section of the inner ring is 
clearly shown but, because of attenuation, the outer ring 
This radar, the 


APS-—45, is subjeci to much greater attenuation than the 


on the opposite side is barely discernible. 


APS-20E, from which figure 1 was obtained, because of 
the shorter wavelength which is used. The inner wall 
loud is shown on the range-height scope as nearly vertical 
ind the echoes extend above 30,000 ft. The diameter of 
; 10 miles at the lower 
evel to about 13 miles at the high levels. 


this inner ‘‘eve’”’ varied from about 


The “clear”? zone between the inner and outer rings on 
igure 2 is confused to some extent on the side nearest the 
ircraft by return from the sea which shows up in the 
The small faint 


cho in the “clear” space on the far side of the inner ring 


iorm of a horizontal fanshaped pattern. 


ould perhaps be associated with the small echo in the 
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FIGURE 3. 
radar taken at 
The aircraft was at the position indicated by the ‘“‘“+” and at 


A photograph of the plan-position scope of the APS—45 


approximately 2215 GMT, September 6, 1960 


an altitude of about 8,300 ft. The 20-mile and 40-mile range 


circles to the southeast of the center and the 60-mile circle in 
the north-northeast sector can be seen faintly. 


same region on figure 1. The time difference of the two 
photographs is only about 30 minutes. 

The radar appearance of the inner ring was more charac- 
teristic of convective type echoes than the outer one and 
it extended to somewhat greater heights. The maximum 
height of the inner ring during the several hours that it 
could be seen on this flight was about 45,000 ft., while 
the clouds forming the outer ring appeared to be more 
stratiform in character and gradually died out on the 
30,000 ft. In 


greater strength and height of the inner echo, it could be 


range-height scope above view of the 
argued that, at this time, the inner radar ring should be 
considered as the true wall cloud of the eve. 

Following the low-level penetration attempt, a climb 
to the 700-mb. level was made and the storm core was 
penetrated from the northeast with relatively little diffi- 
culty. The space between the radar rings was free of 
clouds in a layer which was estimated to extend from 
3,000 ft. to 25,000 ft. 


limited area was visible during the penetration and, un- 


about about However, only a 





) 


fortunately, the exit was made after dark. The stra- 
tocumulus between the two radar rings consisted of well- 
formed concentric horizontal rolls and appeared to be 
unusually low. The ocean surface was not visible through 
the stratocumulus so that we have absolutely no infor- 
mation on the distribution of wind between the inner and 
outer rings. Unfortunately, the research aircraft of the 
National Hurricane Research Project did not enter Donna 
on this day. 

There was a thick layer of cirrus or high altostratus 
extending out from the inner core and covering the area 
between the two rings. The inner core of the storm was 
described by one of the pilots as resembling a giant “‘tastee 
freeze’. Inside the inner ring, the stratocumulus was 
nearly overcast but the sea surface could be seen in spots. 
Cirrus was absent over most of the inner “‘eye’’ and there 
were no clouds above the stratocumulus which, in this 
area, extended to about 6,000 ft. 

A photograph of the APS-45 scope taken during the 
penetration of the inner radar ring (while the radar was on 
horizontal scan) is shown in figure 3. The “clear” ring 
shows up nicely but it is shown as being somewhat nar- 
rower than indicated in figure 1. This difference is not 
surprising in view of the shorter wavelength of this radar 
and the fact that the aireraft was much closer to the 
‘clear”’ ring 


cloud systems. The narrowest portion of the 
and the greatest suggestion of a connection is shown on 
the north side of the center. This is almost directly op- 
posite to what was shown in figure 1 but there is a time 
difference of a little more than 2 hours between these two 
photographs. 

The “double eye” structure was observed throughout a 
5 hours when the afternoon flight of 


The inner 


period of about 
September 6 had the hurricane on its radar. 
ring appeared to be weakening as the aircraft moved away 
from the storm but it was noted on the subsequent flight 
a few hours later. Limited information is available from 
this later flight because the radar was not functioning 


306 MONTHLY WEATHER REVIEW 


SEPTEMBER 11%} 


properly. The eye reorganized during the night and on 
the following day an orthodox eye was found with a 
diameter of about 25 miles. About 12 hours prior to the 
time of figures 1 and 2, there was no evidence of a “‘double 
eye” and the eye diameter was reported as 20 miles. 

A rather poorly documented case of a “double eye” was 
reported in a hurricane of 1947 [3] and the same type of 
structure has been reported in several typhoons in recent 
years. Fortner [1] described typhoon Sarah of 1956 at 
the time when it had ‘fan eye within an eye”’. 
of the spectacular typhoon cases, the inner convective 


In some 


system appears as a giant cumulonimbus centered in a 
large “eye”. The “inner eye” within the cumulonimbus 
vortex has been observed to be as small as 2 miles. In 
these cases, winds well in excess of typhoon force were 
observed in the relatively clear space separating the inner 
and outer wall clouds. In other cases in recent years [2], 
the radar has indicated a very narrow ring of relatively 
weak return just outside a ring-shaped wall cloud but 
these weak areas were not apparent from visual observa- 
tions made from the reconnaissance aircraft. 

To our knowledge, all cases of “double eyes” in intense 
There 


are often radar features, the so-called “false eyes’’, out- 


tropical cyclones have been of the concentric type. 


side the storm core which could be misinterpreted as the 
These, 
however, do not have closed circulations about them and 
are hardly detectable, if at all, in the wind and pressure 
fields. 


eve of the storm by inexperienced observers. 
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DISTRIBUTIONS OF WEEKLY AVERAGES OF DIURNAL TEMPERATURE 
MEANS AND RANGES ABOUT HARMONIC CURVES ' 


CHRISTOPHER BINGHAM 


The Connecticut Agricultural Experiment Station, New Haven, Conn 
(Manuscript received April 11, 1961] 


ABSTRACT 


The distributions of weekly averages of diurnal temperature maxima, minima, means, 


and ranges are found to be non-normal, although the errors in using the normal distribution 


do not impair the usefulness of derived estimates of probability. 
the variance are estimated by harmonic regression. 
abilities for any week from as few as five parameters. 


The central tendency and 
This enables the estimation of prob- 
A three-term harmonic curve fitted 


to individual vears is adequate to describe the course of temperature. 


i, INTRODUCTION 


Useful probability statements about the occurrence of 


diurnal temperature means or cannot be made 
without knowledge of the form of the distribution of the 
values and a practicable means of estimating the param- 


The present paper 


ranges 


eters which specify the distribution. 
explores the distributions of the weekly average of the 
(= (max+min)) the 
weekly average diurnal range, as well as the related dis- 


daily mean temperature and 
tributions of the weekly average diurnal maxima and 


minima. 
2. HARMONIC REGRESSION 


The fundamental tool which will be used in this develop- 


ment is harmonic regression. It is well known that any set 


of data, 7, Io, To, at equally spaced times f,, fo, . 


t), may be exactly fitted by a series of the form 


n 


y=ay+- >) A, cos (pt—@,), t measured in degrees. (1) 
p=1 


This is the sum of cosine curves, each with semi-amplitude 
A, and time of maximum t=4¢,/p. Equation (1) can also be 
written in the form used herein, 
n 
Y¥=Qo7t js i (a 
p=1 


p COS pt- b, sin pt) (2) 


wher 


1,=.1, cos ¢,; 6,=A, sin ¢,; and a2+ 62=A?, p=1,.. .,n. 


(3) 
Such a sum will be called an n-termed Fourier series, and 
a, cos pt+-b, sin pt) will be called the pth term. In many 


rk was supperted in part by funds of a regional project in agricultural clima- 


applications it is found that very few terms are needed 
to give an excellent fit to the data, so that the residuals 
from the curve are of the same magnitude as the basic 
For instance, Craddock [6] ob- 
served that over much of the Northern Hemisphere, a 


errors of observations. 


two-term series provided an adequate fit to the mean 
monthly temperature. 
variability between vears of the coefficients of the two- 
term curve which best fit each year. Without this knowl- 
edge, a proper error term for significance tests for the 
Bliss [4] remedied 
this deficiency, describing a technique paralleling the 


He did not, however, examine any 


reality of given terms is not available. 


standard analysis of variance for the fitting of orthogonal 
polynomials. 

The mathematical model underlying Bliss’s analysis is 
the following. Each observation for the jth unit of time 


in the ith vear, y;;, is considered as a sum 


y (dp + aio) 4- (A, +a) COS t;> (6,+8;,) sin t, 

+-(d2+aj2) COS 2t;+-(b.4+-B,2) sin 2t, 

.. +(a,+a;,) cos rt,+(b,+8,,) sin rt;+e, (4 
where t;=j/k 360°, j=0, 1, 2, . . ., k—1 (if the units are 


weeks, k=52). The e distributed 
normal deviates with zero means and common variance o”. 
«> Gir, Bi, Bis, eo 8 ef re] r} 


2r+1 


are independently 
The vectors (a 0, At, are inde- 
pendent observations from a variate multivariate 
distribution with zero mean vector. This model covers 
both the case of a single curve applicable to every year 
Viag)= ... V(B;,) =0), 


and the case of random variation of the curve from year 


(when the variances V(aw) 
to vear. In the latter case the variance components of the 
For details on 
His 
The 


coefficients enter into the variance of y 
the analysis of variance see Bliss’s [4] bulletin. 
analysis is designed to test several hypotheses. 
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simplest such hypothesis is that a7+6/=0, for particular 
values of j7. One can also test the adequacy of r terms to fit 
the data, and the reality of the between vears variance 
components associated with the coefficients. Significance 
tests are made using the F-test for variance ratios. 

It should be pointed out that a certain arbitrariness 
This is the choice of a 
To illus- 


must occur in the above model. 
zero point from which we measure the time f. 
trate this, let us suppose that 


y y+ ayo) + (A, +a) cos t;+(b;+B,,) sin t, 


+ (€,+aj2) cos 2t,+(b2+Bi2) sin 2t;. 
Then if 0<6< 360°, y can equally well be expressed by 


+-B,,’) sin t,’ 


, , , 
a, ) cost, 4 (b, 


Yis= (Ao TA) TQ 
+ (a,’ +a 2’) COs 2t,’ + (b,’ t Bi2’) sin 2t @ (5) 
where ¢,/=t,+0 and a,’=a, cos 0@—b, sin 9 
b,’ =a, sin 0+4, cos 0 


ly’ =a» cos 20—b, sin 20 


b.’ =a» sin 20+), cos 20, 


with a similar relationship holding for the a’s and the 
3’s. If one fits a curve of the form (5) to data, exactly the 
sume curve is obtained as when the more conventional 
form is fitted. It is easily seen that the variance com- 
ponent contributed by a; will, in general, be different from 
that contributed by a’. However, the sum of the com- 
ponents for a given term will be unchanged. That is, 
Via) + Vib) = Via’) +V(8,’). Further,the semi-ampli- 
tudes A, will be unchanged by the change of origin. Since 
the analysis of variance does not separate the two com- 
ponents of each term, this is sufficient to ensure that the 
analysis is invariant under any choice of origin. 

a,*,...., 6,*, of the coefficients, 


ir 


Estimates, d4*, 
(dg+aw), (Q;+-ay),... , (6,4+8:,), for the ith vear are 


easily made as follows. 


Mi 


(6) 


to 


Sin >> y:, cos mt, 


ih 


— 
* 
bo 
4 
dus 


‘wen Mt, M—i2,...,%75 8 ; 360°. 


Estimates for the mean coefficients dp, a;, ... , 6, are 
obtained by averaging theabovea »* sand b,,,*’s over Vears. 
These are least squares estimates, and hence if the assump- 
tions stated are fulfilled, they are the best unbiased esti- 
mates; further, if the error component is normally distrib- 
uted they are also maximum likelihood estimates. The 
estimates remain unbiased even when the restriction of 


independence of residuals is removed. 
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Using the estimates dy*, ... , @i,*, bat, ... , by 
‘ ” a . 
one can compute an “‘expected” value, 7 ;;, for the jth wee! 
of the 7th vear by 


a 


Yigg +a, cos ty+-b7F, sin t 54 


-a* cos rt;+b*t, sinrt;. (7 
Adequacy of the model is seen in the deviations 


di js=Yus Pes (8) 


nh 
| hese 


deviations will occasionally be referred to simply as d, 


of observed values from their “‘expected”’ values. 
without subscripts. Thus the use of this model produces 
+e A . . . 
two additional arrays, d;, and y,, similar in form to the 
original data. Any operations or computations that may 
be performed on the original data, y;;, may be performed 
A : 
on the dy and on the » For example, to the sample 


variance for the jth week, 
] / l Jj 2 
s(y)=- > ¥3,—+( >) y, (9) 
=F | 2 (2) | 


where f is the number of vears in the sample, corresponds 
the variance of the deviations d,, for the jth week. 


1 TY 1/4, 
87 (¢ = > wn ners ° } 
s(d) =i e> iF 7( 3 l, )| (10 


Bliss applied this technique to a 14-year record of 
monthly mean temperatures at New Haven, Conn. His 
findings confirmed those of Craddock: A two-term series 
accounted for more than 97 percent of the observed sum 
of squares. It is well known that the variance of the 
temperature is higher in winter than in summer. Bliss 
found that a simple sine curve fitted to the log-variance 
of monthly mean temperatures adequately described the 
vearly trend of this variance. Concurrently I applied 
the technique to the monthly averages of the diurnal 
temperature range for the same period at New Haven [3]. 
Both the first and second terms of the regression curve 
were significant. The trend of the range was quite 
unexpected with maxima in May and October and 
absolute and relative minima in January and July, 
respectively. The distributions of both the mean and 
the range were normal for all practical purposes, although 
there were very slight indications of svstematic skewness. 


3. METHOD AND DATA 


To carry out the laborious calculations for the application 
of the technique to weekly averages of temperature 
Iaxima, minima, means, and ranges, electronic cata 
processing equipment was necessary. A program was 
written for the IBM 650 with the following functions. 
From cards containing weekly sums or averages of the 
maximum and minimum temperatures, coefficients 
and b,,*, p=1, 2, 3, of the three-term Fourier curve wich 
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best fitted the data (maximum, minimum, mean, or range, 
pending on a coded instruction card) of each individual 


( 
year were computed. 
needed for moments up to the fourth were accumulated. 


Simultaneously sums of powers 


From these “annual” coefficients the expected values 


y, for each week were computed. In addition, the raw 


moments of the deviations, d),,=yij—yy, were accumu- 
lated, yi, being the observed and y,,; the expected values 
for a given week and year (equation (7), (8)). Finally 
the machine calculated an analysis of variance. Later it 
was found useful to write a program to convert raw 
moments to moments about the mean and compute 

9, =h/ke!? and g.=k,/k3 (11) 
where k,, i=2, 3, 4, are Fisher’s k-statistics [8]. I also 


found it advantageous to write a program to carry out a 
harmonic analysis of unreplicated data and to synthesize 
curves from the estimated coefficients. 

The complete set of programs was carried out on the 
data from only two stations, Storrs, Conn., and Keedys- 
ville, Md. both for the vears 1926-56. In addition the 
mean and the range were studied for Uniontown, Pa., 
Eau Claire, Wis. both for 1926-56, and for Easton, Md. 

1896-1926, and 1926-56. 
These last will be referred to as Easton I and Easton IT, 
The climatological vear starting March 1 


for two independent periods, 


respect ively ‘ 
and omitting February 28 and 29 was used throughout. 


4. CENTRAL TENDENCY 


The mean squares from the analysis of variance for the 
harmonic regression fitted to the maximum, the minimum, 
the range, and the mean at Storrs are given in table 1. 
The sums of squares may be easily calculated using the 
The F-tests were applied 
The 30-vear averages of the mean 


degrees of freedom. correct 
according to Bliss [4]. 
| the r ‘or E her with tl er 
and the range for Easton II together with the average 
harmonie curve are shown in figures 1 and 2 respectively. 
Marimum, Minimum and Mean,—As was expected, the 
first (sine curve) term of the Fourier series accounted for 
approximately 90 percent of the total sums of squares for 


from ana'lusis of variance of harmonic 
averages of diurna! temperature maxima, 


Storrs. Conn., 1926 1956 


Tarte 1.—Mean squares 


rression applied to weekly 
and means, 


minima, ranges, 
Mean squares 
I source var 1 I).} 
Maximum Minimum Range Mean 

tween years 29 *113. 649 *73. 631 TSO. 46 
2 st term ) 173223. 360 1193904. 205 
term 2 12. 24 06. 540 
} term 2 *297. 890 47.425 
Scatter about curve 45 149. 283 T32. 424 
t term X year 5S 51.113 TA6. 064 
i term X year 5s 732.114 332. O81 
J term X year AS *29. 750 $29. 283 
Residual 1305 23. 204 21. 105 
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Ficure 1,—30-year average of mean weekly temperature and the 


fitted harmonic curve, Easton, Md., 1896-1926 (designated I, 


and Easton, Md., 1926-1956 (designated II). Week O begins 
March 1. 
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Figure 2.—30-year average of weekly average diurnal temperature 


range and the fitted harmonie curve, Easton II, 1926-1956 


The fit 
third 
larger than one could expect 


and the mean. 


In a few eases the second or 


the minimum, 


appeared excellent. 


the maximum, 


terms were significant (i.e., 
from chance variation under the null hypothesis) but in 
terms of the percentage of sum of squares accounted for 
(<<1 percent) they were trivial. However, the inter- 
actions of terms by vears all were significant and for the 
second and third terms were more important than the 
This indicates that the 


shape of individual vears cannot in general be adequately 


average effect of these terms. 


described by a simple sine curve, despite the good fit to 
the averages. 

An interesting feature was the consistency over a wide 
area of the departures of the average values for the period 
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1926-1956 of the maximum, minimum, and mean from 
the average Fourier curve. The magnitude and direction 
of the deviations for all four eastern stations were almost 
indistinguishable. In addition, at Keedysville and Storrs, 
the pattern of dev iations of the mean was almost perfectly 
reproduced by the corresponding deviations of the maxi- 
mum and the minimum. “The time of vear when the fit 
was least good and when the bulk of the sum of squares for 
scatter arose was the period from mid-December to mid- 
February. The last two weeks of December were con- 
siderably below the fitted curve, January above it, and 
Although differing in detail 
the same general pattern of winter temperatures was 
observed in both 
Easton, Mad. (see fig. 1). 
cent of the fabled “January thaw.” 

of the “thaw”’ is reportedly well defined at January 20-23 
(11) {14|. This does show up clearly in the Easton I 
record but is not visible in the Easton II record or at other 
This casts doubt 


February again below it. 
independent 30-year samples from 
The warm January is reminis- 
However, the date 


stations analyzed for the later period. 
on the reality of this “singularity,” especially in view of 
the occurrence of maximum variability in January. No 
other such similarities in deviations are apparent in both 
records from Easton. 


As shown in figure 2, the range followed the 


Range 
pattern uacovered in the preliminary analysis of the 
New Haven [3]. We observed pro- 


in early- to mid-May and in mid- 


monthly data at 
nounced maxinia 
September with a summer relative minimum consider- 
ably above the winter minimum. The summer dip was 
least pronounced at Keedysville and Eau Claire but was 
clearly present in all records, including both independent 
Easton records. There were few, if any, recognizable 
similarities in deviations from the fitted curve between 
the stations. However, at all stations, the first two terms 
of the Fourier series were highly significant, with the 
third of lesser importance although still significant. The 
interaction of terms by vears was quite uniformly signifi- 
cant This, of course, indicates that there is considerable 
variation between vears in the shape of the yearly course 
of diuraal range. Of considerable interest is close agree- 
ment of the phase of the first term (fundamental sine 
curve) with that of the sun. This apparently reflects the 
close relationship between the diurnal range and = the 
energy input. That this relationship is not overpowering 
is, however, clearly shown by the spring and fall maxima. 
Despite first guesses that these maxima were corollaries 
of clearer skies in both spring and fall, sunshine and 
cloudiness records showed that in fact the spring period 
tended to be cloudier than the summer, although the high 
daytime gains and nighttime losses of radiation due to 
clear skies in the fall remain an acceptable explanation for 
the autumnal maximum. An examination of the daily 
temperature record for Mt. Carmel, Conn. suggested a 
It was observed that there was a 


possible explanation. 
considerable number of days when the daytime tempera- 
ture went very far above the nighttime minimum but 
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In the late 
spring, the annual course of insolation is considerably 


returned to previous levels after sundown. 


ahead of the course of soil temperature, and hence, in 
general, ahead of the course of nighttime air temperatures, 
This makes it possible when conditions are right, for the 
temperature to rise sharply in the daytime and yet return 
at night, as we observed, to about the same level as the 
previous night. This behavior is corroborated by a com- 
parison of the variances (between vear and within year) 
of the maximum and the minimum temperatures. During 
these spring months the variance of the maximum is 
definitely greater than that of the minimum. 


5. DISTRIBUTIONS 


The analysis as discussed to this point has provided an 
efficient method for estimating a central tendency or 
location statistic for the distributions of temperatures 
and temperature ranges. Instead of 52 individual means, 
we have as few as three coefficients (for a sine curve) 
which give the location of the distribution for every week. 
However, for statements of probability, knowledge of the 
shape of the distributions is required. If it ean be shown 
that the distributions are normal, or Gaussian, then the 
mean and the variance completely specify the distribution. 
On the other hand, if the distributions are non-normal, 
higher moments are necessary to specify or to approxi- 
mate them. 
be fully applicable, several assumptions should be met: 


In addition, for the analysis of variance to 


The deviations of each observation from the theoretical 
regression curve for the year should be (i) normally dis- 
tributed, (ii) homoscedastie (i.e., of equal variance), and 
(iil) independent. To examine the distributions of the 
observations and to test the first and second of these 
assumptions, moments were computed and the following 
procedure carried out. The sample variance for each 


week, s(y), JO, . . . 51, (equation (9)) was computed 


from the observed values under consideration, to be 
referred to as y, the maximum, minimum, range, or mean. 
(equation (11)) as 


In addition, I eal- 


[I also computed g,(y) and go(y) 
measures of skewness and kurtosis. 
culated the same statistics s;?(d), g,(d), g2(d) for the dis- 
tributions of the deviations, d, of the observed values 
from the best fitting annual curves (see equation (10) for 
example). This process gave 52 values, one for each week, 
of the sample statistics just mentioned, for each of the 
two related distributions. 

Variance.—The observed variance of the maximum and 
the minimum for Storrs and Keedysville changed smoothly 
over the year, as did the variance of the mean for all 
stations. Corroborating the previously mentioned tend- 
ency, the variance of these variates was consider: bly 
higher during the winter months than during the summer, 
with a reasonably smooth transition between the extremes 
Unfortunately for strict fulfillment of 
the conditions of the analysis of variance, the variance of 
the deviations, d, about the individual curves showed the 
In the case 


(see figs. 3, 4, 5). 


same pattern, although to a reduced degree. 
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FicUrE 3.—Sample variance of weekly average diurnal maximum 


temperature and transform of the sine curve fitted to the log 


variance, Storrs, Conn. 1926-1956. 
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Figure 4.—Sample variance of weekly average diurnal minimum 
perature and transforms of the sine curve (solid line) and the 
-term Fourier curve (dashed line fitted to the log variance, 
Storrs, Conn., 1926-1956 
of the range there were indications that the variance 


followed a double maximum pattern. Hardly visible at 


Storrs, it was more apparent at Keedysville, and still 
Furthermore, it could be dis- 
(The 


observed variance of the range at Uniontown is shown in 


more so at Uniontown. 


cerned in both independent Easton records. 


figure 6. 


Since the trends in the variance were pronounced, we 
summarized the pattern by applying to the variance the 
same techniques of harmonic regression previously used 
on the temperature variates themselves. In order to 
minimize the non-normality of the distribution of the 
variance, it was transformed to its logarithm before the 


analysis was carried out [2]. As Bliss indicated, the an- 
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e 


nual course of the log variance of the mean could be ap- 
proximated by a simple sine curve. The same was true 
for the log variance of the maximum temperature. In 
both cases neither the second nor the third term was 
significantly different from zero. However, the analysis 
of variance indicated that, for the minimum, the higher 
terms were significant. The second term component was 
quite pronounced at Storrs while the third term was im- 
portant at Keedysville, although at both stations the sine 
curve was clearly the dominating feature. The range 
exhibited a markedly different character. The log vari- 
ance showed a significant second term trend, but little 
evidence of any single wave. For all four variates, when 
the deviations, d, about the annual curves were considered, 
the analysis showed that the F values for all terms de- 
creased although the same relative importance of terms 
seemed to be the rule. 

Correlations Between the Of related interest 


to the course of the variance of the individual variates is 


Variates. 
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Figure 5.—Sample variance of weekly average diurnal mean tem- 


perature and transform of the sine curve fitted to the log variance, 
Storrs, Conn., 1926-1956 
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FIGuRE 6. 
perature range and transform of the two-term Fourier curve 
fitted to the log variance, Uniontown, Pa., 1926-1956 


Sample variance of weekly average of diurnal tem- 
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the course of the correlation between them. Because the 
maximum-minimum and the mean-range coordinate sys- 
tems are orthogonal, the correlation between the average 
maximum and the average minimum temperatures con- 
tains all the available information on the various correla- 
tions. Although not computed directly in our original 
program, this correlation is easily calculated from the 
variances of the maximum, minimum, and range as 


YY 22 
*"max 5 


+S a 
range | (12) 


. “mir 
r (max, min) — 


Qe / 2 
<Smax Smin 


This calculation was carried out for the two stations, 
Keedysville and Storrs, for which all these variances were 
available. At both stations, the correlation r was always 
positive, reaching a maximum in winter and a minimum 
Like the variance of the mean, it followed a 
To obtain 


in summer. 
fairly smooth trend between these extremes. 
a quantitative description of its seasonal course we trans- 
formed r to tanh7'r, in order to stabilize its variance 
and minimize the non-normality of its distribution [9] 
and carried out a harmonic regression on each station. 
The variance about a regression curve fitted to 2 derived 
from samples of 30 is theoretically 1/27=0.037037. We 
can insert this as an added line in the analysis of variance 
with infinite degrees of freedom. This can be used in 
place of the ordinary error row in an analysis of variance 
of replicated data to test the adequacy of the fit and the 
reality of the regression. For both Storrs and Keedysville 
the scatter of z about a sine curve, when tested in this 
manner, was not significantly different from the theoreti- 
cal variance (P20.2, F’<1.2), indicating there was no 
removable systematic variation in the residuals about the 
curve. The sine curve was highly significant (P<0.01 
Thus we may conclude that the transform, 
lation between the weekly averages of maximum and 
minimum temperatures follows a sine curve. The only 
apparent regular departure from the sine curve is in April 
and early May when, at Keedysville, there were eight 
successive weeks when the observed correlation was less 


>, of the corre- 


than the “expected” correlation, calculated from the sine 
curve. At Storrs a similar feature was observed for the 
sume eight weeks, with the exception of one week in 
which the observed correlation slightly exceeded the ‘“‘ex- 
pected.”’ In addition to this, there were other similarities 
between the two stations in the observed correlations. 
This is a further reflection of the previously noted similari- 
ties in the temperature records at quite widely separated 
stations (340 miles). 

One might hope that, by using the “expected” values 
of r(max, min), max, and Spin, Computed from the sine 
curves fitted to their transforms, one could recover the 
eurve deseribing the course of the variance of the range 


by the formula, 


° > ° a, 
a — N x N 
“min <"Smax min- 


g* Re 
range max 
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I hoped in this way to obtain a purely mathematical 
explanation for the behavior of s*ange based on simple 
assumptions concerning the behavior of the other second 
moments. The results, however, are not convincing. 
For Storrs, the curve of s*range computed in this way 
is far too flat, although its only notable feature, a peak 
at week 10, does reflect a similar peak in the observed 
variance. At Keedysville, where the observed variance of 
the range displayed a pronounced double maximum pat- 
tern, I found only a single peak at the time of one ob- 
served maximum, with a point of inflection at the less 
important of the two minima. 

Homoscedasticity and Independence of Residuals.—The 
above analysis has made one thing clear: Our assumption 
of hooscemdasticity is not fulfilled. Although the mathe- 
matical model postulated implies some sort of a periodic 
form for the variance of the actual values, y, the variance of 
the deviations, 7, from the annual curves should be con- 
stant. In the case of the average diurnal range this is 
nearly true. However, the other three variables are far 
from homoscedastic. In any case, the estimates of the 
regression coefficients are still unbiased estimates of the 
population values. It is difficult to assess the effect of 
this upon the tests in the analysis of variance. According 
to Cochran [5], the F-test for non-regressive designs is 
sufficiently robust as to be not misleading, although the 
more general case is not covered. Since, for the mean and 
the two extremes, the sine curve predominates, we can 
clearly accept the adequacy of the overall fit, although the 
tests cannot be considered exact. Since the variability 
of the variance of the range is far less prominent and regu- 
lar, the tests should be less affected by heteroscedasticity. 

Of equal or greater importance, perhaps, although more 
difficult to assess, is the possibility of dependent residuals 
Direct tests for this are available [1] but laborious to 
apply and, as they stand, are not readily applicable to our 
computing scheme. Bliss [4] states that the technique of 
fitting separate curves to each vear will have the effect of 
removing serial correlation between weeks, and leaving 
substantially independent residuals. From inspection of 
individual vears, the observed residuals appear to be ran- 
dom. On the assumption of independent residuals, it 
follows that the moments for each week, as computed 
from the residuals, are independent and thus also ¢ 
and g, computed from these moments are independent. 
Serial correlation in these should tend to inflate the 
significance of trends in these statistics. Thus we would 
expect that if the residuals are independent, the F values 
in tests for the reality of regression curves fitted to the 
annual course of g; and g2 should be lower for the distribu- 
tion of the residuals than for the distribution of the raw 
This, as we shall see, was observed and can 
However, no 


observations. 
be considered evidence of independence. 
conclusive test has been found to clarify this point. 

Tests for Normality.—_We now examine the assumption 
of normal distributions for our variates. Being able to 


work with normal distributions is desirable for three 

















Sr TEMBER 1961 
re:sons. Firstly, such normality is assumed in our use 
of variance ratio tests in the analysis of variance, although 
it has been shown [5] that departures from normality, if 
not extreme, have little effect on the F test. Secondly, 
the normal distribution is known and easy to apply. 
Thirdly, if a distribution is normal all information about 
the values of the parameters defining it is subsumed in 
the sample mean and variance [7]. Furthermore, the 
powerful central limit theorem, applicable because we are 
considering averages of several observations in our 
distributions, tells us that the distributions approach 
normality. 

A further consideration is also important: the use to 
derived from the 
If extreme accuracy at all levels 
of probability is required, for instance, we must be 


which the probability statements 


distribution will be put. 


extremely stringent regarding our tests of the distribution. 
On the other hand, in practical climatological applications 
we are not interested in, for example, the difference 
Even if 


tests show that almost certainly there is some departure 


between a once-in-20 and a once-in-25 event. 


from the normal curve, we may accept a normal approxi- 
mation if the errors in probability estimated from our 
approximation do not impair the usefulness of the 
estimates, 

As measures of departures from normality I chose, as 
k/ko'/? and g.=k,/k2 where 
Departure of these from zero 
Both 


asymptotically normally distributed with mean zero and 


previously mentioned, 4g, 
k, are Fisher's & statistics. 
is indicative of non-normality. Statistics are 
variance depending, in the null case, only on the sample 
size. Skewness or assymmetry is measured by g; while g» 
measures kurtosis. For ease in machine computation, go 
was chosen in preference to 
Sz,—3] 


a 


The g,’s and 
First, compare the 


recommended by Geary and Pearson [10]. 
g's can be tested in two ways. 
observed distributions of the g’s with those expected for 
samples drawn from independent normal populations. 
The ordinary significance test using a standard error is 
Second, examine the vearly course of g; or 
meaningful improbably regular 
pattern is as clear evidence of the presence of non-zero 


of this type. 
g, for patterns. An 
skewness or kurtosis as are high values of g, and qo. 

The chief obstacle to the first type of test is our igno- 
rance of the exact distributions of g, and g) under the 
null hypothesis. It is known that for samples as small as 
30, such as concern us here, the distribution of g, is 
not far from normal while that of g is strongly positively 
skewed. When the cumulative sample frequencies for 
the y’s at Keedysville were plotted on probit paper, the 
curves reflected this expectation. The distributions of 
the y,’s were quite linear indicating approximate normality 
while the go’s exhibited a concave upward curve character- 
istie of positively skewed distributions. 
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TABLE 2.—Number of exceedences of the Ist, 5th, 95th, and 99th 
percentiles in 52 values of gi(y) and g,(d) calculated from weekly 
averages of diurnal temperature maxima, minima, ranges, and 
means, Keedysville, Md., 1926-1956 


ay ald 
>upper | <lower Total >upper <lower Total 
5 percent | 5 percent 5 percent | 5 percent 
Maximum 4 l 5, 1 l 2 
Minimum 2 3 5 2 2 4 
Range 3 l { 2 3 F, 
Mean 1 6 | 0 l 
> upper <lower Total >upper <lower Total 
1 percent | 1 percent I percent | 1 percent 
Maximum 1 0 1 0 0 0 
Minimum 2 2 4 0 0 0 
Range 0 0 0 0 0 0 
Mean l 0 l 0 0 0 


Although the exact distribution of g, is unknown, 
Geary and Pearson [10] have given approximate extreme 
percentage points. A comparison of the observed fre- 
quencies of g; beyond these points seemed the optimal 
procedure. Since these percentage points are not avail- 
able for g, for sample size less than 100, this comparison 
was not possible for g. In table 2 are given the number 
of times the upper and lower 5 percent and 1 percent 
levels for g, are exceeded for the 52 weekly distributions 


We see that 


the agreement is as good as could be expected for the 


of the observations, Y, and the residuals, d. 


distribution of the deviations, d, about the annual curves 
although there are improbably few (i.e., no) values beyond 
the 1 percent point. However, in the distribution of the 
observations, y, themselves, the most striking feature ts 
the appearance of four g; (y)’s beyond the 1 percent point 
among the 52 g,’s derived from the minimum temper- 
Further, the upper 5 percent point of g,(y) for 
Thus 
a certain degree of non-normality is indicated in all the 
There is, however, 


ature. 
the maximum and the mean is exceeded too often. 


original variates except the range. 
no indication that this is true for the deviation from the 
annual curves, d. This last result is important since it is 
assumed in our mathematical model. 

Clearly, the above tests lose much of their validity if 
there is appreciable serial correlation between weekly 
values since this makes the g,’s serially correlated. This 
affects the shape of the observed distribution of the 52 
values of g,(y) for each variate. Thus the second type 
of test referred to above may be more applicable. This 
can best be done by fitting a regression curve on time to 
the computed statistics to uncover any significant pattern 
over the year. If the regression accounts for a significant 
part of the variation, it indicates a real departure from 
the expected values of zero. Because of the general 
robustness of the F-test, the departure of the distributions 
of the g,’s and g’s from normality should not invalidate 
tests of significance of a regression curve. Accordingly, 
as previously mentioned, a three term harmonic curve was 
fitted to each of all available sets of g,’s and g,’s, with 
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somewhat mixed results. For the distribution of d there 
are no more F values significant at the 5 percent level 
than one would expect, with the exception of g, for the 
maximum temperature. For this, the second term for 
Storrs is significant and at Keedysville both the second 
and third terms are significant, although in both cases 
they are quite small. One may conclude, however, that 
there was no appreciable systematic skewness or kurtosis 
in the residual variation, whether one considers the 
mean, range, Maximum, or minimum. 

When we examine the observations y, there is a con- 
sistent pattern among the g,(y)’s. For four out of the 
six records analyzed, including Easton II, the second 
term of the curve fitted to the course of g,(y) for the mean 
temperature was significantly different from zero. In 
the fifth record, Easton I, 
term was more than twice as large as any other, although 


the amplitude of the second 
not quite significant. Eau Claire, on the other hand, 
demonstrated a strong first term and small second and 
third terms. This discrepancy may possibly be traced to 
climatic differences between the Atlantic coastal and the 
Lake States. 


even more apparent in g,(y) for the minimum temper- 


The same double maximum pattern was 
atures. At both stations analyzed, Storrs and Keedys- 
ville, the second term was highly significant, an almost 
identical pattern emerging, There appeared to be a 


tendency toward positive skewness in the fall and 
especially in the spring and negative skewness in the 
summer and especially in the winter. 

We conclude that the distributions of the deviations 
d from annual curves are Gaussian for all four variates, 
minimum and the 


while the observations y of the 


mean have a skewness that changes seasonally. Since 


this pattern of skewness is most pronounced in_ the 
minimum and since the mean, as defined herein, is in 
part derived from the minimum, the primary pattern 
of non-normality is likely to be in the distribution of the 
minimum temperatures. 

Although the tests discussed so far have not indicated 
any real departures from normality in the distributions 
of the range and the maximum, one further test placed 
the normality of the distributions of these 
By the central limit theorem, the dis- 
tribution of averages of independent g,’s or g:’s should 
be approximately normally distributed with zero expec- 


tation and variance 1/n times the variance of a single 


in doubt 


variates, too. 


value. 
of g;,, i=1, 2, we can treat it as a normal deviate with 


Thus if we consider the mean J, of the 52 values 


These then can be compared with 
Any 
means outside, say, the 5 percent level would indicate 


variance 1/52 V(q,). 
the percentage points of the normal distribution. 
a significant average departure from normality. Since 
we have no a priori knowledge of the direction in which 
deviations from the null hypothesis should occur, the 


When 


proper test to use is the two-tailed comparison. 


this test was carried out, no average departures from 
zero were found for either the g,(d)’s or the g.(d)’s of 
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the distributions of the deviations, d, for any of the four 
Neither were there significant non- 
However, both 


variates studied. 
zero average departures in the go(y)’s. 
the maximum and the range displayed definite signs 
of positive average skewness among the y’s. At both 
Keedysville and Storrs 9,(y) surpassed the upper 5 per- 
cent level and, in the case of Storrs, the upper 1 percent 
level. For the range, all except one Storrs, 
with 9,(y)=—0.0041, showed some positive skewness, 
with two records, Easton I and Uniontown, surpassing 


station, 


the 1 percent and 5 percent levels respectively. 

Conclusions Regarding Distributions._-First, the varia- 
tion of all four variates about the annual harmonic curves 
is Gaussian. This is important from a theoretical point 
of view, since it increases our faith in the underyling 
model. However, because we cannot predict the shape 
of the yearly course of the variate, this knowledge is 
of little use from a practical point of view. Second, we 
have quite clear indications of skewness but not kurtosis 
in the distributions of all four variates. In the case of 
the range and the maximum temperature, the average 
skewness is consistently and, in some cases, significantly 
On the other hand, although they display 
no average skewness, the distributions of the mean and 


positive. 


the minimum have skewness that varies systematically 
over the vear. Thus, for completely exact climatological 
statements of probability, normal assumptions will not 
be sufficient. 

In theory one should try to find the exact nature of 
the distribution for each of the four variates. In_ the 
present instance, this could be more misleading than 
any assumptions of normality. First, there are insuf- 
ficient data to establish any distribution as being correct 
without doubt. Second, the variation in skewness, at 
least for the mean and minimum temperatures, suggests 
that the 
This would greatly increase the difficulty of applying 
any distribution. An alternative which is more approach- 
able is the use of the third and fourth moments in fitting 
ad hoe distributions to the variates using either one or 


“correct” distribution may vary seasonally 


several of the Pearson curves or the Edgeworth approxi- 
mation. Since we have no evidence of the departure 
of the fourth moment from the value expected for the 
normal distribution, the third moment, as reflected in 
gi(y), should Further, due to the large 
sampling variation in g,(y) resulting from our use of small 


be sufficient. 


samples, our best estimate of the population skewnesses 
should be the values calculated from the fitted periodic 
regression curve, or in the case of the maximum or the 
range, from the yearly average. Even if such an approxi- 
mation does not provide an exact fit, the degree of change 
from the normal approximation should indicate the 
magnitude of error. 

Tebles of the cumulative distribution functions of the 
first Edgeworth approximation for different values of 7 
the population measure of skewness, are available [13] 
Referring to these I found that for small departures from 
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s)mmetry (|g;|<0.2) use of this approximation does not 
markedly affect probability points given by the normal 
distribution. Assuming a standard deviation, o, of 7° F. 
(near the maximum for any variate) at g,=+0.2 the 
upper and lower 1 percent points are displaced approxi- 
Thus 


it seems safe to say that use of the normal distribution 


mately 1° F. with lesser changes nearer the median. 


for the maximum and the range does not introduce any 
serious errors since the greatest average g,(y) (Storrs 
maximum) is 0.1643. For the minimum temperature at 
Keedysville which displayed the strongest systematic de- 
partures from zero of g,(y), the “expected”’ value of g;(y) 
ranged from about —0.48 to +0.52. It is true that the 
1 percent level is quite distorted by this amount of asym- 
However, the 5 percent level is only about a 
7° F.). 


tion which by the normal distribution would be surpassed 


metry. 
degree off (again letting o Conversely, a devia- 
5 times in 100 years would, under these circumstances, 
Similarly the esti- 
mated 10 percent point would be exceeded about 11 times 
every 100 years. 


occur about 6.5 times in 100 years. 


The Edgeworth approximation should 
be fairly accurate when the shape of a distribution is not 
far removed from the normal form. Since, as mentioned 
earlier, the central limit theorem assures that the distri- 
butions approach normality, the use of this approxima- 
tion in the present case should be appropriate. Hence, 


these figures show that the use of the normal curve 
will be satisfactory on all levels of probability except 
beyond the 5 and 95 percentiles where one should prob- 


ably apply some extreme value distributions. 


6. PROBABILITY ESTIMATION 


The above considerations on the distribution of tem- 
perature variates lead to a satisfactory method of esti- 
mating probabilities. Since we may assume for our 
purposes that the variates are normally distributed, spe- 
cification of the means and variances completely deter- 
mines the distribution. Our harmonic curves fitted to 
the variates provide the first of these parameters, while 
the curves fitted to the log variance provide the second. 
In the case of the maximum and mean temperatures both 
curves are essentially sine curves, with higher terms, 
even when satistically significant, changing the esti- 
mates little. 
is adequately described by a sine curve. 


The course of the minimum temperature 
The higher 
terms of the curve fitted to the variance do not seem 
to make an appreciable difference. For example, for 
Storrs the 5 percentile and the 1 percentile are changed 
a maximum of 0.8° and 1.2° F. by adding a second term. 
Similarly, the probability assigned to a given departure 
from the mean seems to be changed by at most 0.02 to 
0.03. Since this is less than the level of accuracy ordi- 
Hence 
mean, the 
estimation of probabilities for any week of the year re- 
I should point out 


narily desired, the sine curve alone can be used. 


for the maximum, the minimum, and the 


duces to the fitting of six constants. 
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that although a sine curve is adequate to describe the 
courses of these temperatures at the stations studied, the 
possibility is not excluded that at some locations higher 
Each term can, however, be 
specified by only two additional constants. 


terms may be necessary. 


For the range, at least two terms, i.e., five parameters, 
should be used for estimating the central tendency. 
Since at some stations the second term in the curve de- 
scribing the course of log s*,ange is of prime importance, 
five constants seem to be necessary to estimate the vari- 
ance. However, in cases where there is no discernible 
pattern in the variance, it is probably best to use a single 
average variance derived from the row for scatter about 
the average curve. Thus for Storrs one would use 15.7 
(° F.)? as the variance of the range (table 1, row 5). 

Although the above discussion puts forward a relatively 
simple technique for estimating the mean and variance of 
temperatures from a small number of constants, there is 
a further simplification that may produce equivalent re- 
I noticed that 
not only does the log variance of the mean temperature 


sults, at least for the mean temperature. 


follow a simple sine curve but also the time of the maxi- 
mum variance is about 180°, or 6 months, out of phase 
with the temperature. For example, the yearly maximum 
of the mean temperature at Keedysville, computed from 
the best fitting simple sine curve, falls 137°27’ or 139.4 
days after March 4 (the midpoint of week zero) while the 
‘cm 313°48’ or 318.2 days after 
March 4 and 176°21’ or 178.8 days after the maximum 
temperature. 


maximum of log s is at 
This suggests that the log variance and the 
mean may be linearly related with a negative slope. 
When 7 was plotted against log s, this was seen to be 
In fact, the logarithm so little 
changes the relationships among the standard deviations 


substantially the case. 


that an almost equally good fit to a linear relationship is 
obtained if one plots 7 against s. If one assumes that the 
phase difference between the log variance and the mean 
temperature is exactly 180°, which is very close to the 
observed difference, then one can rather easily compute 
the slope of the regression of log s on 7 from a knowledge 
of the amplitudes of the respective sine curves. 


If 


Y=ao+a, sin t+, cos t 


d)+ A, cos (t ¢)) (14) 


where A, is the semi-amplitude and @, is the time in 
degrees of the maximum of the mean as measured from 
the start of the climatological vear, and if 


log s=a)+a; sin t+); cos t 


a,+B, cos (t—¢— 180°), (15) 


then 


™ log s ay, 
cos (t—¢;) cos (t—¢d,— 180") RB ° 
>, 
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Thus 


BB, , B, 
y+ {FoF Go. (16) 


logs=— 


Hence we can give the “regression” of log s on y from a 
knowledge of the sine curves fitted. To test the accuracy 
of this method, both the ordinary least squares line and 
the line computed from equation (16) were calculated 
for the mean temperatures at Storrs and at Easton I and 
Easton II and are given in table 3 together with the 
regression line calculated from equation (16) for the other 
stations studied. We note a remarkable homogeneity 
among the slopes of the regression lines for the three 
records which cover two Maryland stations, with the slope 
at Storrs quite close. However, when one crosses the 
Alleghenies to Uniontown, the relation changes radically, 
Even if the 


observed interstation homogeneity in the three eastern- 


and is even more changed at Eau Claire. 


most stations is fortuitous, the stability of the regression 
lines between two independent records at Easton shows 
that it may be possible to estimate temperature proba- 
bilities for any week of the vear from only five parameters, 
the overall mean, the phase angle, and the amplitude of the 
temperature curve, and the amplitude and mean of the 
log variance. This however, need 


further study, especially with respect to the geographical 


relationship does, 


stability of the ‘‘regression”’ of log s (or s) on y. 

One point in the above discussion may need some expla- 
nation. It is a well-known fact that in drawing samples 
from a single normal population, s*, the estimate of the 
population variance, and 7, the estimate of its mean, are 
How then do we have such a 
the log 


independently distributed. 


marked correlation between the means and 
variances, when we have seen that the distributions are 
not far from normal? The question is, however, meaning- 
less since, although we have drawn 52 samples from normal 
or near normal populations, these populations have not 
been the same. There is neither a priori nor a posteriori 
evidence that the population of mean temperatures for 
week one is the same (i.e., has identical characteristics) 
as the population for any other week. In fact, our in- 
vestigation has demonstrated that these distributions 
are definitely different. Thus one should not expect the 
independence between s* and 7 which would result if the 
populations were the same. A relationship between the 


true values (not merely the estimates) of the averages and 


TARLI 3. Regre sston equations of log 8 on Y computed by least 
squares and by equation (A) 
Station Least squares line Equation (¢ 
Easton I , log s =1.020113—0.00644479 log s = 1.02479 =0.0065439 
Easton II log s =1,018815—0.006057 7 log s=1.02260—0,0061289 
Storrs log « =0.926768 —0.0052399 log « =0.92902 —0.005277 9 
Keedysville log s =1.02814 —0.006067 7 
Uniontow! log «=1.18340—0.0083119 
Eau Claire log s =0.89526 —0 002067 7 


Not com puted 
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variances of the mean temperatures, say, is clearly inci- 
cated by the data at hand. Its reality can be indicated 
statistically but its cause must be physical. 
vein the standard tests for equality of variance, say, 
between two different 
advance of such a test, we have a well established result 


In a similar 


weeks are quite useless. [np 


that the variances are, in general, unequal. 
7. PROBLEMS IN FURTHER APPLICATIONS 


To this point we have spoken only about the probabil- 
ities of average temperatures or average diurnal ranges 
for a given week, without reference to the probabilities 
for derived quantities, or averages over longer periods 
If the mathematical model 
were completely correct, one could, theoretically, compute 


such as months or seasons. 


quite complicated probabilities concerning degree days, 
frost dates, etc., from a knowledge of the joint distribution 
of the coefficients of the curves which describe each vear, 
This knowledge is also needed for efficient use of the 
parameters in relating them, as being in some sense a 
description of a yearly temperature regime, with other 
Before 


proceeding further, we should inquire to what degree the 


elements, climatological, agricultural, or physical. 


estimated coefficients for a given vear do contain the 
salient features of that vear. To this end, I tabulated 
the published departures of the monthly mean tempera- 
tures from the long-term normal for the regions or States 
in which our stations are located. Years were then picked 
by eve which had various characteristies, e.g., warmer than 
normal spring combined with colder than normal winter, 
ete. Then, using the best fitting harmonic curves for 
these vears for the stations involved, the “expected” 
temperature, y, for each week was computed and com- 
pared against the 30-vear average curve. In every case 
the unusual pattern for which that vear had been chosen 
was observed in the fitted harmonic curve. Although 
they lack fine detail, the curves for each vear definitely 
reflect the course of the temperature for that vear. 
Determination of the distribution of the true coefficients 
is far more difficult since they cannot be observed directly. 
A further 
complication is the previously discussed dependence of 
Thus 


a full discussion of this problem would involve determin- 


We can examine only estimated coefficients. 
the coefficients upon our arbitrary choice of origin. 
ing all relationships among the coefficients which are not 


However, despite the 
valuable to examine 


changed by any choice of origin. 
above objections, I considered it 
more carefully the estimates of the coefficients for the 
particular choice of origin we used. To provide a mini- 
mal test of normality in the marginal distributions of th 
separate coefficients, the ranked estimated coefficients for 
the curves fitted to the mean temperature were plotted 
against rankits (expected order statistics for the normal 
distribution [9]) for each of the two Easton samples 
Departure from normality of the estimates should bi 
reflected in non-linearity in the plotted points. In every 
case where any noticeable non-linearity occurred in on 
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Thus I concluded 
that there was no evidence for rejecting the hypothesis of 


rccord, it was not visible in the other. 


normality among the estimates of the coefficients. To 
the degree that non-normality in the distribution of the 
actual coefficients would be reflected in non-normality of 
the estimates, this indicates normality in the marginal 
Hence, they are 
probably distributed in a multivariate normal distribu- 
tion. Further, we computed covariance and correlation 


distributions of the actual coefficients. 


matrices for the seven estimated coefficients for each of 
the two Easton samples, and then combined them to get an 
estimated covariance matrix based on a 60-year sample. It 
can easily be shown that for our choice of origin, the winter 
maximum of the variance of the mean temperature implies 
that there should be positive correlation between a) and 


a, and negative correlation between dp and 6,. Hence 


any test for the significance of these observed correlations 
should be single-tailed, and we should also compute 
We found that there 


were several correlations of statistical significance, i.e., 


the partial correlation r(a;, ;\a9). 


larger than one could expect to occur by chance alone. 
However, their practical significance is negligible since 
they are all less than 0.5 in absolute value. The physi- 
eal or climatological meaning of such correlations is 
not clear, especially in view of the fact that any one may 
be reduced to zero by some choice of an origin in fitting 
the polynomials. 

There is one further problem which would limit the 
application of this method, even were the above problems 
solved. There appears to be considerable difficulty in 
making a meaningful correspondence between the param- 
eter space and the space of “‘vear types.”’ In other words, 
there is no clear way to pick out all those sets of coeffi- 
cients @, d; . .. , d3, 6), . . . , bs, which, when used in 
a harmonic curve, would produce a ‘‘warm spring.” 
Unless one can accomplish this, knowledge of the distri- 
bution of the coefficients cannot be fully applied. 

There are many applications, however, of the average 
curve computed for a location, in which the form of the 
distribution of the vearly coefficients is of minor im- 
portance. Prescott [12] used maps of the phase and 
amplitude of sine curves fitted to monthly temperatures 
for stations in Australia to locate suitable homoclimes for 
Except for the difficulties attached to solving 
polynomial equations involving trigonometric terms, ex- 


hew crops. 


pression of the vearly course of temperature by a periodic 
curve provides a means of locating the time when the 
average temperature is at a Maximum or a minimum, 
If 


Y¥=a)+a, cos t+, sin t+a, cos 2t+6, sin 24+ ..., 
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the relative maxima and minima occur at those values for 
t for which 


—a, sin t+b, cos t—2a, sin 2t+2h. cos 2t— . . .=0. 
Similarly, the day on which the average maximum, 


minimum, or 
value can be computed from the coefficients of the aver- 


mean temperature passes some critical 
age curve. In the case of the mean temperature, of 
course, this type of problem is generally simple since the 
temperature can be considered to follow a sine curve. 
However, for more complex curves, a process of successive 
approximation seems to be necessary. 
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THE WEATHER AND CIRCULATION OF JUNE 1961 


A Hot, Dry Month in the West 


1. RECENT MONTHLY TRENDS 


From March 1961 thru June 1961 there was a noticeable 
trend in the monthly weather anomalies and circulation 
over the United States. 
of the United States tended to become progressively cooler 


Temperatures in the eastern half 


than normal while temperatures in the western half of the 
country continued a warming trend, At the same time, 
there was a gradual shift in the area of heavier than 
normal precipitation from northern and eastern portions 
to the Southern States and a noticeable tendency toward 
drought conditions in parts of the West. 

A corresponding trend was apparent in the circulation 
in the mid-troposphere. Comparison of successive month- 
ly mean 700-mb. charts for March through June shows 
the gradual displacement westward of a long-wave trough 
from the western Atlantic at the beginning of the period 
to the Ohio Valley by June. There were other areas in 
retrograde 


which motion of 


occurred, but in none wes it so pronounced or so important 


apparent major systems 


to the weather in the United States. 


2. MONTHLY MEAN CIRCULATION 


The planetary wave pattern at 700 mb. in June 1961 
consisted of well-defined trough and ridge systems 
While this 


is shown by the height contours alone, it is better indicated 


(fig. 1) 
from the western Pacific eastward to Eurasia. 


by the height departures from normal (dotted lines in 
fig. 1). 

In the Atlantic the circulation was characterized by 
westerlies that averaged considerably faster than normal 
in a jet maximum that extended from Newfoundland to 
the United Kingdom. The strong westerlies resulted from 
the moderately deep Icelandic Low in phase with a strong 
ridge in the Atlantic. The primary Icelandic Low was 
310 feet deeper than normal and was displaced consider- 
ably eastward to the Norwegian Sea from its normal 


position near Baffin Island [1]. 

The Aleutian Low was similarly displaced. 
position in June {1} is in the Bering Sea, and it is usually 
accompanied by a broad cyclonic flow in the central 
Pacifie north of 30° N. 
that the deepening trough in eastern Asia and amplifica- 


Its normal 


The observed flow (fig. 1) shows 
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tion in mid-Pacific contributed to the anomalous position 
of the Aleutian Low in the Gulf of Alaska. Deepening 
of the trough in the Gulf of Alaska was probably assisted 
by local baroclinic processes, suggested by the distribution 
of monthly mean 1000—700-mb. thickness (fig. 2), as well 
as by barotropic transport of vorticity from upstream. 

Growth of the ridge in the central Pacific was a great 
change from May (see fig. 1 of [2]) when a trough was 
located in this area. Heights at 700 mb. increased as 
much as 400 ft. from May to June in the western Aleutians 
as the trough was replaced by a strong ridge in June. 

Over North America the Canadian portion of the ridge 
in the West was somewhat farther east and generally 100 
feet or more stronger than normal. Its companion trough 
downstream was west of its normal position and a little 
deeper than normal. This shortening of the half wave- 
length between the ridge and the trough resulted in a 
vigorous northerly component of the height anomaly field 
from northeastern Canada to Mexico. 

There were pronounced changes in the circulation at 
700 mb. from the first half of June (fig. 3A) to the last 
half of the month (fig. 3B). Elements of the circulation 
over and near North America for June 1-15 included a 
continental ridge and two major troughs, one in the Gulf 
of Alaska and the other in eastern Canada, oriented in a 
United States 


modified omega pattern. Flow over the 


was quite weak except along the northern border. These 
features are shown clearly by the contours of height and 
height anomaly in figure 3A. 

From the first half to the last half of the month the 
circulation was characterized by retrogression and general 
flattening of the pattern at higher latitudes over North 
America (fig. 3). This was associated with weakening of 
the Low in the Gulf of Alaska. At the same time heights 
increased in western United States as much as 200 ft. 
over the Great Basin. Meanwhile, the trough in eastern 
North America retrograded to the Ohio Valley and joined 
with the trough that moved eastward from Baja California. 

In summary, June was a month of relatively low index 
circulation over the United States. This was associated 
not with blocking as was May, but with a strong meridion- 
al flow that was responsible for unseasonable extremes of 
temperature in many parts of the country. 
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Figure 1. 


Mean 700-mb. contours (solid) in tens of feet, and departures from normal (dotted) in 50-ft. intervals, for June 1961. The 


meridional nature of the height anomaly field was indicative of the low index condition over North America and the Pacific 


3. WEATHER RELATED TO THE CIRCULATION 
TEMPERATURE 

Hot, dry weather in the West was associated with many 
hew temperature records principally in Montana, Cali- 
fornia, and while unusually cool conditions 
prevailed over much of the East. 

‘Temperatures in June (fig. 4) averaged about 2° to 4° F. 
cooler than normal from the Southern Plains to the Mid- 
Atlantic States and from the Lower Great Lakes to the 


Oregon, 


Gulf of Mexico. below 
normal heights and cyclonically curved contours on the 
monthly mean 700-mb. chart (fig. 1). The strong ridge 
in western North America steered numerous anticyclones 
southward from northern Canada (Chart IX of [3]) to 


This area was dominated by 


the Northern Plains then southeastward across the Middle 
Atlantic States in a slightly confluent zone. A secondary 
track of Highs extended from western Hudson Bay to 


southern Quebec. In both instances daily anticyclones 
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Figure 2.— Mean thickness (1000-700 mb.) departure from normal 
isopleth interval 50 ft.) for June 1961. Cool air masses 


dominated the East and subsiding air warmed the West. 


bypassed the Great Lakes, a very atypical occurrence in 
June [4]. 

The average thickness from 700 to 1000 mb. (fig. 2) 
was below normal along an axis from Labrador to the 
Gulf of Mexico, generally coincident with the area of 
below normal surface temperatures (fig. 4). Few records 
were established in the East beyond several new daily 
At Norfolk, Va., a new all-time 
was recorded on the 17th; 


minimum temperatures, 
minimum for June of 48° F. 
monthly mean temperatures were reported 
(68.9° F.) and Vicksburg, Miss. 
F.), and the coldest June since 1903 was reported 


record low 
at Louisville, Ky. 
(75.3 
at Shreveport, La. 

In contrast to the East, temperatures were exceptionally 
10° F. 
This region was located under the 


warm in the West. A large area was 6° warmer 
than normal (fig. 4). 
upper level ridge where subsidence was widespread and 
where heights exceeded normal by as much as 170 ft. 
(fig. 1). Many areas in the West reported heat waves 


of varving duration, but one of the most intense was at 


TABLE 1 Vew monthly maximum te m pe rature records in June 19461. 


Maximum 


Locatior temperature 

F 
San Francisco, Calif 106) 
Bakersfield, Calif 113 
Medford, Oreg 109 
Pendleton, Oreg--. 108 
Winslow, Ariz__. 105 
Salt Lake City, Utah 104 
Yakima, Wash 103 
Rapid City, 8. Dak 106 
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FiGurRE 3.— Mean 700-mb. contours (solid) and departures from 
normal (dotted), both in tens of feet for (A) June 1-15, 1961, 
and (B) June 16-30, 1961. 


East and intensification of the ridge in the West were important 


Retrogression of the trough in the 


factors during June. 
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Figure 4.— Departure of average temperature from normal (°F 


Many new temperature records were cstab- 


for June 1961. 
lished, especially in the West. 


(From [5].) 
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Figure 5 Departure of average temperature from normal F. 


for (A) June 1-15, 1961, and (B) June 16-30, 1961. Note that 


B most closely resembles figure s 


Nev., where the daily maximum temperature 
F. from June 17 to 26. Most 
Other 
spectacular records include new maximum temperatures 


Las Vegas, 
ranged from 110° to 113 
ities in Montana had the warmest June of record. 


for June as shown in table 1. 

Above normal temperatures were not confined to the 
United States. The core of warm thickness and above 
normal heights extended northward to the Beaufort Sea. 
Meteorological Branch, 
temperatures 


Reports from the Canadian 
Department of Transport, indicate that 
exceeded normal by 6° to 8° F. from southern Alberta 
to Great Slave Lake and by more than 4° F. 


to 70° N. 


Temperature anomalies discussed above can be examined 


from there 


less grossly by considering the 15-day averages (fig. 5). 
The large changes in the circulation (fig. 3) from the first 
half of June to the last half are clearly reflected in these 
temperature changes. Temperatures for June 1-15 (fig. 
5A) were above normal over much of the United States. 
Maximum departures exceeded 8° F. in Montana and 
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Figure 6.—Percentage of normal precipitation for June 1961 


Only the South and parts of the East received more than 


normal rainfall. (From [5]. 


North Dakota, the which positive height 


departures were largest (fig. 3A). 


region in 


During the last half of June (fig. 5B) large-scale cooling 
in the East caused temperatures to decrease 6° to 10° F. 
from Texas to the Atlantic coast. In the West, influenced 
by the anticyclogenesis aloft (fig. 3B), the area of warmth 
expanded to the southern border and intensified by about 
5° F. in most States 

Widespread North 
occurred when deepening in the Gulf of Alaska reinforced 


warmth over western America 
the ridge that normally lies along the Canadian Rockies. 
Consequently, there was little penetration of western 
North America by Pacific air masses as the average jet 
axis swept northward off the west coast instead of east- 
All aspects 


of the mean circulation that supported the heat in the 


ward along its usual path in middle latitudes. 


West also favored a relatively dry month in this area 


DDLOCIDIT ATION 
PRECIPITATION 


Quite dry conditions prevailed over most of the western 
half of the 
generally heavy in most of the South and in portions of 
the East. 

Rainfall in the West was effectively inhibited by the 
700-mb. flow (fig. 1 in the West and 
its associated ridge were sufficiently strong to present a 


country (fig. 6), while precipitation was 


The ant ievclone 


barrier to the transport of moisture from the Pacific. In 
addition, the desiccating nature of antievelonic conditions 
and the predominantly downslope direction of the anom- 
alous flow contributed to the dry weather. There were 
few monthly records established in the West, but rainfall 
was appreciably below normal and contributed to pro- 
longing a period of light precipitation in that region 

The average 700-mb. flow (fig. 1) was conducive to 
heavier than normal rainfall along and eastward from 
the trough in the central States. Features of the 700-mb. 
chart favoring heavy precipitation in Texas included the 
below normal heights, and easterly 


trailing trough, 
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Figure 7.-—Number of days with fronts in equal area (66,000 
n. mi) quadrilaterals for June 1961. Heavy precipitation 
occurred along the axis of maximum frontal frequency in the 
South; in the North fronts were ‘‘drier’’ due to virtual exclu- 
sion of moisture from Gulf of Mexico by strong northerly flow. 

anomalous wind components. In addition, temperatures 

were fairly cool (figs. 2 and 4) from the surface to 700 mb. 
and the influx of moisture at the surface (Chart XI of 

[3] was moderately strong from the Gulf of Mexico to 

Oklahoma. 

running and heavy precipitation in this area. 


Thus, conditions were optimum for over- 
Record 
amounts for June were observed at Waco and Abilene, 
Tex., and Shreveport, La. 

The distribution of precipitation can also be considered 
as a funetion of frontal activity. Figure 7 shows an 
axis of maximum frontal frequency from Texas to southern 
Alabama to North Carolina. The location of this average 
frontal zone corresponds well to the area of heavy precipi- 
tation, as would be expected. However, the area of ma xi 
mum frontal occurrence from North Dakota to the Lower 
Great Lakes shows a higher frequency than that in the 
South, vet the precipitation “‘fit’’ was relatively poor, 
except in New York and Pennsylvania. The flow at 
all levels west of the trough at 700 mb. contained a 
northerly component, carrying dry air into the country. 
West ineffective in pro- 
ducing rain of any consequence, a major factor in the 


Therefore, fronts in the were 


perpetuation of drought in this area. 
DROUGHT 


A significant shortage of precipitation in the last ten 
months has produced a drought of particular severity 
in portions of the States shown in figure 8A (from [5]). 
The 700-mb. mean circulation of the United 
States and Canada for the period September 1960 to 
June 1961 (fig. 8B) was prepared in order to detect any 
Anom- 


for most 


correspondence between the flow and the drought. 
defined by mean 


this 


alous precipitation is not uniquely 


700-mb. 


components of 


contours, but there are 


PERCENTAGE OF NORMAL PRECIPITATION 
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Ficure 8.—(A) Percentage of normal precipitation (from [5]), and 


(B) mean 700-mb. contours (solid) in tens of feet and depar- 
tures from normal (dotted) in 25-ft. intervals, with centers 


labeled in feet. Both for September 1960-June 1961. 


circulation that seem to augment or inhibit the production 
of precipitation. In this instance the flow itself implies 
little and, in fact, resembles grossly the normal circulation 
that consists of a ridge over the western mountains and 
a trough over eastern North America. However, the 
dotted isopleths in figure 8B representing the anomalous 
components of the 700-mb. heights clarify the situation. 
Note, first, that heights were all above normal in the 
drought area, and, second, that the direction of anomalous 
Although both 


dry weather, the latter is most frequently associated 


flow was northerly. conditions favor 
with a deficiency of precipitation on mean circulations. 

It would be an over-simplification to attempt to ex- 
plain the drought solely in terms of the anomalous 700-mb. 
flow, but there was a significant contribution. Statistical 
evidence [6| indicates that winter precipitation in North 
Dakota is negatively correlated with 700-mb. heights 
in southern Idaho and positively correlated with heights 
in the Middle Atlantic States (both 
ficients exceed 0.5). It follows that anomalous northerly 
flow components of the type portrayed in figure 8B 
should be associated with a dr: regime in the Dakotas. 


correlation coel- 
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Weather Bureau Technical News 


ISOTOPE-POWERED AUTOMATIC WEATHER STATION 


The latest development in the automatic observation of 
weather is the incorporation of a radio-isotope power plant 
into the apparatus of automatic stations. This allows 
unattended operation of the station up to the limit of 
instrument life and eliminates the difficulties met with in 
other automatic stations from lack of a continuous power 
source. The first such station was set up in August on or 
near Graham Island in the Canadian Arctic about half- 
wa\ between the Joint Canadian-United States Arctic 
Weather Stations at Resolute. 
generated by the heat spontaneously produced by radio 
active decay of pellets of strontium-90, in the form of 
strontium titanate, enclosed in a small (5-inch) shielded 


Eureka and Power is 


cylindrical capsule. The heat is transformed directly into 
a continuous flow of electricity by 60 pairs of thermo- 
couples arranged like spokes around the cylindrical heat 
The whole generator is shielded in 4.4 inches of 
The 


overall dimensions are 20 inches long by 18 inches in 


source, 


lead and covered by an outer skin of stainless steel. 


diameter. 

The generator cylinder is the bottom element in a 
larger cylindrical container 8 feet long. The top compart- 
ments of this house the electrical apparatus, recording 
equipment, radio transmitters, and the barometer. Excess 
heat from the generator serves to maintain an interior 
operating temperature of 70° F. The anemometer and 
thermometer are exposed on a mast beside the cylinder, 
the lower 5 feet of which is buried in the ground. 

Observations of wind direction (to nearest 10° heading), 
wind speed (0 to 150 kt. +1 kt., 1 and 8 minute average), 
pressure (28.00 to 32.00 in., + 0.02 in.), and temperature 

75° to +120° F., 1° F.) is transmitted by radio in 
8 bit binary digital form once every 3 hours simultaneously 
on two frequencies, 3.36 and 4.97 megacycles. The trans- 
mitter output is 250 watts on each frequency with range 
The trans- 
mission consists of call letters and data measurements plus 


up to 1500 miles depending on the frequency. 


one repeat. Provision has been made in the package for 


measuring other elements such as precipitation, humidity, 
sky cover, etc., whenever suitable sensors become available. 

The measuring and transmitting equipment are de- 
signed and fabricated to provide reliability consistent with 
the long life of the isotope, to use a minimum of electricity, 
and to produce accurate weather data in usable form. 
Since the generator sustains no wear from moving parts, 
The 


generator could power the station for more than ten years. 


long-term maintenance-free operation is possible. 


Although these stations will always be located in remote 
or completely uninhabited areas such as the present site, 
great care was taken in the development of a safe fuel form 
and adequate, reliable shielding. Strontium becomes a 
biological problem only if it is absorbed by some living 
organism. This danger can be eliminated by using an 
insoluble compound. 
material dense enough to conserve space and relatively 
This 
compound remains stable even beyond its melting point 
of 3000° F., and its rate of solubility in fresh water is so 
low it has not been measured. 


To provide at the same time a fuel 


easy to produce, strontium titanate was selected. 


In salt water its solubility is 
measured in parts per billion. The shielding within the 
generator consists of several layers of an alloy called 
Hastelloy-C. 


immersed in se: 


It would take centuries to corrode, even if 
The outside shielding of the 
generator cylinder is % ton of lead with a final covering of 


water. 


stainless steel. The whole generator is designed so sturdily 
that it could survive a plane crash or an explosion without 
releasing its fuel. 

The station was designed and built by the nuclear di- 
vision of the Martin Company, Baltimore, Md., under a 
the Atomic 
Office of Isotopes Development in cooperation with the 
U.S. Weather Bureau. 


Department of Transport of Canada and the station was 


contract with U.S. Energy Commission, 


The site was suggested by the 


installed by a joint Canadian and United States work 


party. 
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